Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,580,6340CT42.7% 128.2 / 36.1 91E277K (GAA→AAA) hsdSspecificity determinant for hsdM and hsdR
Reads supporting (aligned to +/- strand):  ref base C (10/41);  new base T (38/0);  total (50/41)
Fisher's exact test for biased strand distribution p-value = 3.63e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAA  >  NC_000913/4580492‑4580781
                                                                                                                                                |                                                                                                                                                   
cGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCt                                                                                                                                              <  2:596483/150‑1 (MQ=255)
cGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCt                                                                                                                                              <  1:321993/150‑1 (MQ=255)
 gAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTTGTTCACttttt                                                                                                                                              >  2:1257985/1‑147 (MQ=255)
 gAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCTGa                                                                                                                                            <  1:1401865/151‑1 (MQ=255)
  aGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTTGTTTACttttt                                                                                                                                              >  1:189786/1‑146 (MQ=255)
   gCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCCCtt                                                                                                                                                 >  2:1112307/1‑144 (MQ=255)
      cGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTTGTTCACTTTTTGaa                                                                                                                                           >  1:966366/1‑147 (MQ=255)
      cGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCTGAACAt                                                                                                                                        <  1:379674/150‑1 (MQ=255)
          tAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTTGTTTACTTTTTTAAAAtttt                                                                                                                                     >  1:182651/1‑147 (MQ=255)
          tAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTTGTTTACCTTCTTAACAtttt                                                                                                                                     >  1:250660/1‑147 (MQ=255)
           aaGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGa                                                                                                                                  <  1:8894/151‑1 (MQ=255)
             gTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTTGTTTACTTTTTTTACATTCTAGaa                                                                                                                                 >  1:1401466/1‑150 (MQ=255)
             gTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGaaa                                                                                                                                <  1:394541/151‑1 (MQ=255)
             gTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTGTCTGAACATTCTAGaaa                                                                                                                                <  2:1151295/151‑1 (MQ=255)
                tATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTATATTACTTTCTTAAAATTTTAaaaaac                                                                                                                              >  1:596383/1‑148 (MQ=255)
                tATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAAcc                                                                                                                              <  1:1194469/150‑1 (MQ=255)
                 aTCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTTGTTTACTTTCTGAACATTCTAGaaa                                                                                                                                >  1:531550/1‑147 (MQ=255)
                    aGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTATTTTACTTTTTTAAAAATCTAAAAAACGaa                                                                                                                           >  2:1100282/1‑149 (MQ=255)
                         ataGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTTTTTCACTTTTTTTACATTTTTGaaaaccaaaattggt                                                                                                                    >  2:422005/1‑139 (MQ=255)
                            gCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTTTTTTACTTTTTGAAAATTTTAAAAAACGAAAATAGttttt                                                                                                                 >  1:326244/1‑150 (MQ=255)
                               aaTTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTATTTTACTTTCTTAACATTCTATAAACCAAATAATGttttttt                                                                                                               >  2:310679/1‑147 (MQ=255)
                               aaTTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTTACTTTTTTTAAATTTTAAAAAACGAAAATCGtttt                                                                                                                  >  1:893441/1‑147 (MQ=255)
                                 ttttGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTTTTTAAATTCTAGAAAACGAAAATCGTTTTGAttt                                                                                                             >  1:795244/1‑148 (MQ=255)
                                     gATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACAt                                                                                                         <  2:1115139/150‑1 (MQ=255)
                                     gATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATg                                                                                                        <  1:956642/151‑1 (MQ=255)
                                        gttgtAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGCT‑TTATTTAACTTTTCGAACA‑TCTCAAAACCGCAAATCTGTTTTTTCTCCATgg                                                                                                       >  1:722035/1‑149 (MQ=255)
                                        gttgtAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGcc                                                                                                     <  1:133009/151‑1 (MQ=255)
                                             aaTTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTTGTTTTCTTTTTTAAACTTCTAAAAAAAGAATATTGTTTTTATTTATATTGCCAg                                                                                                   >  2:314096/1‑148 (MQ=255)
                                             aaTTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCAc                                                                                                <  1:1102428/151‑1 (MQ=255)
                                             aaTTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCAc                                                                                                <  2:590155/151‑1 (MQ=255)
                                               ttttttCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTTACTTTTTGAACATTCTAGAAAACAAATATCGGTTTTATTTACATTGgcagca                                                                                                 >  2:104649/1‑148 (MQ=255)
                                                tttttCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTTGTTTACTTTTTGTACATTCTAGAAACCGAATATCGTTTTTTTCTACATGGCCAGCAc                                                                                                >  1:808074/1‑148 (MQ=255)
                                                 ttttCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTa                                                                                             <  2:1449805/150‑1 (MQ=255)
                                                   ttCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACa                                                                                           <  1:890429/150‑1 (MQ=255)
                                                      aaTAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAAc                                                                                       <  1:928513/151‑1 (MQ=255)
                                                       aTAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTTGTTTACTTTTTGAAAATTCTAGAAAACGAATATAGTTTTTATTTACATGGCCAGCCCGGACaaaa                                                                                        >  1:535376/1‑149 (MQ=255)
                                                        tAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTTACTTTTTTAACATTCTAGAAAACGAATATTGTTTTTATTTACATTGCCAGCACGTACAGa                                                                                         >  2:524442/1‑147 (MQ=255)
                                                         aaCCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTTACTTTCTTAAAATTCTTGAAAACCAAAATCGTTTTTTTTTTCATTGCCCGCACGCACAGa                                                                                         >  2:148556/1‑146 (MQ=255)
                                                          aCCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTTTTTTACTTTCTTAACATTCTTGAAACCGAATATCGTTTTTTTCTACATTGCCAGCACGTACAGAACTaa                                                                                    >  1:588960/1‑150 (MQ=255)
                                                          aCCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAAt                                                                                   <  1:1115096/151‑1 (MQ=255)
                                                            ccACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTTTTTTCCTTTCTTCAAATTTTATAAAACGAATATCGTTTTTTTTTACATTGCCAGCACGCACAGAACTaa                                                                                    >  1:1493546/1‑148 (MQ=255)
                                                            ccACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTTTTTTACTTTCTTAAAATTCTTCAAAACCAATATCGTTTTTTTTTTAATTGCCAGCACGTAAAGAACTa                                                                                     >  2:201320/1‑147 (MQ=255)
                                                            ccACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTTGTTTACTTTTTTAACATTTTAAAAAACCAAAATTGTTTTTATTTACATTGCCAGCACggaaaaaaataaatc                                                                                 >  1:1488749/1‑137 (MQ=255)
                                                            ccACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTTGTTTACTTTTTTAAAATTTTAGAAAACGAAAAATGTTTTTATTTACATGGCCAGCAAGTACAGaaaaaa                                                                                    >  2:1528839/1‑144 (MQ=255)
                                                            ccACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTTACTTTCTTTAAATTCTATAAAACAAATATTGTTTTTTTTTAAATAGCCGGCACCCacaaaacaaaac                                                                                  >  2:459860/1‑148 (MQ=255)
                                                             cacaAACACCAACAAATTCTAAACTTCCGTTATAACGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCg                                                                                <  1:74875/151‑1 (MQ=255)
                                                                      caacaaATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTTCTTTTTTTTTTCAACATTGCATAAAAACAATATTGTTTTGATTTACAATGCCCGCACGTTCAGCACTaacgagttgcatt                                                                         >  2:9571/1‑138 (MQ=255)
                                                                        acaaATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTTGTTTACTTTCTTAACATTCTTGAAAACGAATATCGTTTTTATCTACATGGCCAGCACGTACAAAAATAATGCGTAGGAtttggt                                                                     >  1:504807/1‑147 (MQ=255)
                                                                        acaaATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGAt                                                                     <  1:1316748/151‑1 (MQ=255)
                                                                         caaATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTTGTTTACTTTCTTTACATTTTTGAAACCGGAAATCGTTTTTTTTTTCATGGCCCGCACATACAAAAATAATGCGTAGTAtt                                                                         >  1:591405/1‑146 (MQ=255)
                                                                          aaaTTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTTGTTCACTTTTTGTACATTTTTGAAAACCAATATTGTTTTTATTTACATGGGCAGCACGGACAGAAATAAATGGTAGTAtttg                                                                       >  1:1130094/1‑145 (MQ=255)
                                                                          aaaTTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGa                                                                   <  1:775970/151‑1 (MQ=255)
                                                                          aaaTTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGa                                                                   <  2:344331/151‑1 (MQ=255)
                                                                           aaTTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGAc                                                                  <  1:487430/151‑1 (MQ=255)
                                                                              tCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCa                                                                <  2:1407945/150‑1 (MQ=255)
                                                                               cTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTTGTTTTCTTTTTTAAAATTCTATAAACCGAATATTGTTTTTATTTACATTGCCAGCACGTACAAAACTAATACGTAGAAtttg                                                                       >  2:114603/1‑140 (MQ=255)
                                                                               cTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTTACTTTTTTAAAATTCTATAAAACGAATATAGTTTTTATCTACATGGCCAGCACCTACAAAAATAAAGCGGAGTATTGGATGACCa                                                                >  2:882380/1‑149 (MQ=255)
                                                                                  aaCTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAAcacc                                                           <  2:1044539/151‑1 (MQ=255)
                                                                                      tCCGTTATTGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAAcaccaat                                                        <  2:163189/150‑3 (MQ=255)
                                                                                         gTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTTGTTTACTTTTTTAACATTCTAGAAAACGAATATCGTTTTTATTTTCATTGCCAGCACGTACAGAAATAATGCGTAGTATTTGATGAACAAAACCACttt                                                      >  1:599949/1‑149 (MQ=255)
                                                                                         gTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTTACTTTTTGAACATTCTAGAAAAAAAATATCGTTTTGATCTACATGGCCACCCCATACAAAAATAATGCGTAGTTTTTGATGACCaaaaacccctt                                                      >  2:1482638/1‑140 (MQ=255)
                                                                                         gTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCa                                                    <  2:1506392/151‑1 (MQ=255)
                                                                                         gTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCa                                                    <  1:215377/151‑1 (MQ=255)
                                                                                           tataGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTTGTTTAATTTTTTAACATTCTAAAAAACAAAAATAGTTTTTTTTTACATTGCCAACACGTACAAAAACAAATCGTAGTATTTGGTGACCAACACCACttttttt                                                  >  2:1168981/1‑147 (MQ=255)
                                                                                              aGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTg                                                <  1:361086/150‑1 (MQ=255)
                                                                                                 gAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTTGTTTACTTTTTTAAAAATCTAGAAACCGAATATTGTTTTTTTCTTCATGGGCCGCACGGACAAAAATAATACGTAGTATTTGATGACCAACACCACTTTCATTTTGCtt                                            >  1:553073/1‑151 (MQ=255)
                                                                                                 gAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTATTTTAATTTTTGGACATTTTAGAAAACGAAAATACTTTTTTTCTTCATGGCCCGCACGTACAGAACTAATGCGGAGTATTTGAGTATCAAAACCACATTTCTTTGGcgt                                            >  1:294767/1‑149 (MQ=255)
                                                                                                  aGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCttt                                           <  2:707162/151‑1 (MQ=255)
                                                                                                   gTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTTTTTTACTTTTTTAAAATTCTTTAAAACCAAAATCGTTTTTTTTTTCATGGCCAGCCCGTACAAAACTAATGCGTAGTATTTGATGAACAAAACCACTTTCATTTTGCt                                             >  2:514287/1‑148 (MQ=255)
                                                                                                     aaataaaAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTTGCTTCCATTCTTAAACTTTTATAAAACACATAACGTTTTTTTTTTCATAGCCCGGACGTACAGAACTAATGCGTAGGATTTGTTGACCAAAAACCCCTTCTTTTGGtttt                                           >  2:1303147/1‑148 (MQ=255)
                                                                                                       ataaaAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTTCCTTTCTGAAAATTTTTGAAAACCAATATCGTTTTGTTTTACATGGGCAGCACGTACAGAACTAATGCGTAGGATTGGGTGACCAACACCACcttttttttggttt                                           >  1:1254441/1‑135 (MQ=255)
                                                                                                       ataaaAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACGTTCATTTGGCTTtgatga                                      <  2:1223672/151‑1 (MQ=255)
                                                                                                          aaaGATCTCCATCTTGTAATTTGTGGCGG‑G‑TTAGTTTACTTTTTTAAAATTCTTGAAAACGAATAATGGTTTTTTTTTCAAGGTCATCCCGTACCAAACTAAATCGTAGTATTTTATGAACCAAACCCCCTTCTTTTGGGTTTTATGgaa                                    >  2:142303/1‑150 (MQ=255)
                                                                                                            aGATCTCCATCTTGTAATTTGTGGCGG‑TATTTGTT‑TCTTTCTGTAACTTTTAAAAAACCAATATCGTTTTTGTTTAAATGGCGAGCCCGTATACAAATAAAATGTAGTGTTGTATGAACAACAACACTTTTATTTTGCtttttttgaaaa                                  >  2:1477774/1‑141 (MQ=255)
                                                                                                             gATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAATTTACTTTCTTAACATTTTAGAAAACGAATATAGTTTTGATCTACAAGGCCAACACCTACAGAACTAATGCGGAGTATTGGGTGACCCACACCCCCTTTATTTTGCTTTGATGGAAAAc                                 >  1:967401/1‑150 (MQ=255)
                                                                                                                  ccATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTAc                           <  1:887609/151‑1 (MQ=255)
                                                                                                                        tGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGt                         <  1:1461913/147‑1 (MQ=255)
                                                                                                                         gTAATTTGTGGCGG‑T‑TTAGTTTACTTTTTGAACATTCTAGAAAACGAATATCGTTTTGTTTTTCATTGCCCGCACCTACAGAAATAATTCGTAAGATTTGATTACCCACACCCCTTTTATTTTGGTTTGATGAAAAAACAATACGTaaat                     >  2:202841/1‑148 (MQ=255)
                                                                                                                                gtgGCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTaa              <  2:323012/150‑1 (MQ=255)
                                                                                                                                gtgGCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAag             <  1:1304569/151‑1 (MQ=255)
                                                                                                                                  ggCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAagat           <  2:649815/151‑1 (MQ=255)
                                                                                                                                   gCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAagata          <  2:41848/151‑1 (MQ=255)
                                                                                                                                   gCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAagata          <  2:1379768/151‑1 (MQ=255)
                                                                                                                                    cGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAagat           >  1:1530601/1‑149 (MQ=255)
                                                                                                                                     gg‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAagataga        <  2:943154/151‑1 (MQ=255)
                                                                                                                                     gg‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAagat           >  2:693159/1‑148 (MQ=255)
                                                                                                                                      g‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAagatagat       <  2:159223/151‑1 (MQ=255)
                                                                                                                                          ttAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTc     >  2:872171/1‑151 (MQ=255)
                                                                                                                                            aGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCa    >  2:425159/1‑150 (MQ=255)
                                                                                                                                            aGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCa    <  1:1453334/150‑1 (MQ=255)
                                                                                                                                             gTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCaaa  <  1:422022/151‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGG‑T‑TTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAA  >  NC_000913/4580492‑4580781

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑