Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,458,420 T→A 8.0% Y28N (TAT→AAT)  gspE → Type II secretion system protein GspE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,458,4200TA 8.0% 371.1 / 16.9 125Y28N (TAT→AAT) gspEType II secretion system protein GspE
Reads supporting (aligned to +/- strand):  ref base T (58/57);  new base A (5/5);  total (63/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

Only 100 of 130 total aligned reads displayed.

TGCCGGTGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGAGTCGCAACAGATTGC  >  NC_000913/3458271‑3458568
                                                                                                                                                     |                                                                                                                                                    
tggcGGTGCTGGATGAAAACACGTTCAACAGTCACGGCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGAttattt                                                                                                                                                      <  1:1498133/147‑1 (MQ=255)
tGCCGGTGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGAttattt                                                                                                                                                      <  2:333296/150‑1 (MQ=255)
 gCCGGTGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGAttattta                                                                                                                                                     <  2:1443185/150‑1 (MQ=255)
      tGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCg                                                                                                                                                >  1:1621543/1‑150 (MQ=255)
       gCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGc                                                                                                                                               >  1:944756/1‑150 (MQ=255)
        cTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGATTTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGcc                                                                                                                                              <  2:1158612/150‑1 (MQ=255)
        cTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGcc                                                                                                                                              <  1:1417168/150‑1 (MQ=255)
         tGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCg                                                                                                                                             <  2:149845/150‑1 (MQ=255)
         tGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCg                                                                                                                                             <  2:837218/150‑1 (MQ=255)
         tGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATAAATAATCTCGccc                                                                                                                                             <  1:1056301/150‑2 (MQ=255)
         tGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATAAATAATCTCGccc                                                                                                                                             <  2:1073253/150‑2 (MQ=255)
         tGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATAAATAATCTCGccc                                                                                                                                             <  1:972252/150‑2 (MQ=255)
           gATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGAt                                                                                                                                           <  2:1621543/150‑1 (MQ=255)
             tGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATAc                                                                                                                                         <  1:356166/150‑1 (MQ=255)
             tGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATAc                                                                                                                                         >  2:350659/1‑150 (MQ=255)
             tGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATAAATAATCTCGccctctg                                                                                                                                         <  2:403848/150‑6 (MQ=255)
              gAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACg                                                                                                                                        <  2:199701/150‑1 (MQ=255)
               aaaaCACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGc                                                                                                                                       >  1:1285148/1‑150 (MQ=255)
                    acGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTc                                                                                                                                  <  1:268122/150‑1 (MQ=255)
                     cGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGTTACGCCATTCg                                                                                                                                 <  2:944756/150‑1 (MQ=255)
                       ttCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAg                                                                                                                               <  2:149901/150‑1 (MQ=255)
                       ttCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAg                                                                                                                               <  1:415700/150‑1 (MQ=255)
                         cAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCg                                                                                                                             >  1:173686/1‑150 (MQ=255)
                         cAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCg                                                                                                                             >  2:1597009/1‑150 (MQ=255)
                           aCAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGt                                                                                                                           >  2:1320962/1‑150 (MQ=255)
                                cACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCCCTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTc                                                                                                                      >  1:606991/1‑150 (MQ=255)
                                cACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTc                                                                                                                      <  1:1439977/150‑1 (MQ=255)
                                    ccccTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGata                                                                                                                  <  2:483298/150‑1 (MQ=255)
                                        tGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCa                                                                                                              >  2:807674/1‑150 (MQ=255)
                                         gcgcCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCaa                                                                                                             <  1:399124/150‑1 (MQ=255)
                                          cgcCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAAc                                                                                                            <  1:546583/150‑1 (MQ=255)
                                            ccATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATATCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGt                                                                                                          >  2:1515537/1‑150 (MQ=255)
                                             cATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTc                                                                                                         <  2:909388/150‑1 (MQ=255)
                                             cATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTc                                                                                                         >  2:233/1‑150 (MQ=255)
                                               tCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAg                                                                                                       <  1:1096375/150‑1 (MQ=255)
                                                 gTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGt                                                                                                     >  1:988030/1‑150 (MQ=255)
                                                  tCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGtt                                                                                                    <  2:760683/150‑1 (MQ=255)
                                                     cGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGc                                                                                                 <  2:552251/150‑1 (MQ=255)
                                                      ggTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGcc                                                                                                <  2:1477452/150‑1 (MQ=255)
                                                      ggTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGcc                                                                                                <  2:52074/150‑1 (MQ=255)
                                                         gAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGt                                                                                             <  2:468438/150‑1 (MQ=255)
                                                           ggCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGc                                                                                           >  1:851287/1‑150 (MQ=255)
                                                           ggCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGc                                                                                           >  2:1322048/1‑150 (MQ=255)
                                                            gCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGcc                                                                                          <  1:1320962/150‑1 (MQ=255)
                                                              aTTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAg                                                                                        <  2:1285148/150‑1 (MQ=255)
                                                              aTTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAg                                                                                        <  2:527319/150‑1 (MQ=255)
                                                               ttCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGa                                                                                       <  1:1597009/150‑1 (MQ=255)
                                                                  atatGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCa                                                                                    >  1:656685/1‑150 (MQ=255)
                                                                   tatGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCAt                                                                                   >  2:396990/1‑150 (MQ=255)
                                                                   tatGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCAt                                                                                   <  2:1573420/150‑1 (MQ=255)
                                                                     tGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTc                                                                                                         <  1:1603826/126‑1 (MQ=255)
                                                                     tGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTc                                                                                                         >  2:1603826/1‑126 (MQ=255)
                                                                       agaATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGAcc                                                                               <  1:156792/150‑1 (MQ=255)
                                                                        gaATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCa                                                                              >  1:872459/1‑150 (MQ=255)
                                                                            tcactcacCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCtt                                                                          <  2:1257759/150‑1 (MQ=255)
                                                                             cactcacCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTg                                                                         >  1:896813/1‑150 (MQ=255)
                                                                             cactcacCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTg                                                                         >  2:1328100/1‑150 (MQ=255)
                                                                               ctcacCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTc                                                                       >  2:51721/1‑150 (MQ=255)
                                                                                     gTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGc                                                                 <  2:70877/150‑1 (MQ=255)
                                                                                     gTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGc                                                                 >  2:953583/1‑150 (MQ=255)
                                                                                      tACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCt                                                                <  1:920661/150‑1 (MQ=255)
                                                                                       aCCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTg                                                               <  1:1391876/150‑1 (MQ=255)
                                                                                          ccGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGAtt                                                            >  2:884511/1‑150 (MQ=255)
                                                                                           cGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGAttt                                                           >  2:763136/1‑150 (MQ=255)
                                                                                            gCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGAtttt                                                          >  1:678760/1‑150 (MQ=255)
                                                                                              cAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGa                                                        <  2:173686/150‑1 (MQ=255)
                                                                                                  tGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCt                                                    >  2:1094010/1‑150 (MQ=255)
                                                                                                    ggCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCg                                                  >  2:107777/1‑150 (MQ=255)
                                                                                                    ggCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCg                                                  >  2:62578/1‑150 (MQ=255)
                                                                                                      cgcTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGc                                                >  1:1151529/1‑150 (MQ=255)
                                                                                                      cgcTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGc                                                >  2:864320/1‑150 (MQ=255)
                                                                                                         tGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTgg                                             <  2:111249/150‑1 (MQ=255)
                                                                                                           gCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGaa                                           >  1:876474/1‑150 (MQ=255)
                                                                                                            cacaACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAg                                          >  1:368283/1‑150 (MQ=255)
                                                                                                             acaACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGc                                         >  1:1537820/1‑150 (MQ=255)
                                                                                                             acaACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGc                                         <  2:1240913/150‑1 (MQ=255)
                                                                                                              caACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCg                                        >  1:497565/1‑150 (MQ=255)
                                                                                                                 cGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTa                                     >  2:316568/1‑150 (MQ=255)
                                                                                                                 cGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTa                                     >  1:1446839/1‑150 (MQ=255)
                                                                                                                    aTTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTc                                  <  2:1438574/150‑1 (MQ=255)
                                                                                                                       ggTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCAt                               <  1:273428/150‑1 (MQ=255)
                                                                                                                       ggTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCAt                               >  1:31038/1‑150 (MQ=255)
                                                                                                                        gTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATc                              <  2:497565/150‑1 (MQ=255)
                                                                                                                          tATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGACCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAg                            <  2:988030/150‑1 (MQ=255)
                                                                                                                           aTCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGa                           >  2:320041/1‑150 (MQ=255)
                                                                                                                            tctcTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGaa                          >  2:419251/1‑150 (MQ=255)
                                                                                                                            tctcTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGaa                          <  2:939403/150‑1 (MQ=255)
                                                                                                                             ctctATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAAt                         >  2:1115610/1‑150 (MQ=255)
                                                                                                                                tATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATAcc                      >  1:1195523/1‑150 (MQ=255)
                                                                                                                                tATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATAcc                      <  1:107777/150‑1 (MQ=255)
                                                                                                                                  ttCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCgg                    >  1:1534181/1‑150 (MQ=255)
                                                                                                                                       agGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGAGt               <  1:953583/150‑1 (MQ=255)
                                                                                                                                        gggCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCgg                                       <  1:1105692/125‑1 (MQ=255)
                                                                                                                                         ggCGAGATAAATAATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGAGTCg             >  2:432243/1‑150 (MQ=255)
                                                                                                                                         ggCGAGATAAATAATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGAGTCg             >  1:91582/1‑150 (MQ=255)
                                                                                                                                         ggCGAGATAAATAATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGAATTCCATCAGAATACCGGTGAGTCg             >  2:836207/1‑150 (MQ=255)
                                                                                                                                             agaTTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGAGTCGCAAc         <  1:233/150‑1 (MQ=255)
                                                                                                                                             agaTTATTTATCTCGCCGATACGCCATTCGAGCGGGTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGAGTCGCAAc         <  2:606991/150‑1 (MQ=255)
                                                                                                                                              gaTTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGAGTCGCAACa        <  2:1279310/150‑1 (MQ=255)
                                                                                                                                              gaTTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGAGTCGCAACa        >  1:1501430/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGCCGGTGCTGGATGAAAACACGTTCAACAGTCACGCCCCTGCGCCATCGTCACGGTGAGGCATTCATATGAGAATTCACTCACCGTACCCCGCCAGTTGGGCGCTGGCACAACGAATTGGTTATCTCTATTCAGAGGGCGAGATTATTTATCTCGCCGATACGCCATTCGAGCGGTTACTCGATATTCAACGTCAGGTTGGCCAGTGCCAGACCATGACCAGCTTGTCACAGGCTGATTTTGAAGCTCGGCTGGAAGCGGTATTCCATCAGAATACCGGTGAGTCGCAACAGATTGC  >  NC_000913/3458271‑3458568

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