Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 285,202 T→A 15.3% 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.3% 352.8 / 49.2 137intergenic (+34/‑193)yagF/yagGCP4‑6 prophage; dehydratase family protein/CP4‑6 prophage; putative sugar transporter
Reads supporting (aligned to +/- strand):  ref base T (37/79);  new base A (4/17);  total (41/96)
Fisher's exact test for biased strand distribution p-value = 3.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.92e-01

Only 100 of 139 total aligned reads displayed.

GTATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTC  >  NC_000913/285104‑285296
                                                                                                  |                                                                                              
gTATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATAttt                                                                                               >  2:3920466/1‑100 (MQ=255)
 tATTTATGACATTGATAAAATTATCGAGGTAATTAACGCCGGTAAACAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTc                                                                                              <  1:4093427/100‑1 (MQ=255)
 tATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTc                                                                                              >  2:3527496/1‑100 (MQ=255)
 tATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTc                                                                                              <  1:4332220/100‑1 (MQ=255)
  aTTTATGACACCTATAAAATTATCGAGTTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTcc                                                                                             <  1:891336/100‑1 (MQ=255)
   tttATGACACCGATAAAATTATCGAGGTAATTAACGCCGATAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTccc                                                                                            <  2:4133133/100‑1 (MQ=255)
     tATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGc                                                                                          <  2:245447/100‑1 (MQ=255)
      aTGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCa                                                                                         >  2:32077/1‑100 (MQ=255)
       tGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAg                                                                                        >  2:546916/1‑100 (MQ=255)
         acacCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGAc                                                                                      <  1:2438554/100‑1 (MQ=255)
         acacCAATAAAATTATCGAGGTAATTAACGCCGGTAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACg                                                                                     <  2:4176335/100‑1 (MQ=255)
          cacCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGACg                                                                                     <  1:1774496/100‑1 (MQ=38)
          cacCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACg                                                                                     <  2:1514020/100‑1 (MQ=255)
          cacCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACg                                                                                     >  2:3278747/1‑100 (MQ=255)
          cacCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACg                                                                                     >  2:3437230/1‑100 (MQ=255)
          cacCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACg                                                                                     <  2:3918378/100‑1 (MQ=255)
           acCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGACgg                                                                                    <  2:1733920/100‑1 (MQ=38)
           aacGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTGAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGAc                                                                                      <  1:3107839/96‑1 (MQ=37)
            ccGATACCATTATCGAGGTAATTCACGCCGGTAAAACAGCGCTCGGAATTTAATTATTTTAAGAGATATAACCGTCTGCGGAATATTTCCCGCAGAct                                                                                     <  1:1813463/98‑2 (MQ=255)
            ccGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTCATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGc                                                                                          <  1:3647645/93‑1 (MQ=255)
            ccGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGACggt                                                                                   <  2:2289226/100‑2 (MQ=38)
             cGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGACggtt                                                                                  <  1:3926491/100‑3 (MQ=38)
             cGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCt                                                                                  <  2:3455475/100‑1 (MQ=255)
              gagaagaTTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCtt                                                                                 <  1:4044569/94‑1 (MQ=255)
              gATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGACGGttt                                                                                 <  2:2290942/100‑1 (MQ=37)
               aTAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCttt                                                                                <  1:1221180/100‑1 (MQ=255)
                tAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGACGGtttta                                                                               <  1:2600819/100‑2 (MQ=37)
                tAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTCTAAAAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCttcg                                                                               >  2:400874/1‑98 (MQ=255)
                    aTTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGACGGttttatctc                                                                           <  1:3804098/100‑6 (MQ=25)
                     ttATCGAGGTAATTAACGCCAGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTtgttgt                                                                           <  2:3131348/99‑1 (MQ=255)
                       aTCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGtttt                                                                        <  1:2638257/100‑1 (MQ=255)
                        tCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGttttt                                                                       <  2:3060420/100‑1 (MQ=255)
                        tCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGttttt                                                                       <  2:3801855/100‑1 (MQ=255)
                         cGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTg                                                                      <  2:3617050/100‑1 (MQ=255)
                         cGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTg                                                                      <  2:2263113/100‑1 (MQ=255)
                          gAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGa                                                                     <  2:187989/100‑1 (MQ=255)
                           aGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGaa                                                                    <  1:3016523/100‑1 (MQ=255)
                            ggTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGaaa                                                                   >  2:2361334/1‑100 (MQ=255)
                             gTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAAt                                                                  <  1:2320264/100‑1 (MQ=255)
                                 ttAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATttat                                                              <  2:2692597/100‑1 (MQ=255)
                                  tAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCTGCAGACGGCTTTGTTGTTTTTGAAATttatt                                                             <  2:2966807/100‑1 (MQ=255)
                                   aaCGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATttatta                                                            <  1:2035034/100‑1 (MQ=255)
                                         ggTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTaa                                                      <  2:1381340/100‑1 (MQ=255)
                                          gTAAAAAAGTGCTCGGAATTTAATTATTTTAAGAGATAAATCCTTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTaaa                                                     <  1:1101251/100‑1 (MQ=255)
                                           tAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAAc                                                    <  1:506638/100‑1 (MQ=255)
                                             aaaaaGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACaa                                                  >  2:3751503/1‑100 (MQ=255)
                                               aaaGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAAtt                                                <  1:1721636/100‑1 (MQ=255)
                                                aaGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTa                                               >  2:1553520/1‑100 (MQ=255)
                                                  gcgcTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGt                                             <  1:2951362/100‑1 (MQ=255)
                                                   cgcTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGtt                                            <  2:2835116/100‑1 (MQ=255)
                                                    gcTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACCATCAGTTg                                           >  2:4095529/1‑100 (MQ=255)
                                                    gcTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTg                                           <  1:2866685/100‑1 (MQ=255)
                                                     cTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTga                                          <  2:2531085/100‑1 (MQ=255)
                                                     cTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTga                                          >  2:3320710/1‑100 (MQ=255)
                                                     cTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTga                                          <  2:235244/100‑1 (MQ=255)
                                                       cGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTgaga                                        <  1:1569370/100‑1 (MQ=255)
                                                        ggAATTTAATTATTTTAAGAGATAATACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGATATTTATTAATTTAAACAATTAGTTGAGat                                       <  1:1816275/100‑1 (MQ=255)
                                                         gAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGata                                      <  1:3377222/100‑1 (MQ=255)
                                                         gAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGata                                      <  1:4031487/100‑1 (MQ=255)
                                                          aaTTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGata                                      <  1:3960067/99‑1 (MQ=255)
                                                           aTTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGatata                                    >  2:1614331/1‑100 (MQ=255)
                                                           aTTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGatata                                    >  2:1844119/1‑100 (MQ=255)
                                                           aTTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGatata                                    >  2:362193/1‑100 (MQ=255)
                                                           aTTTAATTATTTTAAGAGATAAAACCGTCTCCGGAATATTTCCCGCACACGGCTTTGTTGTTTTTAAAAATTATTAATTTAAACAATTAGTTGAGatata                                    >  2:828427/1‑100 (MQ=255)
                                                            tttAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGatatat                                   <  1:3321070/100‑1 (MQ=255)
                                                            tttAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGatatat                                   <  1:3321261/100‑1 (MQ=255)
                                                            tttAATTATTTCAAGAGATAAAACAGTCTGCGGAATATTTCCCGCAGCCGGCTTTGTTCTCCTTGAAACTTATTACTTTAAACAATTAGTTGAGataca                                    <  1:819293/99‑3 (MQ=255)
                                                               aaTTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGt                                >  2:3906237/1‑100 (MQ=255)
                                                               aaTTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGt                                <  2:4313275/100‑1 (MQ=255)
                                                                attattTTAAGAGATAAAAGCGTCTGCGGAATATTTCCCGCAGACGGCGTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTGGAGATATATCGtt                               <  1:3567/100‑1 (MQ=255)
                                                                attattTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGtt                               >  2:338794/1‑100 (MQ=255)
                                                                 ttattTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGtt                               >  2:931265/1‑99 (MQ=255)
                                                                 ttattTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTg                              <  2:2650279/100‑1 (MQ=255)
                                                                 ttattTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTg                              <  1:546914/100‑1 (MQ=255)
                                                                  tattTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTgg                             <  1:3177046/100‑1 (MQ=255)
                                                                   attTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGc                            >  2:2098104/1‑100 (MQ=255)
                                                                   attTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGc                            <  1:149996/100‑1 (MQ=255)
                                                                    ttttAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACCATTAGTTGAGATATATCGTTGGCg                           >  2:1008303/1‑100 (MQ=255)
                                                                    ttttAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCg                           <  1:4098531/100‑1 (MQ=255)
                                                                    ttttAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCg                           <  1:3600876/100‑1 (MQ=255)
                                                                     tttAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGt                          <  2:1165110/100‑1 (MQ=255)
                                                                      ttAAGAGCTAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTc                         <  1:3711530/100‑1 (MQ=255)
                                                                      ttAAGAGATAAAACCGTCTGTGGAATATTTCACGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTc                         <  1:3145365/100‑1 (MQ=255)
                                                                         agagATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTcaca                      <  2:1889130/100‑1 (MQ=255)
                                                                         agagATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTcaca                      >  2:1125717/1‑100 (MQ=255)
                                                                           caacaaagccGTCTGCGGGAAATATTCCGCAGACGGTTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACaaa                    >  2:72243/9‑100 (MQ=21)
                                                                              ccgatctgTCTGCGGGAAATATTCCGCAGACGGTTTTTTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAg                  <  1:4215814/92‑1 (MQ=25)
                                                                              caaagCCGTCTGCGGGAAATATTCCGCAGACGGATTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGc                 >  2:250640/2‑100 (MQ=25)
                                                                               aatacagTCTGTGGAATATTTCCCGCAGACGGCGTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAg                  <  2:4302522/92‑1 (MQ=255)
                                                                               aaaaCCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCa                >  2:1266966/1‑100 (MQ=255)
                                                                                aagccGTCTGCGGGAAATATTCCGCAGACGGTTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCaa               >  2:568196/4‑100 (MQ=25)
                                                                                  aCCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACATTTAGTTGAGATATATCGTTGGCGTCACAAAAGCaaaa             <  1:2361324/100‑1 (MQ=255)
                                                                                    cGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATa           <  1:1742513/100‑1 (MQ=255)
                                                                                     gTCTGCGGGAAATATTCCGCAGACGGTTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAgagatcgga          >  2:4215845/1‑92 (MQ=25)
                                                                                        tGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGt       <  1:1553514/100‑1 (MQ=255)
                                                                                         gCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTa      >  2:2170215/1‑100 (MQ=255)
                                                                                          cGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTaa     >  2:1766882/1‑100 (MQ=255)
                                                                                            gAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAAtt   >  2:3170624/1‑100 (MQ=255)
                                                                                            gAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAAt    >  1:16594/1‑99 (MQ=255)
                                                                                             aaTATTTCACGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTc  <  1:2122540/100‑1 (MQ=255)
                                                                                                  |                                                                                              
GTATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTC  >  NC_000913/285104‑285296

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