Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 285,204 C→T 16.8% 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.8% 295.3 / 46.7 130intergenic (+36/‑191)yagF/yagGCP4‑6 prophage; dehydratase family protein/CP4‑6 prophage; putative sugar transporter
Reads supporting (aligned to +/- strand):  ref base C (35/73);  new base T (4/18);  total (39/91)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.05e-01

Only 100 of 140 total aligned reads displayed.

TATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATA  >  NC_000913/285105‑285303
                                                                                                   |                                                                                                   
tATTTATGACATTGATAAAATTATCGAGGTAATTAACGCCGGTAAACAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTc                                                                                                     <  1:4093427/100‑1 (MQ=255)
tATTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTc                                                                                                     >  2:3527496/1‑100 (MQ=255)
 aTTTATGACACCTATAAAATTATCGAGTTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTcc                                                                                                    <  1:891336/100‑1 (MQ=255)
 aTTTATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTcc                                                                                                    >  2:3460212/1‑100 (MQ=255)
  tttATGACACCGATAAAATTATCGAGGTAATTAACGCCGATAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTccc                                                                                                   <  2:4133133/100‑1 (MQ=255)
    tATGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGc                                                                                                 <  2:245447/100‑1 (MQ=255)
     aTGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCa                                                                                                >  2:2866702/1‑100 (MQ=255)
     aTGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCa                                                                                                >  2:32077/1‑100 (MQ=255)
      tGACACCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAg                                                                                               >  2:546916/1‑100 (MQ=255)
        acacCAATAAAATTATCGAGGTAATTAACGCCGGTAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACg                                                                                            <  2:4176335/100‑1 (MQ=255)
         cacCGATAAAATTATCGAGTTAATTATCGCCGGTAAAAAAGCGCTCGGAATTTAATTATCTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGACg                                                                                            <  1:686016/100‑1 (MQ=37)
         cacCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAGCCGTCTGCGGGAAATATTCCGCAGACg                                                                                            <  2:1449713/100‑1 (MQ=37)
         cacCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGACg                                                                                            <  1:1774496/100‑1 (MQ=38)
         cacCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACg                                                                                            <  2:3918378/100‑1 (MQ=255)
         cacCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACg                                                                                            >  2:3278747/1‑100 (MQ=255)
         cacCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACg                                                                                            >  2:3437230/1‑100 (MQ=255)
          acCGATAAAATTATCGATGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGACgg                                                                                           <  1:2480078/100‑1 (MQ=38)
          acCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGACgg                                                                                           <  2:1246764/100‑1 (MQ=38)
          acCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGACgg                                                                                           <  1:3987300/100‑1 (MQ=38)
          acCGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACgg                                                                                           >  2:2495086/1‑100 (MQ=255)
          aacGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTGAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGAc                                                                                             <  1:3107839/96‑1 (MQ=37)
           ccGATACCATTATCGAGGTAATTCACGCCGGTAAAACAGCGCTCGGAATTTAATTATTTTAAGAGATATAACCGTCTGCGGAATATTTCCCGCAGAct                                                                                            <  1:1813463/98‑2 (MQ=255)
           ccGATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGACggt                                                                                          <  2:2289226/100‑2 (MQ=38)
             gagaagaTTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCtt                                                                                        <  1:4044569/94‑1 (MQ=255)
             gATAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGACggt                                                                                          <  2:3212944/98‑2 (MQ=38)
              aTAAAATTATCGAGGTAATTAATGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGCGATAAAACCGTCTGCGGGAAATATTCCGCAGACggt                                                                                          <  2:3106006/97‑2 (MQ=37)
              aTAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCttt                                                                                       <  1:1221180/100‑1 (MQ=255)
               tAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGACGGtttta                                                                                      <  1:2600819/100‑2 (MQ=37)
               tAAAATTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTCTAAAAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCttcg                                                                                      >  2:400874/1‑98 (MQ=255)
                aaaaTTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTttgt                                                                                     <  2:2326897/100‑1 (MQ=255)
                  aaTTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTttgttg                                                                                   >  2:1721624/1‑100 (MQ=255)
                  aaTTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTttgttg                                                                                   <  1:3611581/100‑1 (MQ=255)
                   aTTATCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGGAAATATTCCGCAGACGGttttatctc                                                                                  <  1:3804098/100‑6 (MQ=25)
                    ttATCGAGGTAATTAACGCCAGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTtgttgt                                                                                  <  2:3131348/99‑1 (MQ=255)
                      aTCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGtttt                                                                               <  1:3342774/100‑1 (MQ=255)
                      aTCGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGtttt                                                                               <  1:2638257/100‑1 (MQ=255)
                        cGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTg                                                                             <  2:2263113/100‑1 (MQ=255)
                        cGAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTg                                                                             <  2:3617050/100‑1 (MQ=255)
                         gAGGTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGa                                                                            <  2:187989/100‑1 (MQ=255)
                            gTAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAAt                                                                         <  1:2320264/100‑1 (MQ=255)
                             tAATTAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTGTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAAtt                                                                        <  1:3746703/100‑1 (MQ=255)
                               aTTAACGAAAGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATtta                                                                      <  1:1693923/100‑1 (MQ=255)
                                 tAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCTGCAGACGGCTTTGTTGTTTTTGAAATttatt                                                                    <  2:2966807/100‑1 (MQ=255)
                                 tAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATttatt                                                                    >  2:37011/1‑100 (MQ=255)
                                 tAACGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATttatt                                                                    >  2:2540534/1‑100 (MQ=255)
                                  aaCGCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATttatta                                                                   <  1:2035034/100‑1 (MQ=255)
                                     gCCGGTAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAAtt                                                                <  2:2201726/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)
                                          tAAAAAAGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAAc                                                           >  2:1779013/1‑100 (MQ=255)
                                           aaaaaaGCGCTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACa                                                          <  1:4314249/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:3314161/100‑1 (MQ=255)
                                                   gcTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTg                                                  <  1:2866685/100‑1 (MQ=255)
                                                    cTCGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTga                                                 <  2:235244/100‑1 (MQ=255)
                                                      cGGAATTTAATTATTTTAAGAGATATAACATTCTGCGGAATATTGCCCGCAGACGGCATTGTTGTTTTTTAAATTTATTAATTTAAACAATTAGTTgag                                                <  1:4048668/99‑1 (MQ=255)
                                                      cGGAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTgaga                                               <  1:1569370/100‑1 (MQ=255)
                                                        gAATTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGata                                             <  1:3377222/100‑1 (MQ=255)
                                                         aaTTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGata                                             <  1:3960067/99‑1 (MQ=255)
                                                          aTTTAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGatata                                           >  2:1614331/1‑100 (MQ=255)
                                                          aTTTAATTATTTTAAGAGATAAAACCGTCTCCGGAATATTTCCCGCACACGGCTTTGTTGTTTTTAAAAATTATTAATTTAAACAATTAGTTGAGatata                                           >  2:828427/1‑100 (MQ=255)
                                                           tttAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGatatat                                          <  1:3321070/100‑1 (MQ=255)
                                                             tAATTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCg                                        >  2:3803793/1‑100 (MQ=255)
                                                              aaTTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGt                                       >  2:3906237/1‑100 (MQ=255)
                                                              aaTTATTTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGt                                       <  2:4313275/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)
                                                                  attTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGc                                   <  1:149996/100‑1 (MQ=255)
                                                                  attTTAAGAGATAAAACCGTCTGCGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGc                                   >  2:2098104/1‑100 (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)
                                                                       aaGAGATAGAACCGTCTGCGGAATATTTCCCGCAGACTGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTca                               <  1:3059066/99‑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                       <  1:930759/100‑1 (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)
                                                                                         cGGAATATTTCCCGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTaa            >  2:1766882/1‑100 (MQ=255)
                                                                                          ggAATATTTTCCATATATTGCTTTGTTGTTTTTTAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGTCGTCACAAAAGCAAAATAACGTAAt           <  1:1861846/100‑1 (MQ=255)
                                                                                            aaTATTTCACGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTc         <  1:2122540/100‑1 (MQ=255)
                                                                                                   cccGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATa  >  2:2030652/1‑100 (MQ=255)
                                                                                                   cccGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATa  >  2:2107343/1‑100 (MQ=255)
                                                                                                   cccGCAGACGGCTTTGTTGTTTTTGAAATTTATTAATTTAAACAATTAGTTGAGATATATCGTTGGCGTCACAAAAGCAAAATAACGTAATTCGGAAATa  <  1:3560060/100‑1 (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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: