Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,504,880 A→C 8.4% intergenic (+97/+55) frlR → / ← yhfS putative DNA‑binding transcriptional regulator/FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,504,8800AC 8.4% 672.4 / 14.0 248intergenic (+97/+55)frlR/yhfSputative DNA‑binding transcriptional regulator/FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein
Reads supporting (aligned to +/- strand):  ref base A (86/129);  new base C (19/6);  total (107/141)
Fisher's exact test for biased strand distribution p-value = 1.01e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 257 total aligned reads displayed.

CTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTG  >  NC_000913/3504732‑3505029
                                                                                                                                                     |                                                                                                                                                     
 ttctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCaaa                                                                                                                                                     <  1:2252210/150‑1 (MQ=255)
 ttctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAAGTTTACATCCCGTACGTTCCCCGCACCCTAA‑CCCTCTCCCCcaa                                                                                                                                                     >  2:352745/1‑150 (MQ=255)
  tctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAg                                                                                                                                                    <  1:2779306/150‑1 (MQ=255)
   ctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGATATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCAAA‑CCCTCTCCCCAAAgg                                                                                                                                                   >  2:1699164/1‑150 (MQ=255)
    tGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAggg                                                                                                                                                  <  2:1814374/150‑1 (MQ=255)
        gcgcCAATCGCATCACCGGCCCTAGTGAGAATAAAAGACCGGCACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTCA‑CCCTCTCCCCCCCGGGGCGa                                                                                                                                              <  2:1252847/150‑1 (MQ=25)
           ccAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAggg                                                                                                                                           <  2:59571/150‑1 (MQ=255)
           ccAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAggg                                                                                                                                           >  2:3873971/1‑150 (MQ=255)
           ccAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCGCCCCTTTGGGGAGAggg                                                                                                                                           <  1:3701966/150‑1 (MQ=39)
           caaaTCGCATCACCTTAACTATTGATAATAAAAGACATTAACGGTAGTCCGGATATGTTGCGCCAGCATCGAATCGGGCGTTGCTGGCGGGTTGATTTTTACATCCCGTTCGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAggg                                                                                                                                           <  1:2648796/148‑1 (MQ=37)
            cAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCAAA‑CCCACACCCCAAAGGGGCGAgggg                                                                                                                                          >  1:1657554/1‑150 (MQ=255)
            cAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCTCACCCTAA‑CCCGCTACCCCAAGGGGCGAgggg                                                                                                                                          >  1:371974/1‑150 (MQ=255)
            cAATCGCATCACCTTAACAATTGATAATAAAAGACATTAACAGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCTAA‑CCCACACCCCAAAGGGGCGAgggg                                                                                                                                          >  2:2864007/1‑150 (MQ=39)
             aaTCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCGACCCCCAA‑CCCACTCCCCCCAGGGGCggggggg                                                                                                                                         >  1:1877344/1‑149 (MQ=38)
             aaTCGCATCACCTGAACTATTGATAATAAAAGACATTAACCGTAGGCCGGCTAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTGTACATCCCGTCCGGTCCCCGCCCCCAAC‑CCCGCACCCCCCACGGGCGCggggc                                                                                                                                         >  1:2322607/1‑149 (MQ=25)
              aTCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAgggtta                                                                                                                                        <  2:3154067/150‑4 (MQ=255)
              aTCGCATCACCTTAACGATTGATAATACAAGACATTCACCGGAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTGTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCGCCCCTTTGGGGAGAgggtta                                                                                                                                        <  2:2662309/150‑4 (MQ=38)
                  cATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCCCCCTAA‑ACCACTCCCCCAAGGGGCGAGGGGACCGtt                                                                                                                                    >  1:4983848/1‑150 (MQ=255)
                    tCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCa                                                                                                                                  <  2:3481659/150‑1 (MQ=255)
                         ttcaCTCTTGATAATAAAAGACATGCACCGGAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACg                                                                                                                             <  2:4504823/147‑1 (MQ=38)
                          tAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCAAA‑CCCACACCCCCAAAGGGCCAGGGGACCGTGCAATACgg                                                                                                                            >  1:1112763/1‑150 (MQ=38)
                            aCTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGaa                                                                                                                          <  1:4587470/150‑1 (MQ=255)
                              tATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAAtt                                                                                                                        <  2:3695332/150‑1 (MQ=255)
                                ttGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCAAA‑ACCACAACCCCAAGGGGCGAGGGGCCCGGTCAGTACGGAAATTa                                                                                                                      >  1:1945990/1‑150 (MQ=38)
                                ttGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCACAACCAAA‑CCCACTACCCAAATGGGCGAGGGGACCCTTCAGTACGGAATTTa                                                                                                                      >  1:3864422/1‑150 (MQ=38)
                                        aaaaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGGTGAATTTTACATCCCGTACGTTGCCCGCACCCAAA‑ACCACACCCCCAAGGGGCGGGGGGGCCCGTGAGTGAGGAAGTTAGCGTGTCg                                                                                                              >  1:973418/1‑150 (MQ=21)
                                         aaaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCAAA‑CCCTCTCCCCAAACGGGCGAGGGGACCGGTCAGTACGGAATTTAGCTTGACgg                                                                                                             >  2:4141818/1‑150 (MQ=255)
                                         aaaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCTAA‑CCCTCTACCCAAATGGGCGAGGGGACCGGTCAGTAACGAATTTAGCTTGACgg                                                                                                             >  2:673406/1‑150 (MQ=39)
                                           aGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTg                                                                                                           <  2:2363891/150‑1 (MQ=255)
                                            gACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTg                                                                                                           <  2:299478/150‑1 (MQ=255)
                                            gACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGc                                                                                                          <  2:2057266/150‑1 (MQ=255)
                                            gACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGGTGAATTTTACATCCCGTACGGTCCCCACACCCAAA‑CCCACACCCCCAACGGGCCCGGGGCCCGGTCAGTAAGGCAAGTAGCTTGACGGTGc                                                                                                          >  1:3392928/1‑150 (MQ=25)
                                                 tAACCGTCGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACCTCCCGTTCGTTCCCCTCACCCTAA‑CCCTCTCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATcc                                                                                                     <  1:716191/150‑1 (MQ=255)
                                                    ccGTTGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGc                                                                                                  <  1:1308734/150‑1 (MQ=255)
                                                       tAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGTATTTTAAATCCCGTACGGTACCCCCACCCGAA‑ACCAAAACCCAAAGGGGGGAGGGGGCCGGTGCGGACGGAAATTTTCTTGTAGGGGCTATCCAGccat                                                                                               >  1:4714468/1‑147 (MQ=16)
                                                          gCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTaaa                                                                                            <  2:2305820/150‑1 (MQ=255)
                                                             ggATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCTAA‑ACCAATCCCCCAAGGGGCGCGGGGACCGGTCAGTACCGAATTTAGCTTGACGGTGCTATCCAGCATTaaaaaa                                                                                         >  2:742891/1‑147 (MQ=255)
                                                              gATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACt                                                                                        <  2:2680708/150‑1 (MQ=255)
                                                               attcaacgcgcCAGCATCGCCGGATGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCACCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCATGACGGTGCTATCCAGCCGTAAAAACTg                                                                                       >  2:4702036/8‑150 (MQ=25)
                                                                tAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTgg                                                                                      <  2:899307/150‑1 (MQ=255)
                                                                tAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTgg                                                                                      <  2:810866/150‑1 (MQ=255)
                                                                            gcatcgcctcCGGCGATGCTGGCGCGTTGAATTGTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCAc                                                                          <  1:4951324/150‑1 (MQ=255)
                                                                            gcatcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCAc                                                                          <  1:4438183/150‑1 (MQ=255)
                                                                            gcatcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCAc                                                                          <  1:2502572/150‑1 (MQ=255)
                                                                            gcatcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCAAA‑ACCACTCCCCAAACGGGCGAGGGGACCGGTCAGTACGGAAATTATCTTGACGGTGCGAACCAGCAATAAATACAGGCAATaaacgcac                                                                          >  2:1441137/1‑143 (MQ=255)
                                                                            gcatcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCAc                                                                          <  1:1468397/150‑1 (MQ=255)
                                                                              atcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGc                                                                        <  2:2530335/150‑1 (MQ=255)
                                                                              atcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGc                                                                        <  2:2709894/150‑1 (MQ=255)
                                                                              atcgcatcCGGCGATGCTGGCGCGTTGAATTGTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGc                                                                        <  2:3344356/150‑1 (MQ=255)
                                                                                  catcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCACACCCAAA‑CCCACTCCCCAAAGGGGCGCGGGGCCCGGTCAGTACGGAATTTAGCGTGTCGGTGCTATCCAGCATTAAATACTGGCAATAATCTAACGCaaaa                                                                    >  2:4464342/1‑150 (MQ=255)
                                                                                       ggCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCgc                                                               <  1:3063253/150‑1 (MQ=255)
                                                                                       ggCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCAATAAATACTGGCAATACTCTCAAGCAAAATCCgc                                                               >  2:1182410/1‑150 (MQ=255)
                                                                                           aTGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCgcagca                                                           <  1:273938/150‑1 (MQ=255)
                                                                                           aTGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCgcagca                                                           <  2:5014854/150‑1 (MQ=255)
                                                                                            tGCTGGCGCGTTGAAATTTACATCCCGAACGTTCCCCTCACCCAAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACCCAAAATCCGCAGCAc                                                          <  1:1769217/150‑1 (MQ=255)
                                                                                               tGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGt                                                       <  2:4598745/150‑1 (MQ=255)
                                                                                                    cgTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCtt                                                  <  2:5004548/150‑1 (MQ=255)
                                                                                                         aaTTTTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCg                                             <  2:4417009/150‑1 (MQ=255)
                                                                                                           ttttGCATCCCGTACGTTCCCCTCACCCTAA‑CTCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTgg                                                                                      >  1:404618/1‑107 (MQ=255)
                                                                                                           ttttGCATCCCGTACGTTCCCCTCACCCTAA‑CTCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTgg                                                                                      <  2:404618/107‑1 (MQ=255)
                                                                                                            tttACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTg                                          <  1:4179052/150‑1 (MQ=255)
                                                                                                            ggtACATCCGGTACGTTCGGCTCCCGCTCA‑CGCTCGCCCCTCAGGGGCGAGGGGACCGTTCAGTCCGGTATTTCGCTTGACGGTGCTATCGAGCAGTAAATACTGGCTATTCTCTCACGCATAATCCGGTGCTCCGTGTCTTCACAGCtt                                          <  2:3645289/148‑2 (MQ=255)
                                                                                                             ttACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGc                                         <  1:3367103/150‑1 (MQ=255)
                                                                                                                  tCCCGTACGTTCCCCTCACCCAAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTc                                    <  1:1365027/150‑1 (MQ=255)
                                                                                                                    ccGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCgg                                  <  1:8496/150‑1 (MQ=255)
                                                                                                                    ccGTACGTTCCCCACACCCTAA‑ACCACAACCCAAAGGGGCGAGGGGACCGGTCAGTACGGAATTTAGCTTGACGGTGCTAACCAGCATTAAATAAAGGCAATACACTCACCCAAAATACGCAGCACCGTCTCTTCACCGCTGCGATTAgg                                  >  2:3212267/1‑150 (MQ=255)
                                                                                                                     cGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCggg                                 >  2:3430764/1‑150 (MQ=255)
                                                                                                                      gTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGt                                <  2:2167444/150‑1 (MQ=255)
                                                                                                                      gTACGTTCCCCTCACCCTAA‑ACCTCTCCCCAAAGGGGCGAGGGGACCGTTAATTAAGGAATTTTTGTTTAAGGTTCTATAAATCATTAAATAATGTCAATAATATAACGCAAAATACGAAGCAACGTATATTAACCGCTGAGAttatggt                                >  2:1940161/1‑145 (MQ=12)
                                                                                                                         cGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTAGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTaa                             <  2:2804332/150‑1 (MQ=255)
                                                                                                                         cGTTCCCCTCACCCTAA‑CCCACTCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTaa                             <  2:328118/150‑1 (MQ=255)
                                                                                                                          ggtCCCCTCACCCTCA‑CCCTCTCCCCCAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAAt                            <  1:3186872/148‑1 (MQ=255)
                                                                                                                          gTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAAt                            <  1:3509466/150‑1 (MQ=255)
                                                                                                                            tCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGc                          >  2:4345216/1‑150 (MQ=255)
                                                                                                                               ccTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGaa                       >  2:2391040/1‑150 (MQ=255)
                                                                                                                                cTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAAt                      >  1:1253578/1‑150 (MQ=255)
                                                                                                                                 tCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCTCTGCGATTCGGGTTAATGCGAATc                     >  1:2385775/1‑150 (MQ=255)
                                                                                                                                   aCCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGc                   >  2:1134002/1‑150 (MQ=255)
                                                                                                                                     ccTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGcac                 >  2:1275838/1‑150 (MQ=255)
                                                                                                                                     ccTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGcac                 >  2:3377491/1‑150 (MQ=255)
                                                                                                                                     ccTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGcac                 >  2:2434015/1‑150 (MQ=255)
                                                                                                                                         a‑CCCTCTCCCCAACGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAACATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATg             <  2:708057/150‑1 (MQ=255)
                                                                                                                                            ccTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGtt           >  1:4630698/1‑150 (MQ=255)
                                                                                                                                             ctctCCCCATAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTc          <  1:2649268/150‑1 (MQ=255)
                                                                                                                                             ctctCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCTATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTc          <  1:1017135/150‑1 (MQ=255)
                                                                                                                                             ctctCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTc          <  1:2952037/150‑1 (MQ=255)
                                                                                                                                             ctctCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTc          >  1:2110924/1‑150 (MQ=255)
                                                                                                                                             ctctCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTGGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTc          <  2:4118022/150‑1 (MQ=255)
                                                                                                                                              tctcCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCt         >  2:1034794/1‑150 (MQ=255)
                                                                                                                                              tctcCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCt         >  1:3564301/1‑150 (MQ=255)
                                                                                                                                                  cccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGAtt     >  2:4661595/1‑150 (MQ=255)
                                                                                                                                                  cccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGAtt     >  2:3794173/1‑150 (MQ=255)
                                                                                                                                                  cccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGAtt     >  2:2782290/1‑150 (MQ=255)
                                                                                                                                                  cccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGAtt     >  1:4535515/1‑150 (MQ=255)
                                                                                                                                                  cccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGAtt     >  2:180096/1‑150 (MQ=255)
                                                                                                                                                   ccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGTTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtg    >  1:544053/1‑150 (MQ=255)
                                                                                                                                                   ccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtg    >  2:2211929/1‑150 (MQ=255)
                                                                                                                                                    cAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtgt   <  2:1657554/150‑1 (MQ=255)
                                                                                                                                                    cAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtgt   >  1:219311/1‑150 (MQ=255)
                                                                                                                                                     aaaGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtgtg  >  2:407510/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAA‑CCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTG  >  NC_000913/3504732‑3505029

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑