Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,530,076 (A)7→8 coding (155/483 nt) ydcD → putative immunity protein for RhsE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,530,0691.A94.3% 306.4 / 10.5 88coding (148/483 nt)ydcDputative immunity protein for RhsE
Reads supporting (aligned to +/- strand):  ref base . (2/3);  new base A (48/35);  total (50/38)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

GTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGC  >  NC_000913/1529971‑1530162
                                                                                                   |                                                                                              
gTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGa                                                                                               <  2:1850568/100‑1 (MQ=255)
gTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGa                                                                                               >  1:1089458/1‑100 (MQ=255)
  ttATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaaa                                                                                             >  1:1071161/1‑100 (MQ=255)
  ttATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaaa                                                                                             <  2:1116296/100‑1 (MQ=255)
   tATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaaaa                                                                                            >  1:1067589/1‑100 (MQ=255)
    aTTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaaaaa                                                                                           >  2:1662146/1‑100 (MQ=255)
       tgttgtTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaaaaaaaa                                                                                       <  1:1174993/100‑2 (MQ=255)
        gttgttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaaaaaaaat                                                                                      <  2:1613288/100‑3 (MQ=255)
        gttgttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaaaaaaaat                                                                                      <  1:886604/100‑3 (MQ=255)
          tgttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaaaaaaaatat                                                                                    <  1:1698890/100‑5 (MQ=255)
           gttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaaaaaaaatatt                                                                                   >  2:1786914/1‑95 (MQ=255)
           gttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaaaaaaaatatt                                                                                   <  2:1718860/100‑6 (MQ=255)
           gttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaaaaaaaatatt                                                                                   >  1:42110/1‑95 (MQ=255)
              tACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATattatt                                                                                 >  1:952915/1‑100 (MQ=255)
               aCCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATAttattt                                                                                >  2:1714988/1‑100 (MQ=255)
               aCCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATAttattt                                                                                >  2:1715070/1‑100 (MQ=255)
                ccTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATAttattta                                                                               >  1:1644483/1‑100 (MQ=255)
                ccTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATAttattt                                                                                >  2:882215/1‑99 (MQ=255)
                 cgatcttttGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATAttatttat                                                                              <  1:760689/94‑1 (MQ=255)
                 cTTGGGTTTGATAAAAATTCAAAATATGACATCGTATTTTATAGACATAAAAAACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATAttatttat                                                                              <  2:2037596/100‑1 (MQ=255)
                   tGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATaaa                                                                           <  1:1481729/100‑1 (MQ=255)
                      gTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAaat                                                                         <  2:1943834/100‑1 (MQ=255)
                       tttGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTtatagatcg                                                                        >  2:760686/1‑95 (MQ=255)
                         tGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAaataat                                                                      <  1:1677796/100‑1 (MQ=255)
                          gATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaatt                                                                     >  2:243332/1‑100 (MQ=255)
                          gATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaatt                                                                     <  1:993179/100‑1 (MQ=255)
                          gATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaatt                                                                     >  2:115154/1‑100 (MQ=255)
                          gATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaatt                                                                     >  1:1188389/1‑100 (MQ=255)
                            tAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaattaa                                                                   >  2:1523369/1‑100 (MQ=255)
                            tAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaattaa                                                                   >  2:676144/1‑100 (MQ=255)
                            tAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaattaa                                                                   >  2:560978/1‑100 (MQ=255)
                            tAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaattaa                                                                   <  1:1956748/100‑1 (MQ=255)
                            tAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaattaa                                                                   >  2:1980956/1‑100 (MQ=255)
                             aaCAATTCAAATTCTGACATCGTATTTTATGGAGATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaattaat                                                                  >  2:1349923/1‑100 (MQ=255)
                             aaCAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATg                                                                 <  2:1701357/100‑1 (MQ=255)
                              agaataaaaaTTCTGACATCGTATTTTATGTACATAAAACACCAAAGAGCGTTGAGATATATATATATAAAAAAAAATATTATTTATAAAATAATTAATg                                                                 <  2:42110/93‑1 (MQ=255)
                              aCAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATg                                                                 <  2:1067584/100‑1 (MQ=255)
                              aCAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAAAATTCTAAAAAAAAATATTATTTATAAAATAATTAATg                                                                 <  2:1266420/100‑1 (MQ=255)
                                aaTTCAAATACTGACATCGTATTTTATGGACATAAAACACAAAAGTGCGTTAAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGAc                                                               <  2:466654/100‑1 (MQ=255)
                                 aTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGAcc                                                              >  2:1623111/1‑100 (MQ=255)
                                 aTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAAAATTATTTATAAaataatatttttt                                                               >  2:1448521/1‑92 (MQ=255)
                                  ttCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCa                                                             <  1:1753034/100‑1 (MQ=255)
                                    cAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCaaa                                                           >  2:1501049/1‑100 (MQ=255)
                                      aaTTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCaaaaa                                                         <  1:674507/100‑1 (MQ=255)
                                       aTTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCaaaaaa                                                        >  1:646454/1‑100 (MQ=255)
                                        ttCTAACATCGAATTTAATGGAAATAAAAAACAAAAGAGCGTTGAGAAATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAAt                                                       <  1:882218/100‑1 (MQ=255)
                                         tCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAAtt                                                      <  1:1961096/100‑1 (MQ=255)
                                          cTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATta                                                     <  1:1611818/100‑1 (MQ=255)
                                          cTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATta                                                     >  2:618246/1‑100 (MQ=255)
                                           tGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAACATtag                                                    >  1:1704806/1‑100 (MQ=255)
                                           tGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                    >  1:886753/1‑100 (MQ=255)
                                           tGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAACATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                    <  2:1089453/100‑1 (MQ=255)
                                            gACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagt                                                   <  1:1398050/100‑1 (MQ=255)
                                            gACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagt                                                   >  1:309219/1‑100 (MQ=255)
                                              cATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagtag                                                 >  1:72337/1‑100 (MQ=255)
                                               aTCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTaga                                                >  2:1270408/1‑100 (MQ=255)
                                               aTCGTATTTTATGGACATAAAAAACAAAAGAGCTTTGAGATATAACTATATGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagtag                                                 <  2:801808/99‑1 (MQ=255)
                                                  gTATTTGAGGGACATAAAATAACAAAGTGAATTGAGATATATATTTCTGAAAAAAAATATTATTTATAAAATAATTACAAAAAAAAAAATTAGTAGAggg                                             <  2:1186817/100‑1 (MQ=255)
                                                        tATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGt                                       >  1:1154734/1‑100 (MQ=255)
                                                         aTGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTc                                      <  2:1071156/100‑1 (MQ=255)
                                                          tGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCa                                     >  1:1120130/1‑100 (MQ=255)
                                                            gACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt                                   <  2:646452/100‑1 (MQ=255)
                                                            gACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAt                                    >  2:1783903/1‑99 (MQ=255)
                                                             gatctaaaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAttt                                  <  1:1698188/95‑1 (MQ=255)
                                                             aCATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAttt                                  <  2:309218/100‑1 (MQ=255)
                                                               aTAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTta                                <  2:1161434/100‑1 (MQ=255)
                                                               aTAAAACACCAAAGAGCGTTGAGATATATATTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTta                                <  2:105765/100‑1 (MQ=255)
                                                                tAAAACAGCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATCACCAAAAAATTAGTAGAGGGAATGGTCATTTtat                               >  2:1076581/1‑100 (MQ=255)
                                                                tAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACTAAAAAATTAGTAGTGGGAATGGTCATTTtat                               >  2:809636/1‑100 (MQ=255)
                                                                tAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtat                               >  1:1256424/1‑100 (MQ=255)
                                                                tAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtat                               >  2:21970/1‑100 (MQ=255)
                                                                 aaaaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATGTTATTTATAAAATCATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                              >  2:1433710/1‑100 (MQ=255)
                                                                  aaaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtttagatc                             >  2:1698175/1‑95 (MQ=255)
                                                                  aaaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATaa                             >  2:329756/1‑100 (MQ=255)
                                                                   aaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                            <  2:1154728/100‑1 (MQ=255)
                                                                   aaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                            <  1:1623124/100‑1 (MQ=255)
                                                                    acacCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGt                           <  1:1523378/100‑1 (MQ=255)
                                                                    acacCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGt                           >  1:12094/1‑100 (MQ=255)
                                                                    acacCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGt                           <  1:1980974/100‑1 (MQ=255)
                                                                    acacCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGt                           <  1:241562/100‑1 (MQ=255)
                                                                      acCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGAAAAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGta                          <  2:872687/99‑1 (MQ=255)
                                                                          aaGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATaaataatatt                     >  1:1662159/1‑92 (MQ=255)
                                                                          aaGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAAt                      >  2:492275/1‑99 (MQ=255)
                                                                           agagCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGTGAATGATCATTTTATAAGTATAATGGt                   >  2:1623973/1‑100 (MQ=255)
                                                                             agCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGtt                  >  2:1041542/1‑100 (MQ=255)
                                                                               cGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaa                >  1:546904/1‑100 (MQ=255)
                                                                               cGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaa                >  1:1100882/1‑100 (MQ=255)
                                                                                gTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaat               <  1:47504/100‑1 (MQ=255)
                                                                                gTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaat               <  1:329757/100‑1 (MQ=255)
                                                                                 ttGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaata              >  1:906683/1‑100 (MQ=255)
                                                                                  tGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaataa             >  1:35607/1‑100 (MQ=255)
                                                                                     gaTATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGAAATGGTCATTTTATAAGTATAATGGTTAATAATTa          >  2:1953692/1‑100 (MQ=255)
                                                                                     gaTATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTa          >  2:938035/1‑99 (MQ=255)
                                                                                      atatatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTa          <  2:1644470/99‑1 (MQ=255)
                                                                                       tatatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACa        >  2:1121234/1‑100 (MQ=255)
                                                                                        atatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAg       <  1:1349932/100‑1 (MQ=255)
                                                                                          atCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGa     >  2:267839/1‑100 (MQ=255)
                                                                                           tCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGAc    >  1:14976/1‑100 (MQ=255)
                                                                                            cTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACg   >  1:1071302/1‑100 (MQ=255)
                                                                                             tttCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGATGc  >  2:1199537/1‑100 (MQ=255)
                                                                                                   |                                                                                              
GTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGC  >  NC_000913/1529971‑1530162

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