Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,202,274 G→A 6.6% Y220Y (TAC→TAT ymfK ← e14 prophage; putative repressor protein YmfK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,202,2740GA 6.6% 284.4 / 6.8 91Y220Y (TAC→TATymfKe14 prophage; putative repressor protein YmfK
Reads supporting (aligned to +/- strand):  ref base G (41/44);  new base A (3/3);  total (44/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

GGTCTTCGTACGTACCGGGCGGGCTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAGGTGGTT  >  NC_000913/1202125‑1202418
                                                                                                                                                     |                                                                                                                                                
ggTCTTCGTACGTACCGGGCGGGCTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTg                                                                                                                                                  <  1:1201899/150‑1 (MQ=255)
  tCTTCGTACGTACCGGGCGGGCTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTa                                                                                                                                                <  2:1314678/150‑1 (MQ=255)
   cTTCGTACGTACCGGGCGGGCTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTAt                                                                                                                                               <  1:1052641/150‑1 (MQ=255)
       gtacgtacCGGGCGGGCTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTAtcttc                                                                                                                                           >  2:832131/1‑150 (MQ=255)
           gtacCGGGCGGGCTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTgg                                                                                                                                       <  2:340263/150‑1 (MQ=255)
           gtacCGGGCGGGCTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTgg                                                                                                                                       >  2:895766/1‑150 (MQ=255)
                  gcgggcTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCt                                                                                                                                <  2:293084/150‑1 (MQ=255)
                  gcgggcTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCt                                                                                                                                <  1:291448/150‑1 (MQ=255)
                   cgggcTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTg                                                                                                                               <  1:660817/150‑1 (MQ=255)
                   cgggcTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTg                                                                                                                               <  1:110425/150‑1 (MQ=255)
                    gggcTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTga                                                                                                                              >  2:1403456/1‑150 (MQ=255)
                    gggcTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTga                                                                                                                              >  2:1409170/1‑150 (MQ=255)
                     ggcTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTgag                                                                                                                             <  2:1598379/150‑1 (MQ=255)
                     ggcTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTgag                                                                                                                             <  1:1352254/150‑1 (MQ=255)
                           ttCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATc                                                                                                                       >  2:1358534/1‑150 (MQ=255)
                           ttCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATc                                                                                                                       >  1:1569415/1‑150 (MQ=255)
                             cGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATcac                                                                                                                     <  1:1174663/150‑1 (MQ=255)
                                 ggCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTa                                                                                                                 >  2:1269826/1‑150 (MQ=255)
                                 ggCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTa                                                                                                                 <  1:358657/150‑1 (MQ=255)
                                 ggCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTa                                                                                                                 >  1:1000589/1‑150 (MQ=255)
                                  gCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTAc                                                                                                                >  1:1052406/1‑150 (MQ=255)
                                      cctgcctgTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAAt                                                                                                            >  1:1198033/1‑150 (MQ=255)
                                        tgcctgTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTATATATTCTTGGTTTTCCTGAGAAAATCACAGTACCAATta                                                                                                          >  1:165652/1‑150 (MQ=255)
                                         gcctgTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATtat                                                                                                         >  2:1589612/1‑150 (MQ=255)
                                           ctgTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATag                                                                                                       >  2:1619682/1‑150 (MQ=255)
                                                tgttgtTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCaa                                                                                                  >  2:1194917/1‑150 (MQ=255)
                                                tgttgtTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCaa                                                                                                  <  2:703799/150‑1 (MQ=255)
                                                    gttTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTAc                                                                                              >  2:949413/1‑150 (MQ=255)
                                                    gttTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTCATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTAc                                                                                              >  2:1403557/1‑150 (MQ=255)
                                                     tttCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTAcc                                                                                             >  2:864142/1‑150 (MQ=255)
                                                        ctctTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGtt                                                                                          <  2:45425/150‑1 (MQ=255)
                                                          ctTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGa                                                                                        >  1:1289302/1‑150 (MQ=255)
                                                            tGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATc                                                                                      <  2:1162876/150‑1 (MQ=255)
                                                             gggTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCt                                                                                     >  2:781954/1‑150 (MQ=255)
                                                               gTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTa                                                                                   <  1:832131/150‑1 (MQ=255)
                                                                tACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTaa                                                                                  <  1:584558/150‑1 (MQ=255)
                                                                tACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTaa                                                                                  <  1:296208/150‑1 (MQ=255)
                                                                  cATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATg                                                                                <  2:546255/150‑1 (MQ=255)
                                                                    ttATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTa                                                                              <  2:240235/150‑1 (MQ=255)
                                                                         tATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCt                                                                         >  1:1153395/1‑150 (MQ=255)
                                                                          aTCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTATATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTc                                                                        <  2:165652/150‑1 (MQ=255)
                                                                            ctcAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGATCAATTACCGTTTATCTTAATTTAAGGATCAg                                                                      >  2:1486906/1‑150 (MQ=255)
                                                                            ctcAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAg                                                                      <  2:1293474/150‑1 (MQ=255)
                                                                             tcAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAgg                                                                     <  1:781954/150‑1 (MQ=255)
                                                                             tcAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAgg                                                                     >  1:1486677/1‑150 (MQ=255)
                                                                                aaGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCa                                                                  >  1:40756/1‑150 (MQ=255)
                                                                                aaGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCa                                                                  >  1:629219/1‑150 (MQ=255)
                                                                                  ggTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGt                                                                >  2:814441/1‑150 (MQ=255)
                                                                                     aCATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTg                                                             <  2:1617019/150‑1 (MQ=255)
                                                                                      aaTTGTCAAGTCTAAAAAAACCTTCCGAAGCAGGTTCATAAACATTGATTAGGATTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCGATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTgg                                                            <  2:1384003/149‑1 (MQ=255)
                                                                                      aaTTGTCAAGAATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCGATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTgg                                                            <  2:1389721/149‑1 (MQ=255)
                                                                                       aTTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTggg                                                           >  2:1198943/1‑150 (MQ=255)
                                                                                       aTTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTggg                                                           <  1:1148345/150‑1 (MQ=255)
                                                                                       aTTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTggg                                                           >  2:1490359/1‑150 (MQ=255)
                                                                                       aTTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTggg                                                           >  1:136380/1‑150 (MQ=255)
                                                                                        ttGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGt                                                          <  2:1371982/150‑1 (MQ=255)
                                                                                            cAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTaa                                                      >  2:298807/1‑150 (MQ=255)
                                                                                            cAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTATATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTaa                                                                             <  2:1065826/127‑1 (MQ=255)
                                                                                            cAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTATATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTaa                                                                             >  1:1065826/1‑127 (MQ=255)
                                                                                            cAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTATATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTaa                                                                             >  1:1071341/1‑127 (MQ=255)
                                                                                            cAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTATATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTaa                                                                             <  2:1071341/127‑1 (MQ=255)
                                                                                              aGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACg                                                    >  1:205058/1‑150 (MQ=255)
                                                                                               gTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTAGCTTCTTGGTTTTCCTGAGAAAATCCCACTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAAAGc                                                   >  2:1027482/1‑150 (MQ=255)
                                                                                                 ataAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCtt                                                 >  2:221487/1‑150 (MQ=255)
                                                                                                 ataAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAAATATAGAGCAATTACCGTTGAACTTAATGTAAGGCTCAGGCCAGTTTGGGTTAAACGCtt                                                 >  2:17788/1‑150 (MQ=255)
                                                                                                   aaaaaaaCCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTg                                               >  1:1488262/1‑150 (MQ=255)
                                                                                                    aaaaaaCCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTga                                              <  1:1293331/150‑1 (MQ=255)
                                                                                                    aaaaaaCCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTga                                              <  2:210595/150‑1 (MQ=255)
                                                                                                     aaaaaCCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTgag                                             <  1:1358534/150‑1 (MQ=255)
                                                                                                              ccGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCtg                                    >  2:203625/1‑150 (MQ=255)
                                                                                                                  aGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTc                                >  2:1576862/1‑150 (MQ=255)
                                                                                                                  aGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTc                                <  2:1052406/150‑1 (MQ=255)
                                                                                                                    cAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTccc                              >  1:1038902/1‑150 (MQ=255)
                                                                                                                     aGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCa                             <  1:1101702/150‑1 (MQ=255)
                                                                                                                      ggTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCAt                            <  1:1269826/150‑1 (MQ=255)
                                                                                                                          cATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCAtcttc                        <  1:895766/150‑1 (MQ=255)
                                                                                                                          cATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCAtcttc                        <  1:575969/150‑1 (MQ=255)
                                                                                                                             aaaCATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTAt                     <  1:1619682/150‑1 (MQ=255)
                                                                                                                                cATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCaa                  >  2:1625452/1‑150 (MQ=255)
                                                                                                                                  ttGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAAcc                <  1:592758/150‑1 (MQ=255)
                                                                                                                                      ttAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCtttt            <  1:1589612/150‑1 (MQ=255)
                                                                                                                                        aGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATa                                          >  2:1409365/1‑118 (MQ=255)
                                                                                                                                        aGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATa                                          <  1:1409365/118‑1 (MQ=255)
                                                                                                                                          gCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAg        <  1:1512143/150‑1 (MQ=255)
                                                                                                                                           cTTTGATTTTGTATCTTCTTGGTTTTCCTGATAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATATCGCTGTGTCCCATCTTCTATCAACCTTTTGAAgg       <  2:1297656/150‑1 (MQ=255)
                                                                                                                                           cTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAgg       >  2:1271086/1‑150 (MQ=255)
                                                                                                                                            tttGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAggt      <  2:1038902/150‑1 (MQ=255)
                                                                                                                                            tttGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAggt      >  2:988609/1‑150 (MQ=255)
                                                                                                                                              tGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAggtgg    >  1:222348/1‑150 (MQ=255)
                                                                                                                                                aTTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAGGTGGtt  <  1:814441/150‑1 (MQ=255)
                                                                                                                                                aTTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAGGTGGtt  <  2:1289302/150‑1 (MQ=255)
                                                                                                                                                aTTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTAAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAATGTGGtt  <  1:864142/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GGTCTTCGTACGTACCGGGCGGGCTACTTCGTGGGCGTCCTGCCTGTTTGTTGTTTCTCTTGGGTACATTATGTATCTCAAAGGTACATTGTCAAGTATAAAAAAACCTGCCGAAGCAGGTTCATAAACATTGATTAGGCTTTGATTTTGTATCTTCTTGGTTTTCCTGAGAAAATCACAGTACCAATTATAGAGCAATTACCGTTGATCTTAATGTAAGGCTCAGGCCAGTTTGGGTTTAACGCTTTGAGATAACGCTGTGTCCCATCTTCTATCAACCTTTTGAAGGTGGTT  >  NC_000913/1202125‑1202418

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