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

CGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAA  >  NC_000913/4580492‑4580781
                                                                                                                                                |                                                                                                                                                   
cGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCt                                                                                                                                              <  2:596483/150‑1 (MQ=255)
cGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCt                                                                                                                                              <  1:321993/150‑1 (MQ=255)
 gAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTCACttttt                                                                                                                                              >  2:1257985/1‑147 (MQ=255)
 gAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGa                                                                                                                                            <  1:1401865/151‑1 (MQ=255)
  aGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTTACttttt                                                                                                                                              >  1:189786/1‑146 (MQ=255)
   gCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCCCtt                                                                                                                                                 >  2:1112307/1‑144 (MQ=255)
      cGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTCACTTTTTGaa                                                                                                                                           >  1:966366/1‑147 (MQ=255)
      cGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACAt                                                                                                                                        <  1:379674/150‑1 (MQ=255)
          tAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTTACTTTTTTAAAAtttt                                                                                                                                     >  1:182651/1‑147 (MQ=255)
          tAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTTACCTTCTTAACAtttt                                                                                                                                     >  1:250660/1‑147 (MQ=255)
           aaGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGa                                                                                                                                  <  1:8894/151‑1 (MQ=255)
             gTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTTACTTTTTTTACATTCTAGaa                                                                                                                                 >  1:1401466/1‑150 (MQ=255)
             gTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGaaa                                                                                                                                <  1:394541/151‑1 (MQ=255)
             gTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTGTCTGAACATTCTAGaaa                                                                                                                                <  2:1151295/151‑1 (MQ=255)
                tATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTATATTACTTTCTTAAAATTTTAaaaaac                                                                                                                              >  1:596383/1‑148 (MQ=255)
                tATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAAcc                                                                                                                              <  1:1194469/150‑1 (MQ=255)
                 aTCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTTACTTTCTGAACATTCTAGaaa                                                                                                                                >  1:531550/1‑147 (MQ=255)
                    aGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTATTTTACTTTTTTAAAAATCTAAAAAACGaa                                                                                                                           >  2:1100282/1‑149 (MQ=255)
                         ataGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTTTTCACTTTTTTTACATTTTTGaaaaccaaaattggt                                                                                                                    >  2:422005/1‑139 (MQ=255)
                            gCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTTTTTACTTTTTGAAAATTTTAAAAAACGAAAATAGttttt                                                                                                                 >  1:326244/1‑150 (MQ=255)
                               aaTTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTATTTTACTTTCTTAACATTCTATAAACCAAATAATGttttttt                                                                                                               >  2:310679/1‑147 (MQ=255)
                               aaTTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTTACTTTTTTTAAATTTTAAAAAACGAAAATCGtttt                                                                                                                  >  1:893441/1‑147 (MQ=255)
                                 ttttGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTTTTTAAATTCTAGAAAACGAAAATCGTTTTGAttt                                                                                                             >  1:795244/1‑148 (MQ=255)
                                     gATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACAt                                                                                                         <  2:1115139/150‑1 (MQ=255)
                                     gATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATg                                                                                                        <  1:956642/151‑1 (MQ=255)
                                        gttgtAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGCTTTATTTAACTTTTCGAACATCTCAAAACCGCAAATCTGTTTTTTCTCCATgg                                                                                                       >  1:722035/1‑149 (MQ=255)
                                        gttgtAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGcc                                                                                                     <  1:133009/151‑1 (MQ=255)
                                             aaTTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTTTCTTTTTTAAACTTCTAAAAAAAGAATATTGTTTTTATTTATATTGCCAg                                                                                                   >  2:314096/1‑148 (MQ=255)
                                             aaTTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCAc                                                                                                <  1:1102428/151‑1 (MQ=255)
                                             aaTTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCAc                                                                                                <  2:590155/151‑1 (MQ=255)
                                               ttttttCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTTACTTTTTGAACATTCTAGAAAACAAATATCGGTTTTATTTACATTGgcagca                                                                                                 >  2:104649/1‑148 (MQ=255)
                                                tttttCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTTACTTTTTGTACATTCTAGAAACCGAATATCGTTTTTTTCTACATGGCCAGCAc                                                                                                >  1:808074/1‑148 (MQ=255)
                                                 ttttCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTa                                                                                             <  2:1449805/150‑1 (MQ=255)
                                                   ttCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACa                                                                                           <  1:890429/150‑1 (MQ=255)
                                                      aaTAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAAc                                                                                       <  1:928513/151‑1 (MQ=255)
                                                       aTAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTTACTTTTTGAAAATTCTAGAAAACGAATATAGTTTTTATTTACATGGCCAGCCCGGACaaaa                                                                                        >  1:535376/1‑149 (MQ=255)
                                                        tAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTTACTTTTTTAACATTCTAGAAAACGAATATTGTTTTTATTTACATTGCCAGCACGTACAGa                                                                                         >  2:524442/1‑147 (MQ=255)
                                                         aaCCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTTACTTTCTTAAAATTCTTGAAAACCAAAATCGTTTTTTTTTTCATTGCCCGCACGCACAGa                                                                                         >  2:148556/1‑146 (MQ=255)
                                                          aCCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTTTTTACTTTCTTAACATTCTTGAAACCGAATATCGTTTTTTTCTACATTGCCAGCACGTACAGAACTaa                                                                                    >  1:588960/1‑150 (MQ=255)
                                                          aCCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAAt                                                                                   <  1:1115096/151‑1 (MQ=255)
                                                            ccACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTTTTTCCTTTCTTCAAATTTTATAAAACGAATATCGTTTTTTTTTACATTGCCAGCACGCACAGAACTaa                                                                                    >  1:1493546/1‑148 (MQ=255)
                                                            ccACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTTTTTACTTTCTTAAAATTCTTCAAAACCAATATCGTTTTTTTTTTAATTGCCAGCACGTAAAGAACTa                                                                                     >  2:201320/1‑147 (MQ=255)
                                                            ccACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTTACTTTTTTAACATTTTAAAAAACCAAAATTGTTTTTATTTACATTGCCAGCACggaaaaaaataaatc                                                                                 >  1:1488749/1‑137 (MQ=255)
                                                            ccACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTTACTTTTTTAAAATTTTAGAAAACGAAAAATGTTTTTATTTACATGGCCAGCAAGTACAGaaaaaa                                                                                    >  2:1528839/1‑144 (MQ=255)
                                                            ccACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTTACTTTCTTTAAATTCTATAAAACAAATATTGTTTTTTTTTAAATAGCCGGCACCCacaaaacaaaac                                                                                  >  2:459860/1‑148 (MQ=255)
                                                             cacaAACACCAACAAATTCTAAACTTCCGTTATAACGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCg                                                                                <  1:74875/151‑1 (MQ=255)
                                                                      caacaaATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTCTTTTTTTTTTCAACATTGCATAAAAACAATATTGTTTTGATTTACAATGCCCGCACGTTCAGCACTaacgagttgcatt                                                                         >  2:9571/1‑138 (MQ=255)
                                                                        acaaATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTTACTTTCTTAACATTCTTGAAAACGAATATCGTTTTTATCTACATGGCCAGCACGTACAAAAATAATGCGTAGGAtttggt                                                                     >  1:504807/1‑147 (MQ=255)
                                                                        acaaATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGAt                                                                     <  1:1316748/151‑1 (MQ=255)
                                                                         caaATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTTACTTTCTTTACATTTTTGAAACCGGAAATCGTTTTTTTTTTCATGGCCCGCACATACAAAAATAATGCGTAGTAtt                                                                         >  1:591405/1‑146 (MQ=255)
                                                                          aaaTTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTCACTTTTTGTACATTTTTGAAAACCAATATTGTTTTTATTTACATGGGCAGCACGGACAGAAATAAATGGTAGTAtttg                                                                       >  1:1130094/1‑145 (MQ=255)
                                                                          aaaTTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGa                                                                   <  1:775970/151‑1 (MQ=255)
                                                                          aaaTTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGa                                                                   <  2:344331/151‑1 (MQ=255)
                                                                           aaTTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGAc                                                                  <  1:487430/151‑1 (MQ=255)
                                                                              tCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCa                                                                <  2:1407945/150‑1 (MQ=255)
                                                                               cTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTTTCTTTTTTAAAATTCTATAAACCGAATATTGTTTTTATTTACATTGCCAGCACGTACAAAACTAATACGTAGAAtttg                                                                       >  2:114603/1‑140 (MQ=255)
                                                                               cTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTTACTTTTTTAAAATTCTATAAAACGAATATAGTTTTTATCTACATGGCCAGCACCTACAAAAATAAAGCGGAGTATTGGATGACCa                                                                >  2:882380/1‑149 (MQ=255)
                                                                                  aaCTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAAcacc                                                           <  2:1044539/151‑1 (MQ=255)
                                                                                      tCCGTTATTGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAAcaccaat                                                        <  2:163189/150‑3 (MQ=255)
                                                                                         gTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTTACTTTTTTAACATTCTAGAAAACGAATATCGTTTTTATTTTCATTGCCAGCACGTACAGAAATAATGCGTAGTATTTGATGAACAAAACCACttt                                                      >  1:599949/1‑149 (MQ=255)
                                                                                         gTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTTACTTTTTGAACATTCTAGAAAAAAAATATCGTTTTGATCTACATGGCCACCCCATACAAAAATAATGCGTAGTTTTTGATGACCaaaaacccctt                                                      >  2:1482638/1‑140 (MQ=255)
                                                                                         gTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCa                                                    <  2:1506392/151‑1 (MQ=255)
                                                                                         gTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCa                                                    <  1:215377/151‑1 (MQ=255)
                                                                                           tataGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTTAATTTTTTAACATTCTAAAAAACAAAAATAGTTTTTTTTTACATTGCCAACACGTACAAAAACAAATCGTAGTATTTGGTGACCAACACCACttttttt                                                  >  2:1168981/1‑147 (MQ=255)
                                                                                              aGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTg                                                <  1:361086/150‑1 (MQ=255)
                                                                                                 gAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTGTTTACTTTTTTAAAAATCTAGAAACCGAATATTGTTTTTTTCTTCATGGGCCGCACGGACAAAAATAATACGTAGTATTTGATGACCAACACCACTTTCATTTTGCtt                                            >  1:553073/1‑151 (MQ=255)
                                                                                                 gAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTATTTTAATTTTTGGACATTTTAGAAAACGAAAATACTTTTTTTCTTCATGGCCCGCACGTACAGAACTAATGCGGAGTATTTGAGTATCAAAACCACATTTCTTTGGcgt                                            >  1:294767/1‑149 (MQ=255)
                                                                                                  aGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCttt                                           <  2:707162/151‑1 (MQ=255)
                                                                                                   gTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTTTTTTACTTTTTTAAAATTCTTTAAAACCAAAATCGTTTTTTTTTTCATGGCCAGCCCGTACAAAACTAATGCGTAGTATTTGATGAACAAAACCACTTTCATTTTGCt                                             >  2:514287/1‑148 (MQ=255)
                                                                                                     aaataaaAGATCTCCATCTTGTAATTTGTGGCGGTTTTGCTTCCATTCTTAAACTTTTATAAAACACATAACGTTTTTTTTTTCATAGCCCGGACGTACAGAACTAATGCGTAGGATTTGTTGACCAAAAACCCCTTCTTTTGGtttt                                           >  2:1303147/1‑148 (MQ=255)
                                                                                                       ataaaAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTTCCTTTCTGAAAATTTTTGAAAACCAATATCGTTTTGTTTTACATGGGCAGCACGTACAGAACTAATGCGTAGGATTGGGTGACCAACACCACcttttttttggttt                                           >  1:1254441/1‑135 (MQ=255)
                                                                                                       ataaaAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACGTTCATTTGGCTTtgatga                                      <  2:1223672/151‑1 (MQ=255)
                                                                                                          aaaGATCTCCATCTTGTAATTTGTGGCGGGTTAGTTTACTTTTTTAAAATTCTTGAAAACGAATAATGGTTTTTTTTTCAAGGTCATCCCGTACCAAACTAAATCGTAGTATTTTATGAACCAAACCCCCTTCTTTTGGGTTTTATGgaa                                    >  2:142303/1‑150 (MQ=255)
                                                                                                            aGATCTCCATCTTGTAATTTGTGGCGGTATTTGTTTCTTTCTGTAACTTTTAAAAAACCAATATCGTTTTTGTTTAAATGGCGAGCCCGTATACAAATAAAATGTAGTGTTGTATGAACAACAACACTTTTATTTTGCtttttttgaaaa                                  >  2:1477774/1‑141 (MQ=255)
                                                                                                             gATCTCCATCTTGTAATTTGTGGCGGTTTAATTTACTTTCTTAACATTTTAGAAAACGAATATAGTTTTGATCTACAAGGCCAACACCTACAGAACTAATGCGGAGTATTGGGTGACCCACACCCCCTTTATTTTGCTTTGATGGAAAAc                                 >  1:967401/1‑150 (MQ=255)
                                                                                                                  ccATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTAc                           <  1:887609/151‑1 (MQ=255)
                                                                                                                        tGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGt                         <  1:1461913/147‑1 (MQ=255)
                                                                                                                         gTAATTTGTGGCGGTTTAGTTTACTTTTTGAACATTCTAGAAAACGAATATCGTTTTGTTTTTCATTGCCCGCACCTACAGAAATAATTCGTAAGATTTGATTACCCACACCCCTTTTATTTTGGTTTGATGAAAAAACAATACGTaaat                     >  2:202841/1‑148 (MQ=255)
                                                                                                                                gtgGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTaa              <  2:323012/150‑1 (MQ=255)
                                                                                                                                gtgGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAag             <  1:1304569/151‑1 (MQ=255)
                                                                                                                                  ggCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAagat           <  2:649815/151‑1 (MQ=255)
                                                                                                                                   gCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAagata          <  2:41848/151‑1 (MQ=255)
                                                                                                                                   gCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAagata          <  2:1379768/151‑1 (MQ=255)
                                                                                                                                    cGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAagat           >  1:1530601/1‑149 (MQ=255)
                                                                                                                                     ggTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAagataga        <  2:943154/151‑1 (MQ=255)
                                                                                                                                     ggTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAagat           >  2:693159/1‑148 (MQ=255)
                                                                                                                                      gTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAagatagat       <  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)
                                                                                                                                                |                                                                                                                                                   
CGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGATGTTGTAATTTTTTCAATAACCCACAAACACCAACAAATTCTAAACTTCCGTTATAGCGAGTAAATAAAAGATCTCCATCTTGTAATTTGTGGCGGTTTAGTTCACTTTCTGAACATTCTAGAAACCGAATATCGTTTTGATCTACATGGCCAGCACGTACAGAACTAATGCGTAGTATTGGATGACCAACACCACTTTCATTTGGCTTTGATGAAAGACCATTACGTAATTCAGTTAAGATAGATTCAAA  >  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: