Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,646,935 A→G intergenic (+198/+239) flxA → / ← ydfW Qin prophage; protein FlxA/Qin prophage; uncharacterized protein YdfW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,646,9350AG98.2% 391.2 / ‑2.8 112intergenic (+198/+239)flxA/ydfWQin prophage; protein FlxA/Qin prophage; uncharacterized protein YdfW
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (58/52);  total (59/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01

Only 100 of 113 total aligned reads displayed.

ATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATC  >  NC_000913/1646802‑1647070
                                                                                                                                     |                                                                                                                                       
attatCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACt                                                                                                                                   <  2:901152/140‑1 (MQ=255)
attatCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACt                                                                                                                                   <  1:1517440/140‑1 (MQ=255)
attatCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACt                                                                                                                                   >  1:1150421/1‑140 (MQ=255)
       ttGTATCATTTCTGACAGTTAACATGACACGATATAAATAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTgg                                                                                                                            <  1:410577/140‑1 (MQ=255)
       ttGTATCATTTCTGAAAGTTAATATTACAAGATATAAATAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTgg                                                                                                                            <  1:541351/140‑1 (MQ=255)
        tGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTggg                                                                                                                           >  1:814429/1‑140 (MQ=255)
        tGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTggg                                                                                                                           >  1:1166202/1‑140 (MQ=255)
        tGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTggg                                                                                                                           >  1:1084803/1‑140 (MQ=255)
        tGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTGTgggg                                                                                                                          >  1:781570/1‑140 (MQ=255)
          tATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGa                                                                                                                         >  2:1147722/1‑140 (MQ=255)
          tATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGa                                                                                                                         >  2:152628/1‑140 (MQ=255)
           aTCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGa                                                                                                                         >  1:1029098/1‑139 (MQ=255)
           aTCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGa                                                                                                                         <  2:1029098/139‑1 (MQ=255)
           aTCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGAc                                                                                                                        >  2:869182/1‑140 (MQ=255)
           aTCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGAc                                                                                                                        >  2:1184821/1‑140 (MQ=255)
            tCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACa                                                                                                                       <  2:1150421/140‑1 (MQ=255)
              aTTTCTGATAGTTAACATTACAAGATATAAGTAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGt                                                                                                                     >  1:1350116/1‑140 (MQ=255)
              aTTTCTGATAGTTAACATTACAAGATATAAGTAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGt                                                                                                                     >  1:1403630/1‑140 (MQ=255)
                         ttAACATTACAAGATATAAGTAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTg                                                                                                          >  1:1264544/1‑140 (MQ=255)
                                  cAAGATATAAGTAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTcc                                                                                                 >  2:1095751/1‑140 (MQ=255)
                                  cAAGATATAAGTAATCGACGCACTCCCAACTAGTCTAATTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTcc                                                                                                 <  1:411574/140‑1 (MQ=255)
                                   aaGATATAAGTAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCg                                                                                                >  1:224794/1‑140 (MQ=255)
                                   aaGATATAAGTAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCg                                                                                                >  2:238138/1‑140 (MQ=255)
                                     gATATAAGTAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGct                                                                                              >  2:106217/1‑140 (MQ=255)
                                       tataAGTAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGcttc                                                                                            >  2:309007/1‑140 (MQ=255)
                                          aaGTAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGcc                                                                                                        <  1:974259/125‑1 (MQ=255)
                                            gTAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGct                                                                                       >  1:874171/1‑140 (MQ=255)
                                            gTAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGct                                                                                       <  1:1083540/140‑1 (MQ=255)
                                            agaaTGGACGCACTCCCAATTATTCTATTTAAATCTCCATGAGTTTAAATGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGct                                                                                       <  2:526117/138‑1 (MQ=255)
                                               aTGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCtt                                                                                    >  2:1092883/1‑140 (MQ=255)
                                               aTGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACATTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCtt                                                                                    >  2:1252278/1‑140 (MQ=255)
                                                tGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACTAAACTGTCAACTTCCGCTTCTTGCTCTTa                                                                                   >  1:572710/1‑140 (MQ=255)
                                                tGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTa                                                                                   >  1:1314730/1‑140 (MQ=255)
                                                tGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTa                                                                                   >  2:288013/1‑140 (MQ=255)
                                                tGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTa                                                                                   >  1:1181514/1‑140 (MQ=255)
                                                 ggACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAg                                                                                  >  1:286430/1‑140 (MQ=255)
                                                 ggACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAg                                                                                  >  1:279729/1‑140 (MQ=255)
                                                  gACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCACTTCTTGCTCTTAGc                                                                                 >  1:673177/1‑140 (MQ=255)
                                                   gacagcTCCCAATTAGTATATTCAAATCGCCTCGACTTTTACTGACAAACCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCg                                                                                <  2:477397/135‑1 (MQ=255)
                                                   aCGCACTCCCAATTAGTTTATTTAAATCGCCGCGACTTTCACTGACAAGCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCg                                                                                <  2:450883/140‑1 (MQ=255)
                                                   aCGCACTCCCAATTAGTGTATTTAAATCGCCTCGACTTTAACTGACAATCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCg                                                                                <  2:673177/140‑1 (MQ=255)
                                                   aCGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCg                                                                                <  2:1159592/140‑1 (MQ=255)
                                                   aCGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCg                                                                                <  2:1408510/140‑1 (MQ=255)
                                                   aCGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCg                                                                                <  2:279729/140‑1 (MQ=255)
                                                   aCGCACTCCCAATTAGCCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGAACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCg                                                                                <  2:931647/140‑1 (MQ=255)
                                                    cGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCgg                                                                               >  2:1362925/1‑140 (MQ=255)
                                                     gCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGa                                                                              <  2:808236/140‑1 (MQ=255)
                                                       aCTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGct                                                                                       >  1:262672/1‑129 (MQ=255)
                                                       aCTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGct                                                                                       <  2:262672/129‑1 (MQ=255)
                                                       aCTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACa                                                                            <  2:1283530/140‑1 (MQ=255)
                                                         tCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCGACTACCTCTTCTTCCTCTTAGCGGACact                                                                          >  1:526117/1‑138 (MQ=255)
                                                            cAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGc                                                                       <  2:1314730/140‑1 (MQ=255)
                                                               ttAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATa                                                                    <  1:198546/140‑1 (MQ=255)
                                                               ttAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATa                                                                    <  2:1403630/140‑1 (MQ=255)
                                                                 aGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAgg                                                                  <  1:1147722/140‑1 (MQ=255)
                                                                 aGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAgg                                                                  <  2:1264544/140‑1 (MQ=255)
                                                                 aGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAgg                                                                  <  1:106217/140‑1 (MQ=255)
                                                                 aGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAgg                                                                  <  1:309007/140‑1 (MQ=255)
                                                                 aGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAgg                                                                  <  2:286430/140‑1 (MQ=255)
                                                                   tCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCg                                                                                >  2:311261/1‑124 (MQ=255)
                                                                   tCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCAt                                                                     >  2:1410918/1‑135 (MQ=255)
                                                                   tCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCAt                                                                     <  1:1410918/135‑1 (MQ=255)
                                                                   tCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATa                                                                    >  1:1003344/1‑136 (MQ=255)
                                                                   tCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATa                                                                    <  2:1003344/136‑1 (MQ=255)
                                                                   tCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCt                                                                >  2:1386535/1‑140 (MQ=255)
                                                                   tCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAAATTCCGCTTCTTGCTGTTAGCGTACATTAGCAtctgtc                                                                >  1:477397/1‑135 (MQ=255)
                                                                     tATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACa                                                                            <  2:854917/126‑1 (MQ=255)
                                                                     tATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACa                                                                            >  1:854917/1‑126 (MQ=255)
                                                                     tATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTa                                                               <  2:862065/139‑1 (MQ=255)
                                                                      aTTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTAtt                                                             <  1:152628/140‑1 (MQ=255)
                                                                      aTTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTAtt                                                             <  1:1095751/140‑1 (MQ=255)
                                                                       tttAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAATTGCCAACTTCCGCTACTTGCTCTTAGCGGATATTAGCATAGGCTAttt                                                            >  2:667312/1‑140 (MQ=255)
                                                                       tttAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTAttt                                                            <  2:1166202/140‑1 (MQ=255)
                                                                       tttAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTAttt                                                            <  2:112081/140‑1 (MQ=255)
                                                                               gCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAAc                                                    >  1:804510/1‑140 (MQ=255)
                                                                                ccACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACg                                                   <  1:1386535/140‑1 (MQ=255)
                                                                                 cACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGc                                                  >  1:312177/1‑140 (MQ=255)
                                                                                    gAGTTTAACTGACAACCCATGAGCAATTATGAATTGCAATTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTc                                               <  2:637033/140‑1 (MQ=255)
                                                                                     aGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCa                                              >  2:1256188/1‑140 (MQ=255)
                                                                                             tGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACgcgc                                      <  2:312177/140‑1 (MQ=255)
                                                                                              gACAACCCATGATCACTTATGACTTACAAGTATTTCTCTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCa                                     <  1:746135/140‑1 (MQ=255)
                                                                                              gACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCa                                     >  2:158557/1‑140 (MQ=255)
                                                                                              gACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCa                                     >  1:715613/1‑140 (MQ=255)
                                                                                              gACAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCAAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCa                                     <  1:820430/140‑1 (MQ=255)
                                                                                               aCAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCAc                                    >  1:1215392/1‑140 (MQ=255)
                                                                                                cAACCCATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATCAGCCTAGACTATTTACGATAAGGCCTCATTACGCGCAcc                                   >  2:433349/1‑140 (MQ=255)
                                                                                                   cccATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGcc                                <  1:1256188/140‑1 (MQ=255)
                                                                                                   cccATGATCAATTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGcc                                <  2:874171/140‑1 (MQ=255)
                                                                                                             aTTATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAgactgact                      <  2:1124525/140‑1 (MQ=255)
                                                                                                              ttATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCg                                  >  1:864763/1‑127 (MQ=255)
                                                                                                              ttATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTc                     >  1:1173830/1‑140 (MQ=255)
                                                                                                              ttATGAATTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCGCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCg                                  <  2:864763/127‑1 (MQ=255)
                                                                                                                  gAATTGCAACTATTTCTGTAGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAgcg                 <  1:238138/140‑1 (MQ=255)
                                                                                                                   aaTTGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAgcgc                >  2:1075959/1‑140 (MQ=255)
                                                                                                                      tGCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGtt             >  2:160881/1‑140 (MQ=255)
                                                                                                                       gCAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGttt            <  2:1215392/140‑1 (MQ=255)
                                                                                                                        cAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTc           >  1:1369667/1‑140 (MQ=255)
                                                                                                                        cAACTATTTCTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTc           >  1:1075280/1‑140 (MQ=255)
                                                                                                                                 cTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATc  <  1:324525/140‑1 (MQ=255)
                                                                                                                                 cTGTGGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATc  <  1:1075959/140‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                       
ATTATCTTTGTATCATTTCTGATAGTTAACATTACAAGATATAAGTAATGGACGCACTCCCAATTAGTCTATTTAAATCGCCACGAGTTTAACTGACAACCCATGATCAATTATGAATTGCAACTATTTCTGTAGTCACTTTTGTGGGGACAGTCCACAAAACTGCCAACTTCCGCTTCTTGCTCTTAGCGGACATTAGCATAGGCTATTTACCATAACGCCTCATTACGCGCACCGCCCAGACTGACTCAGCGCGTTTCTGGCATATC  >  NC_000913/1646802‑1647070

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