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.A92.2% 320.3 / 20.3 103coding (148/483 nt)ydcDputative immunity protein for RhsE
Reads supporting (aligned to +/- strand):  ref base . (3/5);  new base A (46/49);  total (49/54)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 117 total aligned reads displayed.

TGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCA  >  NC_000913/1529970‑1530163
                                                                                                       |                                                                                               
tGTTTATTTGTTGTTTACCTTGTGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACATCAAAGAGCGTTTAGATATATCTTTCTg                                                                                                  <  2:401389/100‑1 (MQ=255)
tGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTg                                                                                                  >  1:2333795/1‑100 (MQ=255)
tGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTg                                                                                                  <  2:1017355/100‑1 (MQ=255)
tGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTAAGATATATCTTTCTg                                                                                                  >  1:674441/1‑100 (MQ=255)
   ttATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aaa                                                                                             <  1:1215267/100‑1 (MQ=255)
   ttATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTATCTG‑‑aaa                                                                                             <  1:1586304/100‑1 (MQ=255)
      tttGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aaaaaa                                                                                          >  1:1387385/1‑100 (MQ=255)
       ttgttgTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aaaaaaa                                                                                         >  2:1014957/1‑100 (MQ=255)
         tttatttACCTTTGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aaaaaaaa                                                                                        <  1:2845986/95‑2 (MQ=255)
         gttgttTACCTTGGTTTTGATAATAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aaaaaaaat                                                                                       <  1:730440/100‑3 (MQ=255)
            gttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aaaaaaaatatt                                                                                    >  2:2836419/1‑95 (MQ=255)
            gttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aaaaaaaatatt                                                                                    >  1:1348287/1‑95 (MQ=255)
                aCCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATAttattt                                                                                 >  2:2582335/1‑100 (MQ=255)
                  cTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATAttatttat                                                                               <  2:2658408/100‑1 (MQ=255)
                  cTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATAttatttat                                                                               <  2:2658392/100‑1 (MQ=255)
                   ttGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGTGTTGAGATATATCTTTCTGAAAAAAAATATTATTtata                                                                              <  2:1035057/100‑1 (MQ=255)
                     ggttgtGAAAACTATTAGAATTCTGTCATGGTATTTTATGGACAGAAAACACAAAATAACTTTGAGATATATCTTTCTG‑‑AAAAAAACATTATTTATaa                                                                             <  2:2442134/93‑1 (MQ=255)
                       gTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAaat                                                                          >  1:1101481/1‑100 (MQ=255)
                           gATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaatt                                                                      >  1:2679370/1‑100 (MQ=255)
                            aTAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaatta                                                                     <  1:2863685/100‑1 (MQ=255)
                             tAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaattaa                                                                    >  2:2786399/1‑100 (MQ=255)
                             tAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACAACAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaattaa                                                                    >  2:1185436/1‑100 (MQ=255)
                              aaCAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaattaat                                                                   >  2:1401007/1‑100 (MQ=255)
                                cAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGa                                                                 <  1:1183843/100‑1 (MQ=255)
                                 aaTTCAAATTCTGACATTGTATTTTATGGAAATAAAATACAAAATATTGTTGAGATATATCTTTATTA‑‑AAAAAAATTTTATTTATAAAATAATTAATGAc                                                                <  2:1145273/100‑1 (MQ=255)
                                 aaTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGAc                                                                <  2:1238196/100‑1 (MQ=255)
                                 aaTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGAc                                                                <  1:937659/100‑1 (MQ=255)
                                  aTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAAATATTATTTATAAAATAATTAATGAc                                                                <  1:947318/100‑1 (MQ=255)
                                  aTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGAcc                                                               >  2:1318909/1‑100 (MQ=255)
                                  aTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGAcc                                                               >  1:2867311/1‑100 (MQ=255)
                                  aTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAAATATATCTTTATGAAAAAAAATATTATTTATAAAATAATTAATGa                                                                 >  1:135245/1‑98 (MQ=255)
                                    tCAAATTTGGATATCGTATTTTATGTACATAAAAAACGAAAGAGCTTTGATATATATCTTTCTGAAAAAAAATATTATTTTATAaataatattttttttc                                                             <  2:1323703/100‑11 (MQ=255)
                                      aaaTTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCaaaa                                                           <  2:764633/100‑1 (MQ=255)
                                         ttCTGACATCTTATTTTATGGACATAAAATACCAAATAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAAt                                                        <  1:391888/100‑1 (MQ=255)
                                           cTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGATATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATta                                                      <  2:1100647/100‑1 (MQ=255)
                                             gACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagt                                                    <  1:2233407/100‑1 (MQ=255)
                                             gACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagt                                                    <  1:2263790/100‑1 (MQ=255)
                                              aCATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagta                                                   <  1:2494593/100‑1 (MQ=255)
                                               gatgatATTTTATTGACATAAAAGACCAAAGAGCATTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAAAGAAAAAAAAATtagt                                                    <  2:2445465/93‑1 (MQ=255)
                                                aTCGTATTTTATGGACATAACACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTaga                                                 >  1:2864144/1‑100 (MQ=255)
                                                 tggtgttttATGTCCATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAgg                                               >  2:1870065/6‑100 (MQ=255)
                                                 tCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTagag                                                >  1:1650283/1‑100 (MQ=255)
                                                  cGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTagag                                                >  2:2056722/1‑99 (MQ=255)
                                                   gTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTaga                                                 >  1:739230/1‑97 (MQ=255)
                                                   gTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAggg                                              <  2:477847/100‑1 (MQ=255)
                                                    tATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGa                                             <  2:2631294/100‑1 (MQ=255)
                                                     aTTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGaa                                            >  1:2990448/1‑100 (MQ=255)
                                                     aTTTTATAGACATAAAAAACTAAAGAGCGATGAGATATATATTTCTGAAAAAAAATATTATTAATAAAATAAGTAATTACAAAAACAGAAGTAGAGGTaa                                            >  2:724741/1‑100 (MQ=255)
                                                         tATGGACATAAAAAACCAAAGGGCGGTGAGATATATCTGTTGGA‑‑AAAAAAATATTATTTATTAAATAATTAATAAAAAAAAAATTAGTAGAGGGAATgg                                         <  2:2644323/99‑1 (MQ=255)
                                                          aTGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGt                                        <  1:901316/100‑1 (MQ=255)
                                                           tGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCa                                      >  1:1035542/1‑100 (MQ=255)
                                                             gACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt                                    >  2:1541467/1‑100 (MQ=255)
                                                             gACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt                                    >  1:1564408/1‑100 (MQ=255)
                                                              aCATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAttt                                   <  1:1647724/100‑1 (MQ=255)
                                                              aCATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAttt                                   <  2:1680326/100‑1 (MQ=255)
                                                              aCATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAttt                                   >  2:203398/1‑100 (MQ=255)
                                                               cATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATAACCAAAAAATTAGTAGAGGAAATGGTCAtttt                                  >  1:159137/1‑100 (MQ=255)
                                                                 tAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtat                                >  2:2206547/1‑100 (MQ=255)
                                                                 tAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtat                                >  2:146440/1‑100 (MQ=255)
                                                                  aaaaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                               <  2:843839/100‑1 (MQ=255)
                                                                  aaaaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                               <  2:2207938/100‑1 (MQ=255)
                                                                  aaaaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATAACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                               <  2:1641073/100‑1 (MQ=255)
                                                                   aaaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATaa                              >  2:2474961/1‑100 (MQ=255)
                                                                   aaaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATaa                              <  1:3028935/100‑1 (MQ=255)
                                                                    aaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                             <  1:1543467/100‑1 (MQ=255)
                                                                     acacCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                             <  1:2574070/99‑1 (MQ=255)
                                                                      cacCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTATAGGGAATGGTCATTTTATAAGta                           >  2:2697000/1‑100 (MQ=255)
                                                                      cacCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGta                           >  2:3013740/1‑100 (MQ=255)
                                                                      cacCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGta                           >  2:861789/1‑100 (MQ=255)
                                                                      cacCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGta                           >  1:513936/1‑100 (MQ=255)
                                                                      cacCAAAGAGCGTTGAGATATATCTTTCTG‑‑AAAAAAATATTATTTATAAAATAATTAATAACCAAAAAATTAGTAGAGGGAATGGTCATTTGATAAgaa                           >  2:2762847/1‑97 (MQ=255)
                                                                       acCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGtat                          >  2:1685259/1‑100 (MQ=255)
                                                                        ccAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGtata                         >  1:266767/1‑100 (MQ=255)
                                                                          aaaGAGCGTTGAGATATATCTTTATGAAAAAAAATATTATTTATAAAATAATTAATGAAAAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGtata                         <  2:3037229/98‑1 (MQ=255)
                                                                            agagCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTATAGGGAATGGTCATTTTATAAGTATAATgg                     >  2:2080150/1‑100 (MQ=255)
                                                                            agagCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATgg                     <  1:210777/100‑1 (MQ=255)
                                                                              agCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGtt                   <  1:2206570/100‑1 (MQ=255)
                                                                                cGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGtt                   <  2:463296/98‑1 (MQ=255)
                                                                                cGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaa                 >  2:2629136/1‑100 (MQ=255)
                                                                                 gTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGATGGAATGATCATTTTATAAGTATAATGGTta                  >  2:2107633/1‑98 (MQ=255)
                                                                                 gTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTta                  <  2:2969001/98‑1 (MQ=255)
                                                                                 cttGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaa                 <  1:1558783/98‑1 (MQ=255)
                                                                                  ttGAGATATATATTTCTAA‑‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaat                <  2:695723/99‑1 (MQ=255)
                                                                                     agaTATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAAtt            >  2:1111131/1‑100 (MQ=255)
                                                                                      gaTATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTTATGTCCTAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTa           <  1:2474990/100‑1 (MQ=255)
                                                                                      gaTATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACAAAAAAAGTAGTAGAGGGAATTGTCATATTATAAGTAGAATGGTtaataa              >  2:2606507/1‑97 (MQ=255)
                                                                                       atatatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTa           <  1:1234115/99‑1 (MQ=255)
                                                                                       atatatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTAc          <  2:1778917/100‑1 (MQ=255)
                                                                                       atatatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTAc          <  1:1155589/100‑1 (MQ=255)
                                                                                         atatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAg        <  2:2148924/100‑1 (MQ=255)
                                                                                         atatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAg        <  1:2066701/100‑1 (MQ=255)
                                                                                          tatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAAGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAg        <  1:281972/99‑1 (MQ=255)
                                                                                             cTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACg    >  2:2819182/1‑100 (MQ=255)
                                                                                             cTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACg    >  2:1546873/1‑100 (MQ=255)
                                                                                              tttCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGc   >  2:2412299/1‑100 (MQ=255)
                                                                                              tttCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGc   >  1:1472790/1‑100 (MQ=255)
                                                                                              tttCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGc   <  1:167036/100‑1 (MQ=255)
                                                                                              tttCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTAATAAGTATAATGGTTAATAATTACAGGa      >  2:2456037/1‑97 (MQ=255)
                                                                                                ttCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCa  >  2:2717889/1‑100 (MQ=255)
                                                                                                ttCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCa  <  1:550191/100‑1 (MQ=255)
                                                                                                       |                                                                                               
TGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCA  >  NC_000913/1529970‑1530163

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