Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,530,076:1 (A)7→8 100% 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.A96.3% 369.8 / 6.3 109coding (148/483 nt)ydcDputative immunity protein for RhsE
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base A (54/51);  total (55/54)
Fisher's exact test for biased strand distribution p-value = 3.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 122 total aligned reads displayed.

TGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGAC  >  NC_000913/1529970‑1530160
                                                                                                     |                                                                                             
tGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTg                                                                                                >  1:111367/1‑100 (MQ=255)
tGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTg                                                                                                <  1:288101/100‑1 (MQ=255)
 gTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑a                                                                                             <  1:3210064/100‑1 (MQ=255)
         gttgttTACCTTGGGTTTGATAACAATTCACATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aaaaaaaat                                                                                    >  1:1731627/1‑98 (MQ=255)
         gttgttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aaaaaaaat                                                                                    <  2:1043849/100‑3 (MQ=255)
         gttgttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aaaaaaaat                                                                                    <  2:2791320/100‑3 (MQ=255)
           tgttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aaaaaaaatat                                                                                  >  2:3141953/1‑96 (MQ=255)
            gttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aaaaaaaatatt                                                                                 >  1:10431/1‑95 (MQ=255)
            gttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aaaaaaaatatt                                                                                 >  2:1359595/1‑95 (MQ=255)
            gttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aaaaaaaatat                                                                                  >  1:977924/1‑95 (MQ=255)
            gttTACCTTGGGTGTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aaaaaaaatatt                                                                                 >  1:2091088/1‑95 (MQ=255)
             tttACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAAtatta                                                                                 >  1:1793624/1‑100 (MQ=255)
               tACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATattatt                                                                               <  2:1143261/100‑1 (MQ=255)
               tACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATattatt                                                                               <  2:1179680/100‑1 (MQ=255)
                 ccTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATAttattta                                                                             <  2:1209456/100‑1 (MQ=255)
                   ttGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTtata                                                                           >  2:2177899/1‑100 (MQ=255)
                      ggTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATaaaa                                                                        >  2:1506320/1‑100 (MQ=255)
                       gTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAaat                                                                       <  1:2580318/100‑1 (MQ=255)
                       gTTTGATAACAATTCAAATTCTGACAGCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAaat                                                                       >  2:2060618/1‑100 (MQ=255)
                        tttGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAaata                                                                      >  2:3015238/1‑100 (MQ=255)
                           gATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaatt                                                                   >  1:2540984/1‑100 (MQ=255)
                           gATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaatt                                                                   >  2:2669487/1‑100 (MQ=255)
                            aTAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaatta                                                                  >  1:583669/1‑100 (MQ=255)
                             tAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaattaa                                                                 >  1:2080470/1‑100 (MQ=255)
                             tAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaattaa                                                                 >  1:1294474/1‑100 (MQ=255)
                              aaCAATTCAAATTCTGCCATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaattaat                                                                <  2:1524903/100‑1 (MQ=255)
                               aCAATTCAAATTCTGACATCGTATTTTAGGGACATAAAACACCAAAGGGCGTTGAGATATGTCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATg                                                               >  2:2700557/1‑100 (MQ=255)
                                 aaTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGAc                                                             >  2:2440206/1‑100 (MQ=255)
                                 aaTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGAc                                                             <  2:1894129/100‑1 (MQ=255)
                                  aTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAACGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaattaatta                                                              >  1:1870021/1‑96 (MQ=255)
                                  aTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGAcc                                                            <  1:1921523/100‑1 (MQ=255)
                                  aTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGAcc                                                            >  2:1670073/1‑100 (MQ=255)
                                  aTTCAAATTCTGACATCGTATTTTATGGACATAAAACAACAAAGAGCGTTGAGATATATCTATATGAAAAAAAATATTATTTATAAAATAATTAATGAcc                                                            <  1:2854997/100‑1 (MQ=255)
                                     cAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCaaa                                                         >  2:1431231/1‑100 (MQ=255)
                                      aaaTTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATTTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCaaaa                                                        >  1:1967126/1‑100 (MQ=255)
                                      aaaTTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCaaaa                                                        >  1:732912/1‑100 (MQ=255)
                                      aaaTTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCaaaa                                                        <  2:1274615/100‑1 (MQ=255)
                                         ttCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAAt                                                     <  2:2376223/100‑1 (MQ=255)
                                           cTGACATCGTATTTTATGGACATAAAACACCATATAGGGTTGAGATATATTTTTATGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATta                                                   <  1:1367639/100‑1 (MQ=255)
                                            tGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                  <  2:2954100/100‑1 (MQ=255)
                                            tGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                  >  2:1819973/1‑100 (MQ=255)
                                             gACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagt                                                 >  2:2955495/1‑100 (MQ=255)
                                             gACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagt                                                 <  2:289038/100‑1 (MQ=255)
                                             gACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                  >  1:1609284/1‑99 (MQ=255)
                                             gACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagta                                                <  2:810754/100‑1 (MQ=255)
                                              aCATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagta                                                >  2:467444/1‑100 (MQ=255)
                                                aTCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTaga                                              >  1:2909402/1‑100 (MQ=255)
                                                 tCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTagag                                             <  2:226838/100‑1 (MQ=255)
                                                  cGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTagag                                             <  2:2940697/100‑1 (MQ=255)
                                                  cGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTagag                                             >  1:1627593/1‑99 (MQ=255)
                                                   gTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATTAAATAATTAATGACCAAAAAATTAGTAGAggg                                           <  2:3213211/100‑1 (MQ=255)
                                                     aTTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGaa                                         >  1:1233533/1‑100 (MQ=255)
                                                     aTTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGaa                                         >  2:1195663/1‑100 (MQ=255)
                                                      ttttATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAAt                                        <  2:3134720/100‑1 (MQ=255)
                                                        ttATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATgg                                      <  2:1868794/100‑1 (MQ=255)
                                                         tATGGACATAAAACACCAAAGAGCGTTGAAATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGt                                     <  1:528767/100‑1 (MQ=255)
                                                             gACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGtaa                                   >  1:1427114/1‑96 (MQ=255)
                                                              aCATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGTATGGTCAttt                                >  1:1844919/1‑100 (MQ=255)
                                                                 tAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtat                             >  2:2445189/1‑100 (MQ=255)
                                                                  aaaaaaCCAAAGAGCGTTGAGATATATATTAAAAA‑‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                            <  1:2873398/100‑1 (MQ=255)
                                                                  aaaaaaCCAAAGAGCGTTAAAATATATTTTTATGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                            <  1:87684/100‑1 (MQ=255)
                                                                  aaaaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                            <  2:561082/100‑1 (MQ=255)
                                                                  aaaaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                            <  2:3185285/100‑1 (MQ=255)
                                                                   aaaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATaa                           <  1:2418614/100‑1 (MQ=255)
                                                                   aaaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATaa                           >  2:2741179/1‑100 (MQ=255)
                                                                   aaaCACCAAAGAGCGTTGAGATATATCTTTCTG‑‑AAGAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATaa                           <  2:150053/100‑1 (MQ=255)
                                                                    aaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                          <  2:2129100/100‑1 (MQ=255)
                                                                    aaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                          <  2:1313986/100‑1 (MQ=255)
                                                                      cacCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGt                         <  1:618218/100‑1 (MQ=255)
                                                                      cacCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGta                        <  2:1397782/100‑1 (MQ=255)
                                                                       acCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGtat                       >  2:2767795/1‑100 (MQ=255)
                                                                          aaaGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAAt                    <  2:3353715/100‑1 (MQ=255)
                                                                          aaaGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAAt                    <  2:2010087/100‑1 (MQ=255)
                                                                           aaGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATg                   >  2:3157273/1‑100 (MQ=255)
                                                                            agagCGTTGAGATCTATCTTTCAGAAAAAAAATATTATTTATAAAATAATTAATGATCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATgg                  <  1:2404955/100‑1 (MQ=255)
                                                                            agagCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATgg                  <  2:2541004/100‑1 (MQ=255)
                                                                             gagCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATgg                  <  1:2405345/99‑1 (MQ=255)
                                                                               gCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTta               <  2:2091101/100‑1 (MQ=255)
                                                                               gCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTta               >  1:3125783/1‑100 (MQ=255)
                                                                               gCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTta               >  2:686214/1‑100 (MQ=255)
                                                                                cGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaa              >  2:1075407/1‑100 (MQ=255)
                                                                                cGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaa              >  2:2474350/1‑100 (MQ=255)
                                                                                cGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaa              >  2:699277/1‑100 (MQ=255)
                                                                                 gTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACAAAAAAATTAGAAGAGGGAAATATAATTTTATAAGTATAATGGTtaat             >  1:1209000/1‑100 (MQ=255)
                                                                                  ttGGGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaata            <  2:1695393/100‑1 (MQ=255)
                                                                                  ttGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaata            >  1:1084471/1‑100 (MQ=255)
                                                                                    gagaTATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTATAATGAAAGGTCATTTTATAAGTATAAt                    >  1:2835967/1‑90 (MQ=255)
                                                                                    gagaTATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGAACAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTaataat          <  2:2577648/100‑1 (MQ=255)
                                                                                    gagaTATATATTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTaataat          <  1:3138821/100‑1 (MQ=255)
                                                                                      gaTATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATTATCAAAAAATTATTAGAGGTAATGGTAATTTTATAAGTATAATGGTTAATAATTa        >  1:1654569/1‑100 (MQ=255)
                                                                                       atatatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTAc       <  1:1136469/100‑1 (MQ=255)
                                                                                       atatatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTAc       <  2:2125205/100‑1 (MQ=255)
                                                                                        tatatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACa      >  2:2835717/1‑100 (MQ=255)
                                                                                        tatatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACa      >  2:3141133/1‑100 (MQ=255)
                                                                                        tatatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACa      <  2:1014199/100‑1 (MQ=255)
                                                                                        tatatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAAAAACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACa      <  1:2159696/100‑1 (MQ=255)
                                                                                         atatCTTTATAAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACa      <  1:1115374/99‑1 (MQ=255)
                                                                                            tCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGAc  >  2:708910/1‑100 (MQ=255)
                                                                                             cTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTTTAGGGAATGGTTATTTTATAAGTATAATGGTTAATAAtt         >  1:3347602/1‑92 (MQ=255)
                                                                                             cTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAAGGAATGGACATATTATAAGTATAATGGTTAATAATTACAgg    >  1:2002491/1‑97 (MQ=255)
                                                                                                     |                                                                                             
TGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGAC  >  NC_000913/1529970‑1530160

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