Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,464,227 G→A 7.8% intergenic (+157/+25) yfdF → / ← mlaA uncharacterized protein/ABC transporter maintaining OM lipid asymmetry, OM lipoprotein component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,464,2270GA 7.8% 305.3 / 9.1 102intergenic (+157/+25)yfdF/mlaAuncharacterized protein/ABC transporter maintaining OM lipid asymmetry, OM lipoprotein component
Reads supporting (aligned to +/- strand):  ref base G (51/43);  new base A (5/3);  total (56/46)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.96e-01

Only 100 of 103 total aligned reads displayed.

TTTCCAGTTTTTTGATTTTTCAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCTGGAAGTACGCTTCGCGCACCAT  >  NC_000913/2464078‑2464374
                                                                                                                                                     |                                                                                                                                                   
tttCCAGTTTTTTGATTTTTCAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTg                                                                                                                                                     <  2:584637/150‑1 (MQ=255)
  tCCAGTTTTTTGATTTTTCAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCg                                                                                                                                                   >  1:1145639/1‑150 (MQ=255)
       ttttttGATTTTTCAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCa                                                                                                                                              >  1:781982/1‑150 (MQ=255)
             taTTTTTCATCAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTTGGTCATTTTCGGCTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCtttt                                                                                                                                        <  2:1494263/149‑1 (MQ=255)
             taTTTTTCAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACAAAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCtttt                                                                                                                                        <  2:806401/149‑1 (MQ=255)
             gATTTTTCAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCtttt                                                                                                                                        >  2:1374413/1‑150 (MQ=255)
              aTTTTTCAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTCGCAAAACTCACCttttt                                                                                                                                       <  1:1117178/150‑1 (MQ=255)
               tttttCAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTCGCAAAACTCACCtttttt                                                                                                                                      <  2:1063415/150‑1 (MQ=255)
                 tttCAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTCGCAAAACTCACCttttttgt                                                                                                                                    <  2:371887/150‑3 (MQ=255)
                  ttCAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTAtt                                                                                                                                   <  1:1304968/150‑1 (MQ=255)
                   tCAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTCAGTTTAGATACTCACCTTTTTTAttt                                                                                                                                  >  2:110726/1‑150 (MQ=255)
                    cAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTg                                                                                                                                 <  1:1249445/150‑1 (MQ=255)
                        aGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTc                                                                                                                             <  2:372266/150‑1 (MQ=255)
                         gCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCg                                                                                                                                                   >  1:288643/1‑127 (MQ=255)
                         gCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCg                                                                                                                                                   <  2:288643/127‑1 (MQ=255)
                           tGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGTGATAACAACAAAAAAGTTGAGTTTTGCGACTCACATTTTTTATTTGTTTCTTa                                                                                                                          >  1:242265/1‑150 (MQ=255)
                                atCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAg                                                                                                                     >  2:477471/1‑150 (MQ=255)
                                 tCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGa                                                                                                                    >  2:1018357/1‑150 (MQ=255)
                                              gTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCtttt                                                                                                       >  2:692345/1‑150 (MQ=255)
                                              gTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCtttt                                                                                                       >  2:1130563/1‑150 (MQ=255)
                                                gACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCTACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTaa                                                                                                     >  1:1587173/1‑150 (MQ=255)
                                                gACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTaa                                                                                                     <  1:801060/150‑1 (MQ=255)
                                                gACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTaa                                                                                                     >  1:1303513/1‑150 (MQ=255)
                                                gACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTaa                                                                                                     <  1:741343/150‑1 (MQ=255)
                                                gACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTAAGAATCAATATCTTTTaa                                                                                                     >  2:769096/1‑150 (MQ=255)
                                                 aCGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTaaa                                                                                                    >  2:1284401/1‑150 (MQ=255)
                                                  cGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAat                                                                                                   >  2:388401/1‑150 (MQ=255)
                                                       aTCATTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGTCTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATcc                                                                                              <  1:1173678/150‑1 (MQ=255)
                                                       aTCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATcc                                                                                              <  2:886803/150‑1 (MQ=255)
                                                          aGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGa                                                                                           >  2:372928/1‑150 (MQ=255)
                                                          aGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGa                                                                                           >  1:318116/1‑150 (MQ=255)
                                                          aGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGa                                                                                           <  2:424087/150‑1 (MQ=255)
                                                           gTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGaa                                                                                          >  2:1572687/1‑150 (MQ=255)
                                                                tttACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCt                                                                                     <  1:1215328/150‑1 (MQ=255)
                                                                tttACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCt                                                                                     <  1:692345/150‑1 (MQ=255)
                                                                 ttACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAAAAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCtt                                                                                    <  1:845027/150‑1 (MQ=255)
                                                                  tACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTtg                                                                                   <  1:1304467/150‑1 (MQ=255)
                                                                  tACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTtg                                                                                   <  2:781982/150‑1 (MQ=255)
                                                                   aCATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTtgt                                                                                  >  2:611106/1‑150 (MQ=255)
                                                                     aTAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGc                                                                                >  2:582034/1‑150 (MQ=255)
                                                                      tAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCg                                                                               <  1:480765/150‑1 (MQ=255)
                                                                        aTTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGtt                                                                             >  1:52022/1‑150 (MQ=255)
                                                                          tACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCg                                                                           <  1:880569/150‑1 (MQ=255)
                                                                           aCGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGAGTTCgt                                                                          >  1:235207/1‑149 (MQ=255)
                                                                               cTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGAttt                                                                      >  2:37754/1‑150 (MQ=255)
                                                                                 aGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTTTTTCTTATTAAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCTTTCGGATTTTc                                                                    >  1:1293106/1‑150 (MQ=255)
                                                                                 aGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTc                                                                    >  1:1119110/1‑150 (MQ=255)
                                                                                 aGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGATTGTGCGTTCGGATTTTc                                                                    >  2:1267441/1‑150 (MQ=255)
                                                                                 aGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATACTGAATCGCTTGTGCGTTCGGATTTTc                                                                    >  1:534606/1‑150 (MQ=255)
                                                                                        tGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGt                                                             <  2:151526/150‑1 (MQ=255)
                                                                                          tCGGTTTATCTGTTTTTAATTGTGGTTATTGAGTGATATCAATAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGttt                                                           <  1:110726/150‑1 (MQ=255)
                                                                                           cGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACGACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTg                                                          >  1:569263/1‑150 (MQ=255)
                                                                                            ggTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGa                                                         >  1:1660721/1‑150 (MQ=255)
                                                                                            ggTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGa                                                         >  1:100519/1‑150 (MQ=255)
                                                                                            ggTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGa                                                         <  2:36284/150‑1 (MQ=255)
                                                                                                tATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGtta                                                     >  1:1303769/1‑149 (MQ=255)
                                                                                                  tctcCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATAAATATCTTTTAAATCATCATGAATCGCGTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTcgc                                                   >  2:256771/1‑150 (MQ=255)
                                                                                                    tcCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTcgccg                                                 <  2:1660721/150‑1 (MQ=255)
                                                                                                    tcCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTcgccg                                                 >  2:1547512/1‑150 (MQ=255)
                                                                                                     ccTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTcgccgc                                                <  1:391018/150‑1 (MQ=255)
                                                                                                       ttttACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCa                                              <  2:1660337/150‑1 (MQ=255)
                                                                                                         ttACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCAtt                                            >  1:600780/1‑150 (MQ=255)
                                                                                                          tACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTa                                           >  1:393504/1‑150 (MQ=255)
                                                                                                             ttGTGGTAACTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCGGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCg                                        <  2:1163257/150‑1 (MQ=255)
                                                                                                              tgtgGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGa                                       >  1:1695103/1‑150 (MQ=255)
                                                                                                              tgtgGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCATGCGGTTTGAGTTCGCCGCCATTAGCGa                                       >  2:644529/1‑150 (MQ=255)
                                                                                                               gtgGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGAt                                      >  1:1409334/1‑150 (MQ=255)
                                                                                                                 ggTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGa                                    <  1:37754/150‑1 (MQ=255)
                                                                                                                  gTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGAAGCCATTAGCGATGaa                                   >  2:831439/1‑150 (MQ=255)
                                                                                                                   tAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGaaa                                  >  2:1331238/1‑150 (MQ=255)
                                                                                                                     aTTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATc                                <  2:52022/150‑1 (MQ=255)
                                                                                                                      ttGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCa                               <  2:1286430/150‑1 (MQ=255)
                                                                                                                       tGAGGGATAACAACAAAAAAGGTGAGTTTTGCGCCTCACATTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGATTGAGTTCGCAGCCATTAGCGATTAAATCAt                              >  1:1316270/1‑150 (MQ=255)
                                                                                                                       tGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTcc                                                                   <  2:714129/113‑1 (MQ=255)
                                                                                                                       tGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTcc                                                                   >  1:714129/1‑113 (MQ=255)
                                                                                                                         aGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGa                            >  1:477288/1‑150 (MQ=255)
                                                                                                                            gATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCAGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGc                         <  2:722649/150‑1 (MQ=255)
                                                                                                                              tAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCTg                       >  2:1425745/1‑150 (MQ=255)
                                                                                                                              tAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCTg                       <  1:407906/150‑1 (MQ=255)
                                                                                                                                caaataaaaaaGGTGAGTCGCAAAACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCTGGa                     >  1:800198/6‑150 (MQ=38)
                                                                                                                                acaacaAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAACTCATCCTGAATCGCTGGTGCGTTCGGATTTTCCTGCGGTTTGAGTGCGCCGCCATTAGCGATGAAATCATGACGCTGGa                     <  2:1527158/150‑1 (MQ=255)
                                                                                                                                acaacaAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCTGGa                     >  1:651145/1‑150 (MQ=255)
                                                                                                                                 caacaaAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCTGGaa                    <  1:388401/150‑1 (MQ=255)
                                                                                                                                 aaataaaaaaGGTGAGTCGCAAAACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCTGGaa                    >  2:1363775/5‑150 (MQ=39)
                                                                                                                                 aaataaaaaaGGTGAGTCGCAAAACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCTGGaa                    >  2:1168461/5‑150 (MQ=39)
                                                                                                                                  aataaaaaaGGTGAGTCGCAAAACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCTGGAAg                   >  1:10108/4‑150 (MQ=255)
                                                                                                                                  aataaaaaaGGTGAGTCGCAAAACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCTGGAAg                   >  1:3945/4‑150 (MQ=255)
                                                                                                                                  aacaaAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCTGGAAg                   <  2:1145639/150‑1 (MQ=255)
                                                                                                                                  aacaaAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCACAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTACCAGCGGTTTGAGTTCGCCGCCATTAGAGATGAAATCATGACGCTGGAAg                   >  1:1473326/1‑150 (MQ=255)
                                                                                                                                      aaaaaGGTGCGTTTTGCGACTCACCTTTTTTATTTGTTTGTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGACATCATGACGCTGGAAGTACg               >  1:64350/1‑150 (MQ=255)
                                                                                                                                      aaaaaGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCTGGAAGTACg               >  1:769211/1‑150 (MQ=255)
                                                                                                                                        aaaGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGAGTTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCTGGAAGTACGCt             >  2:620874/1‑150 (MQ=255)
                                                                                                                                          aGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCTGGAAGTACGCTTc           <  1:982028/150‑1 (MQ=255)
                                                                                                                                            gTGAGTTTTGGGACTCACCTTTTTTATTTGTTTATTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCTGGAAGTACGCTTcgc         <  2:242265/150‑1 (MQ=255)
                                                                                                                                            gTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCTGGAAGTACGCTTcgc         <  1:1130563/150‑1 (MQ=255)
                                                                                                                                             tGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCTGGAAGTACGCTTcgcg        <  1:946935/150‑1 (MQ=255)
                                                                                                                                                 ttttGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCTGGAAGTACGCTTCGCGCAcc    <  1:1374413/150‑1 (MQ=255)
                                                                                                                                                 ttttGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCTGGAAGTACGCTTCGCGCAcc    <  1:379693/150‑1 (MQ=255)
                                                                                                                                                 ttttGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGGTTGAGTTCGCCGCCATTAGCGAGGAAATCATGACGCTGGAAGTACGCTTCGCGcccc    <  1:379690/150‑1 (MQ=255)
                                                                                                                                                   ttGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCTGGAAGTACGCTTCGCGCACCat  >  2:721021/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TTTCCAGTTTTTTGATTTTTCAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCTGGAAGTACGCTTCGCGCACCAT  >  NC_000913/2464078‑2464374

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