Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 285,202 T→A 15.2% intergenic (+34/‑193) 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,2020TA15.2% 378.5 / 52.2 145intergenic (+34/‑193)yagF/yagGCP4‑6 prophage; dehydratase family protein/CP4‑6 prophage; putative sugar transporter
Reads supporting (aligned to +/- strand):  ref base T (48/75);  new base A (8/14);  total (56/89)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.90e-01

Only 100 of 148 total aligned reads displayed.

TGTATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAA  >  NC_000913/285103‑285300
                                                                                                   |                                                                                                  
tGTATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATAtt                                                                                                    <  1:2796487/100‑1 (MQ=255)
 gTATTTATGACACCGATAAAATTATCGAGGTAATTAACTCAGGTAAAAAAGCGCTCGGGATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATAttt                                                                                                   <  2:2122678/100‑1 (MQ=255)
  tATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTc                                                                                                  >  1:1500528/1‑100 (MQ=255)
  tATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAAACGTCTGCGGAATATTTc                                                                                                  >  2:3162156/1‑100 (MQ=255)
   aTTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTcc                                                                                                 <  2:2660661/100‑1 (MQ=255)
   aTTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTcc                                                                                                 <  2:2301603/100‑1 (MQ=255)
   aTTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTcc                                                                                                 >  1:1699978/1‑100 (MQ=255)
   aTTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTc                                                                                                  >  1:3908581/1‑99 (MQ=255)
   aTTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAAATAAAACCGTCTGCGGAATATTTcc                                                                                                 >  1:2255715/1‑100 (MQ=255)
    tttATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTccc                                                                                                <  2:1040744/100‑1 (MQ=255)
    tttATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTccc                                                                                                <  2:2493229/100‑1 (MQ=255)
    tttATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTccc                                                                                                <  2:3189426/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)
       aTGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCa                                                                                             <  2:259874/100‑1 (MQ=255)
         gACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGa                                                                                           <  2:425289/100‑1 (MQ=255)
          acacCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGAc                                                                                          <  1:2482083/100‑1 (MQ=38)
          acacCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGAc                                                                                          <  1:873738/100‑1 (MQ=255)
           cacCGATAAAATTATCGATGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAGCCGTCTGCGGGAAATATTCCGCAGACg                                                                                         <  1:413438/100‑1 (MQ=37)
           cacCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGACg                                                                                         <  2:1122125/100‑1 (MQ=38)
           cacCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACg                                                                                         >  1:147831/1‑100 (MQ=255)
           cacCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAACGCTCGGAATTTAATTATTTTAAGAGATAAAACCCTCTGAGGAATATTTCTCGCAGACg                                                                                         >  1:2080345/1‑100 (MQ=255)
            acCGATAAAATTATCGAGGTAATTAACGACGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGACg                                                                                         <  2:899520/99‑1 (MQ=38)
             ccGATAAAATTATCGATGTAATTATCTCCGGTAAAAAAACGCTCGTAATCTAATTATTTTAAGAGATACAACCGTCCGCGGCATATTTCCCGCAGACGGc                                                                                       <  2:12628/100‑1 (MQ=255)
             ccGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGc                                                                                       >  1:3182280/1‑100 (MQ=255)
              cgtcaataTTATGTAGGTAATTTACGCTGGTAGAGAGGGTCTGTGTATTTAATTTTTTTAATAGATAAAACCGTCTGCGGGAAATATTCCGCAGACg                                                                                         <  2:3617555/90‑1 (MQ=2)
              cGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGACggt                                                                                       <  2:69511/99‑2 (MQ=38)
                aTAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGCAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGACGGtttt                                                                                    <  2:2742101/100‑1 (MQ=37)
                 tAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCAttg                                                                                   >  1:3572674/1‑100 (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)
                        aTCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGtttt                                                                            <  1:3829676/100‑1 (MQ=255)
                           gAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGa                                                                         >  1:3924030/1‑100 (MQ=255)
                           gAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGa                                                                         >  1:3005075/1‑100 (MQ=255)
                           gAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGa                                                                         <  2:3247404/100‑1 (MQ=255)
                            aGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGaa                                                                        <  1:3494492/100‑1 (MQ=255)
                            aGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGaa                                                                        <  1:2672468/100‑1 (MQ=255)
                              gTAATTAATGCCGGTAAAAAAGCGTTCGGTATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAAt                                                                      <  1:317625/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)
                                 aTTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACAGCTTTGTTGTTTTTGAAATtta                                                                   >  1:799729/1‑100 (MQ=255)
                                   tAAAGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGATATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATttatt                                                                 <  1:2128071/100‑1 (MQ=255)
                                    aaCGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATttatta                                                                <  1:76440/100‑1 (MQ=255)
                                    aaCGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATttatta                                                                >  1:2973480/1‑100 (MQ=255)
                                    aaCGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATttatta                                                                <  1:3319106/100‑1 (MQ=255)
                                     aCGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACTGCTTTGTTGTTTTTGAAATTTATTaa                                                               <  2:3748827/100‑1 (MQ=255)
                                     aCGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTaa                                                               >  1:1379163/1‑100 (MQ=255)
                                       gCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAAtt                                                             >  1:1446020/1‑100 (MQ=255)
                                         cGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGCTAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTa                                                           <  2:2635866/100‑1 (MQ=255)
                                           gTAAAAAAGCGCTCTGAATTTAATTATTTTAAGATATAAAATTGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTaa                                                          <  1:3159896/99‑1 (MQ=255)
                                             aaaaaaGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACa                                                       <  2:3700892/100‑1 (MQ=255)
                                              aaaaaGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGAAGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACaa                                                      >  1:1001440/1‑100 (MQ=255)
                                                aaaGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAAtt                                                    >  1:2592519/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)
                                                      cTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTga                                              >  1:2832459/1‑100 (MQ=255)
                                                       tCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACAGATTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTga                                              >  1:2110287/1‑99 (MQ=255)
                                                        cGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCTGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTgag                                             <  1:2766779/99‑1 (MQ=255)
                                                         ggAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGat                                           <  2:2975853/100‑1 (MQ=255)
                                                          gAATTTAATTATTTTAAGAGATAAAACTGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGata                                          <  2:475664/100‑1 (MQ=255)
                                                               tAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCg                                     <  1:230877/100‑1 (MQ=255)
                                                               tAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTATTTTTGAAATATATTAATTTAAAAAATTAGTTGAGatata                                        >  1:506637/1‑97 (MQ=255)
                                                                 attattTTAAGTGTTAAAACCGTAGGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCg                                     <  2:418948/98‑1 (MQ=255)
                                                                  ttattTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTg                                  >  1:3306844/1‑100 (MQ=255)
                                                                   tattTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTgg                                 >  1:364354/1‑100 (MQ=255)
                                                                    attTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGc                                <  2:1521093/100‑1 (MQ=255)
                                                                    attTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGc                                <  2:511508/100‑1 (MQ=255)
                                                                     ttttAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGATGTTTTTTTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCg                               <  1:1007699/100‑1 (MQ=255)
                                                                     ttttAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCg                               <  2:1447983/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)
                                                                             gaTAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACaaaa                       >  1:3816559/1‑100 (MQ=255)
                                                                              acaaagccGTCTGCGGGAAATATTCCGCAGACGGTTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAg                      >  1:2131356/7‑100 (MQ=21)
                                                                              aTAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAg                      <  2:2973454/100‑1 (MQ=255)
                                                                               tAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGc                     >  1:2729075/1‑100 (MQ=255)
                                                                               tAAAACCATCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTATACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGc                     <  2:2097084/100‑1 (MQ=255)
                                                                                aaaaCCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGACACAAAAGCa                    >  1:1904952/1‑100 (MQ=255)
                                                                                     cGTCTGCGGGAAATATTCCGCAGACGGTTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTAAGATATATCATTAAAATCACAAAAGCAAAATa               >  1:853631/1‑100 (MQ=21)
                                                                                     cGTCTGCGGGAAATATTCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATa               >  1:2212236/1‑100 (MQ=38)
                                                                                     cGTCTGCGGGAAATATTCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATa               >  1:827827/1‑100 (MQ=38)
                                                                                     cGTCTGCGGGAAATATTCCGCAGACGGCTTTGATGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATa               >  1:3721473/1‑100 (MQ=38)
                                                                                     cGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATa               >  1:709968/1‑100 (MQ=255)
                                                                                      gTCTGCGGGAAATATTCCGCAGACGGTTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATaa              >  1:1044377/1‑100 (MQ=38)
                                                                                        cTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATGTCGTTGGCGTCACAAAAGCAAAATAACg            <  2:2592495/100‑1 (MQ=255)
                                                                                        cTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACg            >  1:2428118/1‑100 (MQ=255)
                                                                                         tGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGt           <  2:1070468/100‑1 (MQ=255)
                                                                                         tGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCAAAAAAGCAAAATAACGt           >  1:2325913/1‑100 (MQ=255)
                                                                                              aaTTTTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTc      <  2:2041019/100‑1 (MQ=255)
                                                                                                 atTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCgg    <  2:1689823/99‑1 (MQ=255)
                                                                                                  tttCCCGCAGACGGCTTTGTTTTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCTAAATAACGTAATTCGGaa  <  2:2325891/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
TGTATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAA  >  NC_000913/285103‑285300

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: