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.A97.6% 422.8 / 2.7 125coding (148/483 nt)ydcDputative immunity protein for RhsE
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base A (65/57);  total (67/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 145 total aligned reads displayed.

TGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCA  >  NC_000913/1529970‑1530163
                                                                                                    |                                                                                                
tGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTg                                                                                                   <  1:293219/100‑1 (MQ=255)
tGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTg                                                                                                   <  2:1879953/100‑1 (MQ=255)
 gTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑a                                                                                                <  2:628982/100‑1 (MQ=255)
  tttATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aa                                                                                               <  2:675142/100‑1 (MQ=255)
    tATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aaaa                                                                                             >  1:81859/1‑100 (MQ=255)
     aTTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aaaaa                                                                                            >  2:731156/1‑100 (MQ=255)
      tttGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aaaaaa                                                                                           >  2:2629581/1‑100 (MQ=255)
      tttGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aaaaaa                                                                                           <  2:1137105/100‑1 (MQ=255)
        tgttgtTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aaaaaaaa                                                                                         >  2:2202187/1‑99 (MQ=255)
        tgttgtTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aaaaaaaa                                                                                         <  1:1183801/100‑2 (MQ=255)
         tttgttTACCTGTAATTTAATAACAATACAAATTCTGAAATCCAATTTTATAGACATAAAACACCAAAGAGCGGCCAGATAGATCATTCTG‑‑aaaaaa                                                                                           <  1:1425798/96‑1 (MQ=255)
         gttgttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aaaaaaaat                                                                                        >  2:1930799/1‑98 (MQ=255)
          ttgttTATATTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACAACAAAGAGAGTTGAGATATATCTTTATG‑‑aaaaaaaat                                                                                        <  1:1864911/99‑3 (MQ=255)
            gttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aaaaaaaatatt                                                                                     >  2:1421645/1‑95 (MQ=255)
               tACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGGGCGTTGAGATATATCTTTCTGAAAAAAAATattatt                                                                                  <  2:672915/100‑1 (MQ=255)
                aCCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATattatt                                                                                  <  1:3388182/99‑1 (MQ=255)
                aCCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATAttattt                                                                                 <  2:2358778/100‑1 (MQ=255)
                aCCTTAGGTTTGATAACAATTCAAATTCTGACATATTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATAttattt                                                                                 <  1:3509177/100‑1 (MQ=255)
                 ccTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATAttattta                                                                                <  2:640256/100‑1 (MQ=255)
                 ccTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATAttattta                                                                                >  2:3180708/1‑100 (MQ=255)
                  cTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATAttatttat                                                                               <  2:1873253/100‑1 (MQ=255)
                   ttGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTtata                                                                              >  1:188264/1‑100 (MQ=255)
                     gggTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATaaa                                                                            >  2:1585493/1‑100 (MQ=255)
                     gggTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATaaa                                                                            <  2:3082840/100‑1 (MQ=255)
                      ggTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATaaaa                                                                           >  2:1544867/1‑100 (MQ=255)
                       gTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAaat                                                                          <  2:4140990/100‑1 (MQ=255)
                       gTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAaat                                                                          >  2:4246025/1‑100 (MQ=255)
                         ttGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAaataa                                                                        >  1:1698251/1‑100 (MQ=255)
                           gATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaatt                                                                      >  2:2709988/1‑100 (MQ=255)
                           gATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaatt                                                                      <  2:3359636/100‑1 (MQ=255)
                           gATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaatt                                                                      <  2:365115/100‑1 (MQ=255)
                            aTAACAATTCAAATTCTGACATCGTATTTTATAAAAATAAAACACCAAAGAGCGTTGAGATATATATTTAAGAAAAAAAATATTATTTATAAAAtaatt                                                                      <  1:668843/99‑1 (MQ=255)
                            aTAACAATTCAAATTCTGACATCGTATATTATGAACATAAAACACCAAAGAGCTTTAAGATATATCAAAAAAAAAAAAAATATTATTTATAAAAtaatta                                                                     <  1:1084778/100‑1 (MQ=255)
                            aTAAAAATTCAAATTCTTAAATCCTATTTTATAGACATAAAACAACAAAGAGTTTTGAGATATAAATTAATAAAAAAAAATATTATTTATAAaataat                                                                       <  1:527984/98‑1 (MQ=255)
                             tagcAATTCAAATTCTTAAATCGTATTTTATATACATAAAACACCAAAGAGCTTTTAGATAGAACTTTCTGAAAAAAAATATTATTTATAAAAtaattaa                                                                    <  1:2063342/97‑1 (MQ=255)
                             tAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaattaa                                                                    >  2:138163/1‑100 (MQ=255)
                              aaCAATTCAAATTCTGACATCTTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaattaat                                                                   <  1:273195/100‑1 (MQ=255)
                              aaCAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaattaa                                                                    >  2:1285449/1‑99 (MQ=255)
                               aCAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATg                                                                  >  2:1434645/1‑100 (MQ=255)
                                cAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGa                                                                 <  2:906537/100‑1 (MQ=255)
                                 aaTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGAc                                                                <  2:2271874/100‑1 (MQ=255)
                                  aTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGAcc                                                               <  2:1449943/100‑1 (MQ=255)
                                   ttCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCa                                                              >  2:925704/1‑100 (MQ=255)
                                       aaTTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCaaaaa                                                          >  1:3923995/1‑100 (MQ=255)
                                         ttCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAAt                                                        <  2:1650853/100‑1 (MQ=255)
                                         ttCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAAt                                                        <  2:1301228/100‑1 (MQ=255)
                                         ttCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAAt                                                        <  2:3850269/100‑1 (MQ=255)
                                           cTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATta                                                      >  2:1982108/1‑100 (MQ=255)
                                            tGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                     >  2:4221826/1‑100 (MQ=255)
                                            tGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                     >  2:3288671/1‑100 (MQ=255)
                                            tGACATCGTATTTTATGAACATAAAACACCAAAGAGCGTTGAGATATAAAAAAAAAAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                     <  1:1143992/100‑1 (MQ=255)
                                              aCATAGAATTTAATGGATATAAAAAACCAAAGAGAGTTTATATATAAAAATAAAAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATta                                                      <  1:229142/97‑1 (MQ=255)
                                                ttggtgttttATGTCCATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTaga                                                 >  2:2409187/7‑100 (MQ=255)
                                                aTCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTaga                                                 >  1:1978546/1‑100 (MQ=255)
                                                aTCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTaga                                                 >  2:2029139/1‑100 (MQ=255)
                                                aTCGTATTTTATGGACATAAAACACAAAAGAGCGTTGAGATATAAAAAAAAAAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTaga                                                 <  1:3636600/100‑1 (MQ=255)
                                                 tCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTagag                                                >  1:2215573/1‑100 (MQ=255)
                                                 tCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTagag                                                <  2:528851/100‑1 (MQ=255)
                                                  cGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAgg                                               >  1:954786/1‑100 (MQ=255)
                                                   gTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATATATAAAATAATTAATGACCAAAAAATTAGTAGAggg                                              >  2:2999726/1‑100 (MQ=255)
                                                     aTTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGaa                                            <  1:4246003/100‑1 (MQ=255)
                                                      ttttATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAACAATTAGTAGAGGGAAt                                           <  2:2326428/100‑1 (MQ=255)
                                                      ttttATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAAt                                           >  2:1026792/1‑100 (MQ=255)
                                                      ttttATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAAt                                           >  2:3129507/1‑100 (MQ=255)
                                                         tATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGt                                        >  1:4147187/1‑100 (MQ=255)
                                                            ggACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTc                                       >  1:1186701/1‑98 (MQ=255)
                                                            ggACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCa                                      >  1:2942548/1‑99 (MQ=255)
                                                            ggACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAt                                     >  1:799226/1‑100 (MQ=255)
                                                            ggACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAt                                     >  2:1733335/1‑100 (MQ=255)
                                                             gACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGAAGAGGGAATgg                                         >  1:1718655/1‑95 (MQ=255)
                                                              aCATAAAACACCAAAGAGCGTTAAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt                                    <  1:1928733/99‑1 (MQ=255)
                                                               cATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtttt                                  <  2:2619849/100‑1 (MQ=255)
                                                                 tAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGCCCAAAAAAATAGAAGAGGGAAAGGTCATTTtat                                >  2:3485377/1‑100 (MQ=255)
                                                                 tAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtat                                <  2:72038/100‑1 (MQ=255)
                                                                  aaaaCATTAAAGAGCGTAAAGATATATCTTTATGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                               <  1:1774651/100‑1 (MQ=255)
                                                                   aaaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATaa                              <  2:2655484/100‑1 (MQ=255)
                                                                    aaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                             >  1:1121928/1‑100 (MQ=255)
                                                                    aaCACCAAAGAGCGTTGAGATATAAATATATAAAAAAAAATATTATTTATAAAATAATTAATGAAAAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                             <  1:3366667/100‑1 (MQ=255)
                                                                       acCAAAGAGCGTTGAGATATAAAAAAAAAAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGta                           <  1:691207/99‑1 (MQ=255)
                                                                        ccAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGtata                         >  2:2217845/1‑100 (MQ=255)
                                                                           aaGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACAAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAAt                       >  1:4240585/1‑99 (MQ=255)
                                                                            agagCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATgg                     >  1:800167/1‑100 (MQ=255)
                                                                              agCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGtt                   <  1:1304632/100‑1 (MQ=255)
                                                                                cGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaa                 >  2:3365292/1‑100 (MQ=255)
                                                                                cGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaa                 <  1:4051694/100‑1 (MQ=255)
                                                                                cGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTta                  <  1:3628747/99‑1 (MQ=255)
                                                                                cGTTGAGATATATCTTTCCGAAAAAAAATCTTATTTATAAAATAATTAATGAACAAAAAATTAGTAGACGGAATGGTCATTTTATAAGTACAATGGTta                  >  1:3974926/1‑99 (MQ=255)
                                                                                 gTTGAGATATATCTTTCTG‑‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaata               >  2:1245972/1‑100 (MQ=255)
                                                                                  ttGAGATATATCTTTCTGAAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaat                <  2:1701969/100‑1 (MQ=255)
                                                                                   tGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaat                <  1:3634318/98‑1 (MQ=255)
                                                                                   tGAGATATATCTTTCTG‑‑AAAAAAACTATTATTCATAAAATAATTAATGACCAAAACATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtcataa              <  2:194793/100‑1 (MQ=255)
                                                                                     agaTATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaataan             <  1:2127804/99‑2 (MQ=255)
                                                                                     agaTATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAAtt            >  1:2786819/1‑100 (MQ=255)
                                                                                      gaTATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTa           >  1:3844162/1‑100 (MQ=255)
                                                                                         atatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATAAGTATTATAAATGATCATTTTATAAAAAAAATGGTTAATAATTa           >  1:3203880/1‑97 (MQ=255)
                                                                                          tatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAgg       >  2:1978285/1‑100 (MQ=255)
                                                                                            tCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGAc     >  1:71288/1‑100 (MQ=255)
                                                                                              tttCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGc   <  2:537709/100‑1 (MQ=255)
                                                                                               ttCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAACAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGc   <  2:2380279/99‑1 (MQ=255)
                                                                                               ttCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCa  >  2:3839736/1‑100 (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: