Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,472,437 T→C 17.6% N33D (AAT→GAT)  pyrL ← pyrBI operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,472,4370TC17.6% 419.5 / 82.0 193N33D (AAT→GAT) pyrLpyrBI operon leader peptide
Reads supporting (aligned to +/- strand):  ref base T (87/72);  new base C (25/9);  total (112/81)
Fisher's exact test for biased strand distribution p-value = 5.53e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.63e-01

Only 100 of 202 total aligned reads displayed.

CATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACA  >  NC_000913/4472338‑4472536
                                                                                                    |                                                                                                   
cATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGAtt                                                                                                     <  2:3969045/100‑1 (MQ=255)
    gcgACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAgg                                                                                                 >  2:1929043/1‑100 (MQ=255)
    gcgACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAgg                                                                                                 >  1:3141287/1‑100 (MQ=255)
    gcgACTAAGGTAGTGTATGGAAATGGTATGCTTCTGATATAGCGGTTTAGTCATCTTTTATCTCCTGAGACCTAAGCAAAAAAAAGCCCCTCGATTGAgg                                                                                                 <  1:975031/100‑1 (MQ=255)
     cgACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAggg                                                                                                >  1:101532/1‑100 (MQ=255)
     cgACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAggg                                                                                                <  1:1744089/100‑1 (MQ=255)
     cgACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAggg                                                                                                >  1:466473/1‑100 (MQ=255)
     cgACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAggg                                                                                                <  2:4027053/100‑1 (MQ=255)
     cgACTAAGGTAGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCGGGGCAAAAAAAAGCCCCTCGATTGAgg                                                                                                 <  1:1882054/99‑1 (MQ=255)
       aCTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGc                                                                                              >  2:1601404/1‑100 (MQ=255)
       aCTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGc                                                                                              >  1:743370/1‑100 (MQ=255)
        cTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCt                                                                                             >  2:696570/1‑100 (MQ=255)
        cTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAGCCCCTCAATCGAGGGGCtt                                                                                            <  2:3176163/100‑2 (MQ=255)
          aaGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCAATCGAGGGGCttt                                                                                           <  1:1560276/100‑3 (MQ=255)
           aGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTggg                                                                                          <  1:670986/100‑1 (MQ=255)
           aGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTggg                                                                                          <  1:2959390/100‑1 (MQ=255)
              tCGTTTATGGAAATGATATGTTTCTGATATAGCGGAGTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAAGCCCCTCAATCGAGGGGCtttttt                                                                                        <  2:3101762/100‑6 (MQ=255)
               cGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCAATCGAGGGGCttttttt                                                                                       <  2:4199959/99‑7 (MQ=39)
                gTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCAATCGAGGGGCttttttttg                                                                                     <  2:827396/100‑9 (MQ=38)
                 tttATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAACAAGCCCCTCGATTGAGCGGCTGGGAAt                                                                                       >  1:4211900/1‑97 (MQ=255)
                 tttATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATggg                                                                                    <  2:66836/100‑1 (MQ=255)
                 tttATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATggg                                                                                    <  1:583048/100‑1 (MQ=255)
                   tacggAAATCATATGTTTCTGATATAGCGAATTAGCCATCTTTTATCTAAAAAAAAAAAAAAAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTg                                                                                  <  1:3683589/97‑1 (MQ=255)
                   tATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTg                                                                                  <  1:1351573/100‑1 (MQ=255)
                    aaaggaaTGATATGTTTATGATATAGCGGAATAGCCAGCTTGTATCTCCTGACGCAGGAGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTg                                                                                  <  1:2300073/94‑1 (MQ=255)
                    aTGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGa                                                                                 <  2:3746941/100‑1 (MQ=255)
                     tGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAGAGCCCCTCGATTGAGGGGCTGGGAATGGGTGAt                                                                                >  2:3913670/1‑100 (MQ=255)
                      ggAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATc                                                                               >  1:349830/1‑100 (MQ=255)
                      ggAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATc                                                                               <  2:1996016/100‑1 (MQ=255)
                      ggAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATc                                                                               >  2:2536642/1‑100 (MQ=255)
                       gAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCa                                                                              >  1:84796/1‑100 (MQ=255)
                        aaaTGTTATGATTCTGATATAGCGGAGGAGCCATCTTTTATCTCCTGACGCCTGTGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCa                                                                              <  1:1196509/99‑1 (MQ=255)
                          aTGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACg                                                                           >  2:3715565/1‑100 (MQ=255)
                          aTGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACg                                                                           <  2:3467118/100‑1 (MQ=255)
                            gATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACggn                                                                         <  2:1045276/100‑2 (MQ=255)
                                tGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCCCGATTGAGGGGCTGGGAATGGGTGATCAACGGgaaga                                                                     >  2:2415554/1‑100 (MQ=255)
                                  tttCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGaaa                                                                   >  2:2458654/1‑100 (MQ=255)
                                  tttCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGaaa                                                                   >  2:3809873/1‑100 (MQ=255)
                                       gATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACgg                                                              >  2:169285/1‑100 (MQ=255)
                                         tataGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCa                                                            >  2:1309406/1‑100 (MQ=255)
                                          ataGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAg                                                           >  2:4224690/1‑100 (MQ=255)
                                          ataGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAg                                                           <  2:925499/100‑1 (MQ=255)
                                           taGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAgg                                                          <  2:2405404/100‑1 (MQ=255)
                                            aGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGc                                                         <  2:2949075/100‑1 (MQ=255)
                                            aGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGc                                                         >  1:3590358/1‑100 (MQ=255)
                                            aGCGGATGAGCCATTTTTTATTTCCTGAGGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGc                                                         <  1:2551679/100‑1 (MQ=255)
                                              cGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCa                                                       >  2:1631061/1‑100 (MQ=255)
                                                gATTAGCCATCTTTTATCTGCTGACGAATGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGc                                                         <  1:537690/96‑1 (MQ=255)
                                                 aTTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCg                                                    <  2:3780352/100‑1 (MQ=255)
                                                     gCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCTCTTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAACACGGCAGGCCAGCGTCtt                                                >  1:3437592/1‑100 (MQ=255)
                                                     gCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCtt                                                >  2:3374557/1‑100 (MQ=255)
                                                     gCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCtt                                                >  1:2459590/1‑100 (MQ=255)
                                                       cATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCtttt                                              >  2:2371131/1‑100 (MQ=255)
                                                          cTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTc                                           >  2:1099287/1‑100 (MQ=255)
                                                          cTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTc                                           >  2:4082694/1‑100 (MQ=255)
                                                          cTTTTATCTCCTGACGCCTGAGCAAAAAAAAACCCATCGATTGAGGGGCTGGCAATGGGTAAAAAAAAAAAAAAAAAACGGCAGGCCAGCGTCTTTTTTc                                           <  1:4053632/100‑1 (MQ=255)
                                                             ttATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGa                                        >  2:1521274/1‑100 (MQ=255)
                                                               aTCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGAcg                                      >  2:1170859/1‑100 (MQ=255)
                                                               aTCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGAcg                                      <  2:953745/100‑1 (MQ=255)
                                                                 ctcCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGAcgcg                                    <  2:2204333/100‑1 (MQ=255)
                                                                  tcCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCgg                                   >  2:268770/1‑100 (MQ=255)
                                                                    cggACGTCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTa                                 <  1:1228263/98‑1 (MQ=255)
                                                                    cTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGAcgctggc                                 <  2:349831/100‑5 (MQ=255)
                                                                    cTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTa                                 <  2:1199936/100‑1 (MQ=255)
                                                                      gACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAg                               <  2:159022/100‑1 (MQ=255)
                                                                       aCGCCTAGGCAAAAAAAAGCCCCTCGATTGAGGCGCTGGTACTGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTac                                <  1:2734248/98‑2 (MQ=255)
                                                                        cGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGAc                             >  1:2998120/1‑100 (MQ=255)
                                                                           cTGGGCAAAAAAAAGCCCCTCGAGTGAGGGGCCAGGCAGACGTCATCCACGGGCAGAACAACGGCAGGTCAGCGGCATGTTTCAGACGCGGTAAGACa                            >  1:2060953/1‑98 (MQ=255)
                                                                              ggCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATg                       >  2:3165012/1‑100 (MQ=255)
                                                                               gCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGt                      >  2:4099150/1‑100 (MQ=255)
                                                                                cccattcccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATg                       >  1:4276566/10‑99 (MQ=38)
                                                                                cAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTc                     >  2:3990063/1‑100 (MQ=255)
                                                                                cAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTc                     >  2:275198/1‑100 (MQ=255)
                                                                                cAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTc                     >  2:177151/1‑100 (MQ=255)
                                                                                  aaaaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCg                    <  2:213100/100‑1 (MQ=255)
                                                                                   attcccaGCCCCTCAATCGAGGGGCTGGGAATGGATGATCAACGAGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGt                      >  1:512759/7‑97 (MQ=38)
                                                                                   aaaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGa                   >  1:1167169/1‑100 (MQ=255)
                                                                                   aaaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGa                   >  2:2161482/1‑100 (MQ=255)
                                                                                     tcccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAac                 >  1:3806337/5‑100 (MQ=255)
                                                                                     aaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAac                 <  2:1500277/100‑1 (MQ=255)
                                                                                      cccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAaca                >  2:131836/4‑100 (MQ=255)
                                                                                      cccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAaca                >  1:2635748/4‑100 (MQ=255)
                                                                                       ccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCATACGTGGTAAGACAAAATGTCGAacac               >  1:2216465/3‑100 (MQ=255)
                                                                                       ccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacac               >  2:2851851/3‑100 (MQ=255)
                                                                                       ccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacac               >  2:495376/3‑100 (MQ=255)
                                                                                       ccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacac               >  1:3988805/3‑100 (MQ=255)
                                                                                       aaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacac               <  1:1479194/100‑1 (MQ=255)
                                                                                        caGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacac               >  2:485673/2‑99 (MQ=255)
                                                                                        aaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacaca              >  1:1716005/1‑100 (MQ=255)
                                                                                         aGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacacac             >  2:311893/1‑100 (MQ=255)
                                                                                           ccccTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTg           <  1:70038/100‑1 (MQ=255)
                                                                                           ccccTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTg           >  1:2616910/1‑100 (MQ=255)
                                                                                           ccccTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTg           <  1:936899/100‑1 (MQ=255)
                                                                                           ccccTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTg           >  2:3327390/1‑100 (MQ=255)
                                                                                             ccTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTTTTACTTTTTAAGACAAAATGTCGAACACACTGaa         >  1:431352/1‑100 (MQ=255)
                                                                                             ccTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGaa         >  1:2791410/1‑100 (MQ=255)
                                                                                             ccTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGaa         >  2:455860/1‑100 (MQ=255)
                                                                                              cTCGATTGAGGGGCTGGGACTGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTCTTCAGACGCGGTAAGACAAAATGTCGAACACACTGacc        <  1:1468079/100‑3 (MQ=255)
                                                                                               tCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacagatcggaag       >  2:720073/1‑91 (MQ=255)
                                                                                                    tGATGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACa  <  1:3073169/100‑1 (MQ=255)
                                                                                                    |                                                                                                   
CATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACA  >  NC_000913/4472338‑4472536

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: