Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 752,068 G→C 6.4% R34R (CGC→CGG ybgQ ← putative fimbrial usher protein YbgQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913752,0680GC 6.4% 373.9 / 12.6 110R34R (CGC→CGGybgQputative fimbrial usher protein YbgQ
Reads supporting (aligned to +/- strand):  ref base G (46/57);  new base C (3/4);  total (49/61)
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.61e-01

Only 100 of 113 total aligned reads displayed.

TGATGCATGGAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGATGATG  >  NC_000913/751920‑752213
                                                                                                                                                    |                                                                                                                                                 
  aTGCATGGAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCgcgc                                                                                                                                                <  2:188872/150‑1 (MQ=255)
     cATGGAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCAtt                                                                                                                                             >  2:1410425/1‑150 (MQ=255)
         gAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGAtt                                                                                                                                         >  1:1176365/1‑150 (MQ=255)
         gAATGGTTTTCTCACCCTTTTTTTGCCAGTAAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACTCTAACAAAATATTCACCGGGATCAATTACTCCTTTTTCCTTTAATAATTAAATATCAATGCGCGACCGCATTGAta                                                                                                                                         <  2:1044750/150‑2 (MQ=255)
           aTGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGCGACCGCATTGAtatt                                                                                                                                       <  2:941759/150‑4 (MQ=255)
            tGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGCGACCGCATTGAtattt                                                                                                                                      <  2:115239/150‑5 (MQ=255)
             ggTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATc                                                                                                                                     <  1:1279038/150‑1 (MQ=255)
              gTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCg                                                                                                                                    <  2:758379/150‑1 (MQ=255)
              gTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCg                                                                                                                                    <  1:125675/150‑1 (MQ=255)
                 ttGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGTGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCgaga                                                                                                                                 <  2:1539494/150‑1 (MQ=255)
                 ttGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCgaga                                                                                                                                 <  1:501522/150‑1 (MQ=255)
                  tGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGaa                                                                                                                                >  1:966415/1‑150 (MQ=255)
                  tGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATCCGGTCGCGCATTGATTTATCGAGaa                                                                                                                                <  1:626136/150‑1 (MQ=255)
                     cACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACAt                                                                                                                             >  1:743465/1‑150 (MQ=255)
                      aCCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCga                                                                                                                                   <  2:1051554/143‑1 (MQ=255)
                      aCCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCga                                                                                                                                   >  1:1051554/1‑143 (MQ=255)
                       cccTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTc                                                                                                                           >  1:1112164/1‑150 (MQ=255)
                         cTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAg                                                                                                                         <  2:743465/150‑1 (MQ=255)
                          tttttttGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGa                                                                                                                        >  1:138164/1‑150 (MQ=255)
                           ttttttGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGAt                                                                                                                       >  2:745523/1‑150 (MQ=255)
                            tttttGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGAtt                                                                                                                      <  2:1663137/150‑1 (MQ=255)
                                tGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTc                                                                                                                           <  1:978133/141‑1 (MQ=255)
                                 gCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGGCTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAAtt                                                                                                                 >  1:1555558/1‑150 (MQ=255)
                                  ccAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTc                                                                                                                <  1:342375/150‑1 (MQ=255)
                                    aGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGa                                                                                                              <  2:1003485/150‑1 (MQ=255)
                                     gTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGAc                                                                                                             <  2:874870/150‑1 (MQ=255)
                                         aTTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGcc                                                                                                         >  2:216536/1‑150 (MQ=255)
                                         aTTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGcc                                                                                                         >  2:795256/1‑150 (MQ=255)
                                         aTTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGcc                                                                                                         >  1:121772/1‑150 (MQ=255)
                                                gCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTgcgc                                                                                                  >  2:416239/1‑150 (MQ=255)
                                                  ccATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCaa                                                                                                >  2:1266826/1‑150 (MQ=255)
                                                   cATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAg                                                                                               <  1:397773/150‑1 (MQ=255)
                                                    aTTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAgg                                                                                              <  1:829095/150‑1 (MQ=255)
                                                        cTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATc                                                                                          >  2:1583766/1‑150 (MQ=255)
                                                            ttttGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCa                                                                                      <  1:304332/150‑1 (MQ=255)
                                                               tgttgtTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGa                                                                                   >  1:615561/1‑150 (MQ=255)
                                                                gttgttATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGAc                                                                                  <  1:108924/150‑1 (MQ=255)
                                                                gttgttATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGAc                                                                                  <  1:98736/150‑1 (MQ=255)
                                                                 ttgttATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGAcc                                                                                 >  1:563218/1‑150 (MQ=255)
                                                                 ttgttATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGAcc                                                                                 <  1:1300558/150‑1 (MQ=255)
                                                                  tgttATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACca                                                                                >  1:36078/1‑150 (MQ=255)
                                                                    ttATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACcagg                                                                              <  1:713091/150‑1 (MQ=255)
                                                                    ttATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACcagg                                                                              <  1:1410425/150‑1 (MQ=255)
                                                                     tATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACcaggc                                                                             >  1:458267/1‑150 (MQ=255)
                                                                     tATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACcaggc                                                                             >  1:1657679/1‑150 (MQ=255)
                                                                      aTTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACcaggca                                                                            >  2:325415/1‑150 (MQ=255)
                                                                      aTTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACcaggca                                                                            >  2:315315/1‑150 (MQ=255)
                                                                      aTTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACcaggca                                                                            >  2:1123382/1‑150 (MQ=255)
                                                                       ttCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACcaggcag                                                                           >  1:914063/1‑150 (MQ=255)
                                                                       ttCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTAACGAGAACAATCACAATGAATTCGACCCCCAATGCGCAAGCCAACGCCATGACcaggcag                                                                           >  1:908133/1‑150 (MQ=255)
                                                                        tCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACcaggcagg                                                                          <  2:264689/150‑1 (MQ=255)
                                                                        tCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACcaggcagg                                                                          >  1:785382/1‑150 (MQ=255)
                                                                          aCCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGTAGGAt                                                                        <  1:500161/150‑1 (MQ=255)
                                                                          aCCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGAt                                                                        <  2:1176365/150‑1 (MQ=255)
                                                                            ggCAACGCTAACAAAATATTCAGCGGGAGCAATGACTCCTTTTTGCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATAc                                                                      <  2:740396/149‑1 (MQ=255)
                                                                             gCAAAGCTAACAAAATATTCACCGGTAGCAATGAGTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACa                                                                     <  1:261106/150‑1 (MQ=255)
                                                                                   cTAACAAAATATTGACAGGTAGCAATGACTGATTTTTCCTTTAATAATGAAATATCAATGCGGTCGAGCATTGATTTATCGAGATCATTCAGATTGAATTCTACCGGCATTGCGCAAGGCTTCGCCATGACCAGGCAGGATACAAAagag                                                               <  2:1363987/150‑1 (MQ=255)
                                                                                   cTAACAAAATATTCACGGGGAGCAATGAGTCCTTTTTCCTTTTATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGATCATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAagag                                                               <  2:891370/150‑1 (MQ=255)
                                                                                   cTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAagag                                                               <  1:416239/150‑1 (MQ=255)
                                                                                   cTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAagag                                                               <  2:388291/150‑1 (MQ=255)
                                                                                    tAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAagaga                                                              <  1:745523/150‑1 (MQ=255)
                                                                                     aaCAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAagagag                                                             <  1:896503/150‑1 (MQ=255)
                                                                                      aCAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGt                                                            >  1:220761/1‑150 (MQ=255)
                                                                                       cAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACcaggc                                                                             >  2:712939/1‑132 (MQ=255)
                                                                                       cAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACcaggc                                                                             <  1:712939/132‑1 (MQ=255)
                                                                                       aaaaaTATTCACCGGGATCAATGACTCCTTTTTCGTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGATAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCTCCATGACCAGGCAGGATACAAAAGAGAGTc                                                           <  2:136092/149‑1 (MQ=255)
                                                                                         aaaTATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGa                                                         >  1:1578188/1‑150 (MQ=255)
                                                                                         aaaTATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGa                                                         >  2:728957/1‑150 (MQ=255)
                                                                                            tatTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATaa                                                      <  1:542947/150‑1 (MQ=255)
                                                                                                cACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAatat                                                  >  1:653438/1‑150 (MQ=255)
                                                                                                    gggAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCAc                                              <  1:1266826/150‑1 (MQ=255)
                                                                                                    gggAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCAc                                              >  1:1087461/1‑150 (MQ=255)
                                                                                                     ggAGCAATGACTCCTTTTTCGTTTAATAATGAAATATCAATGCGGTCGCTCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTGACg                                             <  1:1248802/150‑1 (MQ=255)
                                                                                                     ggAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACg                                             <  1:6833/150‑1 (MQ=255)
                                                                                                      gAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACgg                                            <  2:615561/150‑1 (MQ=255)
                                                                                                       aGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAagaga                                                              <  2:277186/131‑1 (MQ=255)
                                                                                                        gCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGgtg                                          >  1:270929/1‑150 (MQ=255)
                                                                                                           aTGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTcc                                       <  2:772509/150‑1 (MQ=255)
                                                                                                           aTGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTcc                                       >  1:887856/1‑150 (MQ=255)
                                                                                                            tGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCa                                      <  2:835922/150‑1 (MQ=255)
                                                                                                              aCTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATa                                    <  1:325415/150‑1 (MQ=255)
                                                                                                               cTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATAc                                   >  1:738334/1‑150 (MQ=255)
                                                                                                                 ccTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCt                                 <  1:301744/150‑1 (MQ=255)
                                                                                                                 ccTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCt                                 >  2:1044310/1‑150 (MQ=255)
                                                                                                                     tttCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATa                             <  2:785382/150‑1 (MQ=255)
                                                                                                                       tCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGATTTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATaaa                           >  1:1040398/1‑150 (MQ=255)
                                                                                                                          tttAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAatat                        <  1:1245588/150‑1 (MQ=255)
                                                                                                                            taataaTGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATAttt                      >  2:171851/1‑150 (MQ=255)
                                                                                                                            taataaTGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATAttt                      >  2:166111/1‑150 (MQ=255)
                                                                                                                             aataatGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATAtttt                     <  1:457561/150‑1 (MQ=255)
                                                                                                                             aataatGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATAtttt                     <  2:1224549/150‑1 (MQ=255)
                                                                                                                              ataatGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTa                    <  2:966415/150‑1 (MQ=255)
                                                                                                                               taatGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTAt                   <  1:932824/150‑1 (MQ=255)
                                                                                                                                 atGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGa                 >  1:212703/1‑150 (MQ=255)
                                                                                                                                 atGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGa                 >  2:175177/1‑150 (MQ=255)
                                                                                                                                   gataaatCAATGCGCGACCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGaaa               >  1:1178232/6‑150 (MQ=255)
                                                                                                                                           aaTGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATCCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCga       <  2:1657679/150‑1 (MQ=255)
                                                                                                                                              gCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCgatga    >  1:740533/1‑150 (MQ=255)
                                                                                                                                                ggTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGatgatg  <  1:795256/150‑1 (MQ=255)
                                                                                                                                                ggTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGatgatg  <  1:216536/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TGATGCATGGAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGATGATG  >  NC_000913/751920‑752213

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