Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,776,881 G→A 5.4% R2K (AGA→AAA)  ypjJ → uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,776,8810GA 5.4% 229.4 / 2.2 74R2K (AGA→AAA) ypjJuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (33/37);  new base A (2/2);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-01

CCCAGCCAGGCAGATAAGACTATAACGCAGCGTCTGGTGCAGGCGCTTCAGCTCGTTGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCACCGGAAAATATCAGTGGCACGGTAGCGTCTGTCATTACACCGGCAGGGATGTTCCGGATATCACAGGAGTCCTGGC  >  NC_000913/2776737‑2777028
                                                                                                                                                |                                                                                                                                                   
cccAGCCAGGCAGATAAGACTATAACGCAGCGTCTGGTGCAGGCGCTTCAGCTCGTTGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAAt                                                                                                                                                  >  8:4671/1‑148 (MQ=255)
 ccagccagGCAGATAAGACTATAACGCAGCGTCTGGTGCAGGCGCTTCAGCTCGTTGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATTAGAAt                                                                                                                                                  >  4:77220/1‑147 (MQ=255)
    gccagGCAGATAAGACTATAACGCAGCGTCTGGTGCAGGCGCTTCAGCTCGTTGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATc                                                                                                                                              <  4:155586/148‑1 (MQ=255)
    gccagGCAGATAAGACTATAACGCAGCGTCTGGTGCAGGCGCTTCAGCTCGTTGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATc                                                                                                                                              <  8:275604/148‑1 (MQ=255)
    gccagGCAGATAAGACTATAACGCAGCGTCTGGTGCAGGCGCTTCAGCTCGTTGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATc                                                                                                                                              <  6:147188/148‑1 (MQ=255)
        ggCAGATAAGACTATAACGCAGCGTCTGGTGCAGGCGCTTCAGCTCGTTGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTa                                                                                                                                          <  5:260991/148‑1 (MQ=255)
         gCAGATAAGACTATAACGCAGCGTCTGGTGCAGGCGCTTCAGCTCGTTGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTaa                                                                                                                                         <  7:216062/148‑1 (MQ=255)
         gCAGATAAGACTATAACGCAGCGTCTGGTGCAGGCGCTTCAGCTCGTTGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTaa                                                                                                                                         <  8:390046/148‑1 (MQ=255)
             aTAAGACTATAACGCAGCGTCTGGTGCAGGCGCTTCAGCGCGTTGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGc                                                                                                                                     <  2:303923/148‑1 (MQ=255)
                aGACTATAACGCAGCGTCTGGTGCAGGCGCTTGAGCTCGTTGATTTCCGTGTCGCTGACCATCTTATTGTTGGGGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCg                                                                                                                                   <  6:23458/147‑1 (MQ=17)
                 gACTATAACGCAGCGTCTGGTGCAGGCGCTTCAGCTCGTTGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCggg                                                                                                                                 >  2:4357/1‑148 (MQ=255)
                  aCTATAACGCAGCGTCTGGTGCAGGCGCTTCAGCTCGTTGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCggg                                                                                                                                 >  3:257491/1‑147 (MQ=255)
                  aCTATAACGCAGCGTCTGGTGCAGGCGCTTCAGCTCGTTGATATCCGTGTCCCTGAACATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCAggg                                                                                                                                 >  3:365070/1‑147 (MQ=255)
                   cTATAACGCAGCGTCTGGTGCAGGCGCTTCAGCTCGTTGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCa                                                                                                                               <  2:324221/148‑1 (MQ=255)
                    tataACGCAGCGTCTGGTGCAGGCGCTTCAGCTCGTTGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAAAATTATCAGTAAACGCCGGGCaa                                                                                                                              >  7:239404/1‑148 (MQ=255)
                     ataACGCAGCGTCTGGTGCAGGCGCTTCAGCTCGTTGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCa                                                                                                                                             <  3:41462/132‑1 (MQ=255)
                     ataACGCAGCGTCTGGTGCAGGCGCTTCAGCTCGTTGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCa                                                                                                                                             >  4:41462/1‑132 (MQ=255)
                     ataACGCAGCGTCTGGTGCAGGCGCTTCAGCTCGTTGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAAt                                                                                                                             >  8:378979/1‑148 (MQ=255)
                     ataACGCAGCGTCTGGTGCAGGCGCTTCAGCTCGTTGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAAt                                                                                                                             >  8:149328/1‑148 (MQ=255)
                      taACGCAGCGTCTGGTGCAGGCGCTTCAGCTCGTTGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAAt                                                                                                                             >  7:51582/1‑147 (MQ=255)
                      taACGCAGCGTCTGGTGCAGGCGCTTCAGCTCGTTGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATg                                                                                                                            <  7:378979/148‑1 (MQ=255)
                        aCGCAGCGTCTGGTGCAGGCGCTTCAGCTCGGTGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGa                                                                                                                           <  4:221401/147‑1 (MQ=255)
                         cGCAGCGTCTGGTGCAGGCGCTTCAGCTCGTTGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAAGGACg                                                                                                                         <  7:365734/148‑1 (MQ=255)
                           cAGCGTCTGGTGCAGGCGCTTCAGCTCGTTGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGat                                                                                                                       <  3:91037/148‑1 (MQ=255)
                           cAGCGTCTGGTGCAGGCGCTTCAGCTCGTTGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGat                                                                                                                       <  7:149328/148‑1 (MQ=255)
                           cAGCGTCTGGTGCAGGCGCTTAAGCTAGTTGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGATCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGat                                                                                                                       <  1:142157/148‑1 (MQ=255)
                                 cTGGTGCAGGCGCTTCAGCTCGTTGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAAAATTa                                                                                                                                                >  7:288164/1‑117 (MQ=255)
                                 cTGGTGCAGGCGCTTCAGCTCGTTGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAAAATTa                                                                                                                                                <  8:288164/117‑1 (MQ=255)
                                        atgcgcTTCAGCTCGTTGTTATCCGTGTGCATGACAATCTGATTGTTGTTGGCAGTCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGTATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCAt                                                                                                          <  4:365070/146‑1 (MQ=14)
                                           cgcTTCAGCTCGTTGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCt                                                                                                       >  2:173615/1‑148 (MQ=255)
                                            gcTTCAGCTCGTTGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAAAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTg                                                                                                      <  8:239404/148‑1 (MQ=37)
                                                      gTTGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGaa                                                                                             <  3:181793/147‑1 (MQ=39)
                                                        tGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATgaga                                                                                                                                                    >  2:194722/1‑90 (MQ=38)
                                                        tGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATgaga                                                                                                                                                    <  1:194722/90‑1 (MQ=38)
                                                        tGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCt                                                                                          >  5:279326/1‑148 (MQ=38)
                                                               cttgtCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTa                                                                                   <  7:4671/146‑1 (MQ=37)
                                                                      cTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCAcc                                                                            >  3:21490/1‑148 (MQ=37)
                                                                      cTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCAcc                                                                            >  6:224076/1‑148 (MQ=37)
                                                                       tGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGAAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCttt                                                                                        >  7:122817/1‑135 (MQ=37)
                                                                       tGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGAAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCttt                                                                                        <  8:122817/135‑1 (MQ=37)
                                                                        gACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGc                                                                                                            <  3:377742/114‑1 (MQ=39)
                                                                        gACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGc                                                                                                            >  4:377742/1‑114 (MQ=39)
                                                                        gACAATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGTTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATAACGCCAGCATCCTGATTCCCGAATCTTTCGCTAATGCACCg                                                                           >  2:170756/1‑147 (MQ=11)
                                                                            aTCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCACCGGaaaa                                                                      >  4:33503/1‑148 (MQ=39)
                                                                             tCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGa                                                                                                                           >  4:172413/1‑94 (MQ=40)
                                                                             tCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGa                                                                                                                           <  3:172413/94‑1 (MQ=40)
                                                                              cTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCACCGGAAAata                                                                    >  5:63100/1‑148 (MQ=39)
                                                                                 aTTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCACCGGAAAat                                                                     <  3:77220/144‑1 (MQ=39)
                                                                                 aTTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCACCGGAAAATATc                                                                  <  1:27882/147‑1 (MQ=40)
                                                                                    gttgGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCACCGGAAAATATCAGt                                                               <  4:257491/147‑1 (MQ=255)
                                                                                             aggcaAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCACCGGAAAATATCAGTGGCACGGTAg                                                     <  1:173615/148‑1 (MQ=39)
                                                                                              ggcaAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCACCGGAAAATATCAGTGGCACGGTAGc                                                    >  4:173663/1‑148 (MQ=39)
                                                                                              ggcaAATCTATTCGTTCGCAGAACACGGTCTGCGTTGAGGTATTACATGAGAAGTATCAGTAAACGCCGGGCAATGCCGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCGGCGGCAAATATCAGTGGCACGGTAGc                                                    >  7:305586/1‑148 (MQ=11)
                                                                                                  aaTGTATTCGTTCGCAGAACAGGGTCTGCGTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCAAGCTGAGTCCCGAATCATTCGCTACTGCACCGGAAAATATAAGTGGCACGGTAGCGTCt                                                >  8:281541/1‑148 (MQ=14)
                                                                                                     cTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGGGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCACCGGAAAATATCAGTGGCACGGTAGCGTCTGt                                              >  4:102942/1‑147 (MQ=255)
                                                                                                        ttcgttcgCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCACCGGAAAATATCAGTGGCACGGTAGCGTCTGTCAtt                                          <  1:4357/148‑1 (MQ=39)
                                                                                                         tcgttcgCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCACCGGAAAATATCAGTGGCACGGTAGCGTCTGTCATTa                                         >  8:310517/1‑148 (MQ=39)
                                                                                                            ttcgCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCACCGGAAAATATCCGTGGCACGGTAGCGTCTGTCATTacac                                      <  4:299160/148‑1 (MQ=255)
                                                                                                                 aGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCACCGGAAAATATCAGTGGCACGGTAGCGTCTGTCATTACACCGGCa                                 <  3:33503/148‑1 (MQ=39)
                                                                                                                 aGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCACCGGAAAATATCAGTGGCACGGTAGCGTCTGTCATTACACCGGCa                                 <  6:279326/148‑1 (MQ=39)
                                                                                                                   aaCACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCACCGGAAAATATCAGTGGCACGGTAGCGTCTGTCATTACACCGGCAgg                               >  5:189244/1‑148 (MQ=39)
                                                                                                                   aaCACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCACCGGAAAATATCAGTGGCACGGTAGCGTCTGTCATTACACCGGCAgg                               <  3:173663/148‑1 (MQ=39)
                                                                                                                    acacGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCACCGGAAAATATCAGTGACACGGTAGCGTCTGTCATTACACCGGCAggg                              <  6:63100/148‑1 (MQ=255)
                                                                                                                     cacGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCACCGGAAAATATCAGTGGCACGGTAGCGTCTGTCATTACACCGGCAGGGa                             <  3:232060/148‑1 (MQ=39)
                                                                                                                             gCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCACCGGAAAATATCAGTGGCACGGTAGCGTCTGTCATTACACCGGCAGGGATGTTCCg                      >  2:18793/1‑147 (MQ=255)
                                                                                                                                  gAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCt                                                                                          <  6:141334/74‑1 (MQ=39)
                                                                                                                                  gAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCt                                                                                          >  5:141334/1‑74 (MQ=39)
                                                                                                                                      tATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCCGAGTCCCGAATCTTTCGCTACTGCACCGGAAAATATCAGTGGCACGGTAGCGTCTGTCATTACACCGGCAGGGATGTTCCGGATATCACAg            >  4:11573/1‑148 (MQ=255)
                                                                                                                                       aTTACATGAGAATTATCAGTAAACGCCGGGCAATGCCGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCACCGGAAAATATCAGTGGCACGGTAGCGTCTGTCATTACACCGGCAGGGATGTTCCGGATATCACAg            >  5:364577/1‑147 (MQ=255)
                                                                                                                                       aTTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCACCGGAAAATATCAGTGGCACGGTAGCGTCTGTCATTACACCGGCAGGGATGTTCCGGATATCACAg            >  1:148679/1‑147 (MQ=255)
                                                                                                                                         tACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCg                                                                                      <  7:339532/71‑1 (MQ=39)
                                                                                                                                         tACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCg                                                                                      >  8:339532/1‑71 (MQ=39)
                                                                                                                                         tACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCACCGGAAAATATCAGTGGCACGGTAGCGTCTGTCATTACACCGGCAGGGATGTTCCGGATATAACAGGAg         >  5:286956/1‑148 (MQ=255)
                                                                                                                                            aTGGGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCACCGGAAAATATCAGTGGCACGGTAGCGTCTGTCATTACACCGGCAGGGATGTTCCGGATATCACAGGAGTcc      <  3:102942/148‑1 (MQ=255)
                                                                                                                                            aTGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCACCGGAAAATATCAGTGGCACGGTAGCGTCTGTCATTACACCGGCAGGGATGTTCCGGATATCACAGGAGTcc      <  7:211883/148‑1 (MQ=255)
                                                                                                                                            aTGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCACCGGAAAATATCAGTGGCACGGTAGCGTCTGTCATTACACCGGCAGGGATGTTCCGGATATCACAGGAGTcc      >  2:55950/1‑148 (MQ=255)
                                                                                                                                              gagaATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCACCGGAAAATATc                                                                  <  2:322537/86‑1 (MQ=255)
                                                                                                                                              gagaATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCACCGGAAAATATc                                                                  >  1:322537/1‑86 (MQ=255)
                                                                                                                                              gagaATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCACCGGAAAATATCAGTGGCACGGTAGCGTCTGTCATTACACCGGCAGGGATGTTCCGGATATc                >  7:253410/1‑136 (MQ=255)
                                                                                                                                              gagaATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCACCGGAAAATATCAGTGGCACGGTAGCGTCTGTCATTACACCGGCAGGGATGTTCCGGATATc                <  8:253410/136‑1 (MQ=255)
                                                                                                                                               agaATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCACCGGAAAATATCAGTGGCACGGTAGCGTCTGTCATTACACCGGCAGGGATGTTCCGGATATCACAGGAGTCCTgg   >  8:51240/1‑148 (MQ=255)
                                                                                                                                                gaATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCACCGGAAAATATCAGTGGCACGGTAGCGTCTGTCATTACACCGGCAGGGATGTTCCGGATATCACAGGAGTCCTGGc  >  5:20967/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CCCAGCCAGGCAGATAAGACTATAACGCAGCGTCTGGTGCAGGCGCTTCAGCTCGTTGATATCCGTGTCCCTGACCATCTGATTGTTGGTGGCAGGCAAATCTATTCGTTCGCAGAACACGGTCTGCTTTGAGGTATTACATGAGAATTATCAGTAAACGCCGGGCAATGACGATATACCGCCAGCATCCTGAGTCCCGAATCTTTCGCTACTGCACCGGAAAATATCAGTGGCACGGTAGCGTCTGTCATTACACCGGCAGGGATGTTCCGGATATCACAGGAGTCCTGGC  >  NC_000913/2776737‑2777028

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: