Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,599,760 T→G 6.5% intergenic (‑102/+169) 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,7600TG 6.5% 392.9 / 16.0 124intergenic (‑102/+169)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 T (49/67);  new base G (4/4);  total (53/71)
Fisher's exact test for biased strand distribution p-value = 7.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-01

Only 100 of 127 total aligned reads displayed.

GCCAGCGCGATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGT  >  NC_000913/3599612‑3599908
                                                                                                                                                    |                                                                                                                                                    
gCCAGCGCGATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATActc                                                                                                                                                     >  1:935425/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:1134079/150‑1 (MQ=255)
          tACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACa                                                                                                                                           <  1:216438/150‑1 (MQ=255)
          tACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACa                                                                                                                                           <  1:1295763/150‑1 (MQ=255)
           aCATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGAGTATGCTGCTAAAGCac                                                                                                                                          <  2:28342/150‑2 (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                                                                                                                              <  1:1303417/150‑1 (MQ=255)
                       tAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTa                                                                                                                              <  1:1283243/150‑1 (MQ=255)
                       tAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTa                                                                                                                              <  1:1288558/150‑1 (MQ=255)
                        aaCGCTTTACCCTTTACGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTGGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTAt                                                                                                                             <  1:1208385/150‑1 (MQ=255)
                         aCGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATg                                                                                                                            >  1:730340/1‑150 (MQ=255)
                           gCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGat                                                                                                                          >  2:292725/1‑150 (MQ=255)
                           gCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGat                                                                                                                          >  1:234421/1‑150 (MQ=255)
                            cTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGata                                                                                                                         <  1:1405553/150‑1 (MQ=255)
                               tACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTg                                                                                                                      >  2:1503755/1‑150 (MQ=255)
                                aCCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTgg                                                                                                                     >  1:737960/1‑150 (MQ=255)
                                aCCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTgg                                                                                                                     <  1:201359/150‑1 (MQ=255)
                                 cccTTTATGTTCATCCTGATAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGa                                                                                                                    <  1:473857/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)
                                  ccTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGaa                                                                                                                   >  2:674713/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                                                                                                                                         >  2:1146087/1‑121 (MQ=38)
                                       aTGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGTTTTAGCAGAGTATGCTGCTAAAGCacg                                                                                                                                         <  1:1146087/123‑3 (MQ=38)
                                          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)
                                              tCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGTAATAGCTATTTTGa                                                                                                       <  2:200710/150‑1 (MQ=255)
                                                cttagtaTCCCCATTTTTCTGGTTATTACGTTTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTTTGGATGTCTACCCGTGGTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACa                                                                                                     <  2:1405929/144‑1 (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:1509175/1‑150 (MQ=255)
                                                 tGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAg                                                                                                    <  1:1200568/150‑1 (MQ=255)
                                                    gaATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGttt                                                                                                 <  2:397926/150‑1 (MQ=255)
                                                     aaTCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTtta                                                                                                <  2:8127/150‑1 (MQ=255)
                                                      aTCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTttat                                                                                               <  1:1022494/150‑1 (MQ=255)
                                                      aTCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTttat                                                                                               <  1:375522/150‑1 (MQ=255)
                                                         cccATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTaa                                                                                            >  1:1269041/1‑150 (MQ=255)
                                                            aTTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACaa                                                                                         <  1:1263331/150‑1 (MQ=255)
                                                             ttcttcTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAAt                                                                                        <  2:1454218/150‑1 (MQ=255)
                                                              tcttctGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGATTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATc                                                                                       <  1:921533/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)
                                                                            cGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTgg                                                                         <  2:1345418/150‑1 (MQ=255)
                                                                             gtgtgtTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATc                                                                                       <  1:1338483/135‑1 (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)
                                                                                          ttgttTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAgg                                                           >  1:1474611/1‑150 (MQ=255)
                                                                                           tgttTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAggg                                                          <  1:710023/150‑1 (MQ=255)
                                                                                           tgttTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAggg                                                          >  2:1485949/1‑150 (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)
                                                                                                    gCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGat                                                 >  2:101669/1‑150 (MQ=255)
                                                                                                     cACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGatt                                                <  2:1116728/150‑1 (MQ=255)
                                                                                                           aTTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACa                                          >  2:275013/1‑150 (MQ=255)
                                                                                                           aTTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACa                                          >  1:1006311/1‑150 (MQ=255)
                                                                                                              tCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGatta                                               <  1:413318/142‑1 (MQ=255)
                                                                                                               cGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAgg                                      >  1:867220/1‑150 (MQ=255)
                                                                                                                gTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGc                                     >  1:188939/1‑150 (MQ=255)
                                                                                                                                ggtatgttttAGCAGAGTATGCTGCTAAAGCATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAAcc                     >  1:1004764/8‑150 (MQ=255)
                                                                                                                                   gTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCg                  >  2:328393/1‑150 (MQ=255)
                                                                                                                                   gTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCg                  <  1:1369502/150‑1 (MQ=255)
                                                                                                                                    tGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAaaaaacaaa                         <  2:1364454/142‑1 (MQ=255)
                                                                                                                                      ttttAGCAGAGTATGCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGAc               >  1:162020/2‑150 (MQ=255)
                                                                                                                                      cTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCTGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGAc               <  2:1040542/150‑1 (MQ=255)
                                                                                                                                      cTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGAc               >  2:1070693/1‑150 (MQ=255)
                                                                                                                                          agcagcATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACtctc           >  1:1316652/1‑150 (MQ=255)
                                                                                                                                             agcaTACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATc        >  1:487631/1‑150 (MQ=255)
                                                                                                                                                aTACTCTGCTAAAACATACCCGATTTTTATGATATTGGAACAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAg     <  2:1056637/150‑1 (MQ=255)
                                                                                                                                                 tACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAgg    >  1:1159456/1‑150 (MQ=255)
                                                                                                                                                 tACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAgg    >  2:953985/1‑150 (MQ=255)
                                                                                                                                                  aCTCTTCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAggg   <  1:347502/150‑1 (MQ=255)
                                                                                                                                                   ctctGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGt  <  1:1439650/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCCAGCGCGATACATCCTGCCAGTAACGCTTTACCCTTTATGTTCATCCTGAGAATCCCCATTCTTCTGGTTATTACGTGTGTTGTGATGTTGTTTTTCAGCACTTTATTTCGTTTTATGCATGACTACCCGTGCTTTAGCAGCATACTCTGCTAAAACATACCCGATTTTTATGATATTGGAATAGCTATTTTGACAGTTTATTAACAATCTGCGTGGGGATTGGCGTTTTGCCGGAGGGGAAATTGATTATTACAGAGGCCCAAAAAACAAAAACCCCGGACTCTCATCCAGGGT  >  NC_000913/3599612‑3599908

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