Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 285,204 C→T 16.2% intergenic (+36/‑191) yagF → / → yagG CP4‑6 prophage; dehydratase family protein/CP4‑6 prophage; putative sugar transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913285,2040CT16.2% 299.0 / 42.0 135intergenic (+36/‑191)yagF/yagGCP4‑6 prophage; dehydratase family protein/CP4‑6 prophage; putative sugar transporter
Reads supporting (aligned to +/- strand):  ref base C (43/69);  new base T (8/14);  total (51/84)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.07e-03

Only 100 of 147 total aligned reads displayed.

TATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATA  >  NC_000913/285105‑285303
                                                                                                   |                                                                                                   
tATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTc                                                                                                     >  1:1500528/1‑100 (MQ=255)
tATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAAACGTCTGCGGAATATTTc                                                                                                     >  2:3162156/1‑100 (MQ=255)
 aTTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTcc                                                                                                    <  2:2660661/100‑1 (MQ=255)
 aTTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTcc                                                                                                    >  1:1699978/1‑100 (MQ=255)
 aTTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTcc                                                                                                    <  2:2301603/100‑1 (MQ=255)
 aTTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAAATAAAACCGTCTGCGGAATATTTcc                                                                                                    >  1:2255715/1‑100 (MQ=255)
  tttATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTccc                                                                                                   <  2:2493229/100‑1 (MQ=255)
   ttATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCg                                                                                                  <  1:384234/100‑1 (MQ=255)
    tATGACACCGATAAAATTATCGAGGTAATTAACGCCGTTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGc                                                                                                 <  1:3162186/100‑1 (MQ=255)
    tATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGc                                                                                                 <  2:359876/100‑1 (MQ=255)
    tATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGc                                                                                                 <  2:1194447/100‑1 (MQ=255)
    tATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGc                                                                                                 <  1:3082045/100‑1 (MQ=255)
    tATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGc                                                                                                 >  1:2804899/1‑100 (MQ=255)
     aTGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCa                                                                                                <  2:259874/100‑1 (MQ=255)
       gACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGa                                                                                              <  2:425289/100‑1 (MQ=255)
        acacCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAATCGCTCGAAATTTAATTATTTTAATAGATAAAACCGTCTGCGGGAAATATTCCGCAGAc                                                                                             <  2:211389/100‑1 (MQ=38)
        acacCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGAc                                                                                             <  1:2482083/100‑1 (MQ=38)
        acacCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGAc                                                                                             <  1:873738/100‑1 (MQ=255)
         cacCGATAAAATTATCGATGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAGCCGTCTGCGGGAAATATTCCGCAGACg                                                                                            <  1:413438/100‑1 (MQ=37)
         cacCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACg                                                                                            >  1:147831/1‑100 (MQ=255)
         cacCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAACGCTCGGAATTTAATTATTTTAAGAGATAAAACCCTCTGAGGAATATTTCTCGCAGACg                                                                                            >  1:2080345/1‑100 (MQ=255)
          acCGATAAAATTATCGAGGTAATTAACGTCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGACgg                                                                                           <  1:38335/100‑1 (MQ=38)
          acCGATAAAATTATCGAGGTAATTAACGACGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGACg                                                                                            <  2:899520/99‑1 (MQ=38)
            cgtcaataTTATGTAGGTAATTTACGCTGGTAGAGAGGGTCTGTGTATTTAATTTTTTTAATAGATAAAACCGTCTGCGGGAAATATTCCGCAGACg                                                                                            <  2:3617555/90‑1 (MQ=2)
            cGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGACggt                                                                                          <  2:69511/99‑2 (MQ=38)
              aTAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGTGGAATATTTCCCGCAGACGGCttt                                                                                       >  1:3641495/1‑100 (MQ=255)
               tAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGACGGtttta                                                                                      <  2:1096635/100‑2 (MQ=37)
               tAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATATTAAGAGATAAAACCGTCTACGGAATATTTCCCGCAGACGGCtt                                                                                        >  1:267991/1‑98 (MQ=255)
                 aaaTTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTttgtt                                                                                    >  1:2893724/1‑100 (MQ=255)
                 aaaTTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTttgtt                                                                                    >  1:1925596/1‑100 (MQ=255)
                 aaaTTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCCttgtt                                                                                    >  1:1546318/1‑100 (MQ=255)
                   aaaatcaatgTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGCCGGCTTtgttgt                                                                                  <  2:718595/91‑1 (MQ=255)
                   aTTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTttgtt                                                                                    <  2:2255695/98‑1 (MQ=255)
                    ttATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAGCCGTCTGCGGGAAATATTCCGCAGACGGttttatctct                                                                                 <  1:3302035/100‑7 (MQ=18)
                       tCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAGCCGTCTGCGGGAAATATTCCGCAGACGGttttatctcttaa                                                                              <  1:2940604/100‑10 (MQ=18)
                         gAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGa                                                                            >  1:3005075/1‑100 (MQ=255)
                         gAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGa                                                                            <  2:3247404/100‑1 (MQ=255)
                          aGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGaa                                                                           <  1:2672468/100‑1 (MQ=255)
                          aGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGaa                                                                           <  1:3494492/100‑1 (MQ=255)
                            gTAATTAATGCCGGTAAAAAAGCGTTCGGTATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAAt                                                                         <  1:317625/100‑1 (MQ=255)
                            gTAATTAACGCCGGTAAAAAAGCGCTTGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAAt                                                                         <  1:517511/100‑1 (MQ=255)
                             tAATTAACGCCGGTAAAAGAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAAtt                                                                        >  1:1341945/1‑100 (MQ=255)
                             tAATTAACGCCGGTAAAAAAGCGCTCTGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGaaa                                                                          <  2:138914/98‑1 (MQ=255)
                             tAATTAACGCCGGTAAAAAAGCGCGCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAAtt                                                                        >  2:783149/1‑100 (MQ=255)
                              aaTTAGCGCCGCTAAAAAAGCGCTCGGAATTTAGTTATGTTAAGAGAGAAAATCGTCTGTGGAATATTTCCCGTAGTCGGCTTTGTTGTTTTTGAAAttt                                                                       <  2:3826718/100‑1 (MQ=255)
                              aaTTAACGACGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAATCCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAAttt                                                                       <  2:2560311/100‑1 (MQ=255)
                                ttAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATttat                                                                     <  2:1621238/100‑1 (MQ=255)
                                 tAAAGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGATATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATttatt                                                                    <  1:2128071/100‑1 (MQ=255)
                                  aaCGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATttatta                                                                   >  1:2289260/1‑100 (MQ=255)
                                   aCGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACTGCTTTGTTGTTTTTGAAATTTATTaa                                                                  <  2:3748827/100‑1 (MQ=255)
                                   aCGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTaa                                                                  >  1:1379163/1‑100 (MQ=255)
                                     gCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAAtt                                                                <  1:391534/100‑1 (MQ=255)
                                     gCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAAtt                                                                >  1:1446020/1‑100 (MQ=255)
                                       tggTGAAAAAGCGCTCGGAAGTTAATTACTTTAAGAGATAAATCCGTGTGCGGAATATTTCTCGCAGAGGGCTTTGTTGTTTTTGAACTTTATTAAtt                                                                <  2:2642500/97‑1 (MQ=255)
                                           aaaaaaGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACa                                                          <  2:3700892/100‑1 (MQ=255)
                                           aaaaaaGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACa                                                          <  2:3202073/100‑1 (MQ=255)
                                            aaaaaGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGAAGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACaa                                                         >  1:1001440/1‑100 (MQ=255)
                                               aaGCGCTGGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAAtt                                                       <  2:3716703/99‑1 (MQ=255)
                                                aGCGTTCAGAATTTAATTATTTTAAGAGATAAAATCGGCTGCGGAATATTTGCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTa                                                      <  2:3305391/99‑1 (MQ=255)
                                                   gcTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTg                                                  >  1:1448611/1‑100 (MQ=255)
                                                   gcTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTg                                                  >  1:3165719/1‑100 (MQ=255)
                                                     tctgAATTTAATTATTTTAAGAGATAAAACCTTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTg                                                  <  1:3730016/95‑1 (MQ=255)
                                                     tCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTgag                                                <  2:2885599/100‑1 (MQ=255)
                                                     tCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACAGATTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTga                                                 >  1:2110287/1‑99 (MQ=255)
                                                      cGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCTGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTgag                                                <  1:2766779/99‑1 (MQ=255)
                                                        gAATTTAATTATTTTAAGAGATAAAACTGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGata                                             <  2:475664/100‑1 (MQ=255)
                                                        gAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGat                                              <  1:2586101/99‑1 (MQ=255)
                                                          aTTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGatata                                           >  1:2344834/1‑100 (MQ=255)
                                                          aTTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGata                                             <  1:3925517/98‑1 (MQ=255)
                                                               attattTTAAGTGTTAAAACCGTAGGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCg                                        <  2:418948/98‑1 (MQ=255)
                                                               attattTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGtt                                      <  2:1446005/100‑1 (MQ=255)
                                                                ttattTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTg                                     >  1:3306844/1‑100 (MQ=255)
                                                                 tattTTTAGAGATAAAACCGTCTGCGGAATATTTCTCGCAGACGGCTTTGTTGTTTTTGAAATTGATTAATTTAAACAATTAGTTGAGATATATCGTTgg                                    <  1:3437606/100‑1 (MQ=255)
                                                                 tattTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTgg                                    <  2:2984767/100‑1 (MQ=255)
                                                                 tattTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTgg                                    >  1:364354/1‑100 (MQ=255)
                                                                  attTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGc                                   <  2:511508/100‑1 (MQ=255)
                                                                  attTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGc                                   <  2:1521093/100‑1 (MQ=255)
                                                                   ttttAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGATGTTTTTTTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCg                                  <  1:1007699/100‑1 (MQ=255)
                                                                   ttttAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCg                                  >  1:391692/1‑100 (MQ=255)
                                                                       aaGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAAAAATTAGTTGAGATATATCGTTGGCGTcac                              >  1:604519/1‑100 (MQ=255)
                                                                        agagATAAAACCGTCTGCTGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTcaca                             <  2:1796492/100‑1 (MQ=255)
                                                                          caacaaagccGTCTGCGGGAAATATTCCGCAGACGGTTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACaaa                           >  1:2174230/9‑100 (MQ=21)
                                                                           gaTAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACaaaa                          >  1:3816559/1‑100 (MQ=255)
                                                                           aacaaagccGTCTGCGGGAAATATTCCGCAGACGGTTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACaaaa                          >  1:85253/8‑100 (MQ=21)
                                                                            acaaagccGTCTGCGGGAAATATTCCGCAGACGGTTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAg                         >  1:2131356/7‑100 (MQ=21)
                                                                            aTAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAg                         <  2:2973454/100‑1 (MQ=255)
                                                                              aaaaCCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGACACAAAAGCa                       >  1:1904952/1‑100 (MQ=255)
                                                                                 aCCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTAAGATATATAGTTGGCGTCACAAAAGAaaaa                    >  1:963800/1‑100 (MQ=255)
                                                                                   cGTCTGCGGGAAATATTCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATa                  >  1:827827/1‑100 (MQ=38)
                                                                                   cGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATa                  >  1:709968/1‑100 (MQ=255)
                                                                                   cGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAAATATATAGTTGGCGTCACAAAAGCAAAATa                  >  1:797062/1‑100 (MQ=255)
                                                                                    gTCTGCGGGAAATATTCCGCAGACGGTTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATaa                 >  1:1044377/1‑100 (MQ=38)
                                                                                     tCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAAc                <  2:1210340/100‑1 (MQ=255)
                                                                                     tCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAAc                <  2:1901208/100‑1 (MQ=255)
                                                                                      cTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATGTCGTTGGCGTCACAAAAGCAAAATAACg               <  2:2592495/100‑1 (MQ=255)
                                                                                       tGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGt              <  2:1070468/100‑1 (MQ=255)
                                                                                       tGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCAAAAAAGCAAAATAACGt              >  1:2325913/1‑100 (MQ=255)
                                                                                            aaTTTTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTc         <  2:2041019/100‑1 (MQ=255)
                                                                                                tttCCCGCAGACGGCTTTGTTTTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCTAAATAACGTAATTCGGaa     <  2:2325891/100‑1 (MQ=255)
                                                                                                   cccGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATa  >  1:765129/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATA  >  NC_000913/285105‑285303

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