Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,599,756 A→G 6.5% intergenic (‑98/+173) livJ ← / ← rpoH branched chain amino acid/phenylalanine ABC transporter periplasmic binding protein/RNA polymerase, sigma 32 (sigma H) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,599,7560AG 6.5% 370.6 / 11.8 122intergenic (‑98/+173)livJ/rpoHbranched chain amino acid/phenylalanine ABC transporter periplasmic binding protein/RNA polymerase, sigma 32 (sigma H) factor
Reads supporting (aligned to +/- strand):  ref base A (49/65);  new base G (4/4);  total (53/69)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

Only 100 of 125 total aligned reads displayed.

TGAATGCCAGCGCGATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAG  >  NC_000913/3599607‑3599905
                                                                                                                                                     |                                                                                                                                                     
 gAATGCCAGCGCGATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCAt                                                                                                                                                      >  1:1454218/1‑150 (MQ=255)
     gCCAGCGCGATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATActc                                                                                                                                                  >  1:935425/1‑150 (MQ=255)
      ccAGCGCGATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATActct                                                                                                                                                 >  1:200710/1‑150 (MQ=255)
        aGCGCGATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGc                                                                                                                                               <  1:1455982/150‑1 (MQ=255)
         gcgcGATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCt                                                                                                                                              >  2:445157/1‑150 (MQ=255)
             gATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACACGTGCTTTAGCAGCATACTCTGCTaaaa                                                                                                                                          >  1:1405929/1‑150 (MQ=255)
              aTACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAAc                                                                                                                                         <  2:1125186/150‑1 (MQ=255)
              aTACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAAc                                                                                                                                         <  2:1134079/150‑1 (MQ=255)
               tACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACa                                                                                                                                        <  1:216438/150‑1 (MQ=255)
               tACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACa                                                                                                                                        <  1:1295763/150‑1 (MQ=255)
                 cATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATa                                                                                                                                      >  1:906100/1‑150 (MQ=255)
                  aTCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATAc                                                                                                                                     <  2:788541/150‑1 (MQ=255)
                   cccTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGGGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATAcc                                                                                                                                    <  2:181936/149‑1 (MQ=255)
                      tGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGa                                                                                                                                 >  2:1022494/1‑150 (MQ=255)
                         cAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGAttt                                                                                                                              <  2:821021/150‑1 (MQ=255)
                         cAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTAAGCACTTATTTCGTTTTATGCATGACTACCCGTGCTTTACCAGCATACTCTGCTAAAACATACCCGAtttt                                                                                                                             >  2:206823/1‑150 (MQ=255)
                            tAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTa                                                                                                                           <  2:1306951/150‑1 (MQ=255)
                            tAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTa                                                                                                                           <  1:1283243/150‑1 (MQ=255)
                            tAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTa                                                                                                                           <  1:1303417/150‑1 (MQ=255)
                              aCGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATg                                                                                                                         >  1:730340/1‑150 (MQ=255)
                                gCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGat                                                                                                                       >  2:292725/1‑150 (MQ=255)
                                    tACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTg                                                                                                                   >  2:1503755/1‑150 (MQ=255)
                                     aCCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTgg                                                                                                                  <  1:201359/150‑1 (MQ=255)
                                     aCCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTgg                                                                                                                  >  1:737960/1‑150 (MQ=255)
                                     aCCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACATGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTgg                                                                                                                  <  1:1436342/150‑1 (MQ=255)
                                      cccTTTATGTTCATCCTGATAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGa                                                                                                                 <  1:473857/150‑1 (MQ=255)
                                      cccTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGa                                                                                                                 <  1:1203401/150‑1 (MQ=255)
                                      cccTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGa                                                                                                                 >  2:782534/1‑150 (MQ=255)
                                      cccTTTATGTTCATCCTGAGAATCACCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGa                                                                                                                 <  1:1008117/150‑1 (MQ=255)
                                       ccTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGaa                                                                                                                >  2:721669/1‑150 (MQ=255)
                                         tttATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATa                                                                                                              <  1:891124/150‑1 (MQ=255)
                                          ttATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAg                                                                                                             <  1:849311/150‑1 (MQ=255)
                                           caTGGTCATCCTCAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTGTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGc                                                                                                            <  2:31805/149‑1 (MQ=255)
                                            aTGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGTTTTAGCAGAGTATGCTGCTAAAGCacg                                                                                                                                      <  1:1146087/123‑3 (MQ=38)
                                            aTGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCt                                                                                                           <  2:737960/150‑1 (MQ=255)
                                               ttCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTAtt                                                                                                        <  1:1349032/150‑1 (MQ=255)
                                                tCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTAttt                                                                                                       <  2:935425/150‑1 (MQ=255)
                                                 cATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTAtttt                                                                                                      <  1:1045821/150‑1 (MQ=255)
                                                  atgcTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTg                                                                                                     <  2:321775/147‑1 (MQ=255)
                                                  aTCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTg                                                                                                     <  1:764818/150‑1 (MQ=255)
                                                    ccTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGAc                                                                                                   >  2:1061938/1‑150 (MQ=255)
                                                     cttagtaTCCCCATTTTTCTGGTTATTACGTTTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTTTGGATGTCTACCCGTGGTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACa                                                                                                  <  2:1405929/144‑1 (MQ=255)
                                                     cTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACa                                                                                                  >  2:1439650/1‑150 (MQ=255)
                                                     cTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACa                                                                                                  >  1:1040542/1‑150 (MQ=255)
                                                      tGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGAttttt                                                                                                                            >  1:947649/1‑123 (MQ=255)
                                                      tGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGAttttt                                                                                                                            <  2:947649/123‑1 (MQ=255)
                                                      tGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAg                                                                                                 <  1:1200568/150‑1 (MQ=255)
                                                      tGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAg                                                                                                 <  2:317521/150‑1 (MQ=255)
                                                      tGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAg                                                                                                 >  1:1509175/1‑150 (MQ=255)
                                                         gaATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGttt                                                                                              <  2:397926/150‑1 (MQ=255)
                                                          aaTCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTtta                                                                                             <  1:131426/150‑1 (MQ=255)
                                                           aTCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTttat                                                                                            <  1:375522/150‑1 (MQ=255)
                                                           aTCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTttat                                                                                            <  1:1022494/150‑1 (MQ=255)
                                                             ccccATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTttatta                                                                                          <  2:893531/150‑1 (MQ=255)
                                                              cccATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTaa                                                                                         >  1:1269041/1‑150 (MQ=255)
                                                              cccATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTaa                                                                                         >  1:1271014/1‑150 (MQ=255)
                                                                 aTTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACaa                                                                                      <  1:1263331/150‑1 (MQ=255)
                                                                  ttcttcTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAAt                                                                                     <  2:1454218/150‑1 (MQ=255)
                                                                   tcttctGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATc                                                                                    <  1:449201/150‑1 (MQ=255)
                                                                           ttattaCGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTggg                                                                            <  2:730340/150‑1 (MQ=255)
                                                                             attaCGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGa                                                                          >  2:1286139/1‑150 (MQ=255)
                                                                             attaCGTGTGTTGTGATGTTGTTTTTAAGCACTTTATTTCGTTTTATGCATGACTACCCGTTCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGa                                                                          <  1:1485699/150‑1 (MQ=255)
                                                                                 cGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTgg                                                                      >  1:108431/1‑150 (MQ=255)
                                                                                  gtgtgtTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATc                                                                                    >  2:1338483/1‑135 (MQ=255)
                                                                                     tgttgtGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGtt                                                                  <  2:234421/150‑1 (MQ=255)
                                                                                      gttgtGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGttt                                                                 <  2:706535/150‑1 (MQ=255)
                                                                                        tgtgATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGGTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTg                                                               >  2:522298/1‑150 (MQ=255)
                                                                                          tgATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGttt                                                                 <  2:107790/145‑1 (MQ=255)
                                                                                          tgATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGttt                                                                 >  1:107790/1‑145 (MQ=255)
                                                                                              gttgttTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAg                                                         >  1:1116728/1‑150 (MQ=255)
                                                                                               ttgttTTTCAGCACTTTATTTCGTTTTATGCATTACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAgg                                                        <  1:275463/150‑1 (MQ=255)
                                                                                               ttgttTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAgg                                                        <  1:1503755/150‑1 (MQ=255)
                                                                                                tgttTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAggg                                                       <  1:710023/150‑1 (MQ=255)
                                                                                                tgttTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAggg                                                       >  2:691578/1‑150 (MQ=255)
                                                                                                 gttTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAgggg                                                      >  1:931344/1‑150 (MQ=255)
                                                                                                  tttttCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGa                                                     <  2:1509175/150‑1 (MQ=255)
                                                                                                    tttCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGaaa                                                   >  2:1156069/1‑150 (MQ=255)
                                                                                                      tCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAAtt                                                 <  1:853657/150‑1 (MQ=255)
                                                                                                        aGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGa                                               >  1:1340933/1‑150 (MQ=255)
                                                                                                          cACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGatt                                             <  2:1116728/150‑1 (MQ=255)
                                                                                                            cTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTGTTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGattat                                           <  1:974396/150‑1 (MQ=255)
                                                                                                                aTTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACa                                       >  2:275013/1‑150 (MQ=255)
                                                                                                                aTTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACa                                       >  1:1006311/1‑150 (MQ=255)
                                                                                                                 tttCTTTTTATGGTTGTCTACCGGTGCTTTAGCAGAATACTATGATAAAAGATACAAGATTTTTATGATATTGGATTATCTATTTTGACAGTTTATTAACAATCTGGGTGGGGATTGGAGTTTTGCCGGATGGGAAATTGATTATTACag                                      <  2:228243/150‑1 (MQ=255)
                                                                                                                 tttCGTTTTATGCATGACTACCCGAGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACag                                      >  2:1518843/1‑150 (MQ=255)
                                                                                                                   tCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGatta                                            <  1:413318/142‑1 (MQ=255)
                                                                                                                   tCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGatta                                            >  2:413318/1‑142 (MQ=255)
                                                                                                                    cGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAgg                                   >  1:867220/1‑150 (MQ=255)
                                                                                                                     gTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGc                                  >  1:188939/1‑150 (MQ=255)
                                                                                                                       tttATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGccc                                <  1:340689/150‑1 (MQ=255)
                                                                                                                            gCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATAGTGGAATAGCTATTTTGACAGTTTAGTAACAATCGGCGTGGGGATTGGCGTGTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCaaaaa                           >  1:1502789/1‑150 (MQ=255)
                                                                                                                                     ggtatgttttAGCAGAGTATGCTGCTAAAGCATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAAcc                  >  1:1004764/8‑150 (MQ=255)
                                                                                                                                        gTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCg               <  1:1369502/150‑1 (MQ=255)
                                                                                                                                         tGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAaaaaacaaa                      >  1:1364454/1‑142 (MQ=255)
                                                                                                                                         tGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAaaaaacaaa                      <  2:1364454/142‑1 (MQ=255)
                                                                                                                                          gttttAGCAGAGTATGCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGa             >  2:676266/3‑150 (MQ=255)
                                                                                                                                           cTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCTGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGAc            <  2:1040542/150‑1 (MQ=255)
                                                                                                                                               agcagcATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACtctc        >  1:1316652/1‑150 (MQ=255)
                                                                                                                                                  agcaTACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATc     >  1:487631/1‑150 (MQ=255)
                                                                                                                                                     aTACTCTGCTAAAACATACCCGATTTTTATGATATTGGAACAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAg  <  2:1056637/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGAATGCCAGCGCGATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAG  >  NC_000913/3599607‑3599905

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