Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,472,437 T→C 16.1% 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,4370TC16.1% 408.6 / 69.5 191N33D (AAT→GAT) pyrLpyrBI operon leader peptide
Reads supporting (aligned to +/- strand):  ref base T (75/83);  new base C (19/11);  total (96/95)
Fisher's exact test for biased strand distribution p-value = 1.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.13e-02

Only 100 of 192 total aligned reads displayed.

ATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACA  >  NC_000913/4472339‑4472536
                                                                                                  |                                                                                                   
aTCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCTGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTg                                                                                                    <  1:431612/100‑1 (MQ=255)
 tCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGa                                                                                                   >  1:2154845/1‑100 (MQ=255)
   gcgACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAgg                                                                                                 <  2:1393550/100‑1 (MQ=255)
   gcgACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAgg                                                                                                 >  2:649550/1‑100 (MQ=255)
      aCTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCAATCTTTTATCTCCTGACGCCTGGGAAAAAAAAGCCCATCGATTGAggg                                                                                                >  2:3393939/1‑97 (MQ=255)
         aaGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTgg                                                                                           >  2:1106898/1‑100 (MQ=255)
          aGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTggg                                                                                          <  2:3655598/100‑1 (MQ=255)
          aGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCAATCGAGGGGCtttt                                                                                          <  2:1567586/100‑4 (MQ=255)
          aGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGACCCTCGATTGAGGGGCTggg                                                                                          >  1:918082/1‑100 (MQ=255)
           ggTCTTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCATGGGCAAAAAAAAGCCCCTCAATCGAGGGGCtttt                                                                                          <  1:1537232/99‑4 (MQ=255)
           ggTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAACGCCCCTCGATTGAGGGGCTGGGa                                                                                         <  2:1219104/100‑1 (MQ=255)
               gtgtATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATTTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGaa                                                                                        <  1:1666266/94‑1 (MQ=255)
               gTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAAt                                                                                       <  2:2457948/98‑1 (MQ=255)
                tttATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATggg                                                                                    <  1:1131037/100‑1 (MQ=255)
                tttATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATggg                                                                                    <  1:3617693/100‑1 (MQ=255)
                tttATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCATGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATgg                                                                                     <  1:483643/99‑1 (MQ=255)
                 ttATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGt                                                                                   <  2:1247202/100‑1 (MQ=255)
                   tctgAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTg                                                                                  <  2:1443084/96‑1 (MQ=255)
                   aTGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGa                                                                                 <  2:3912358/100‑1 (MQ=255)
                         aTGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACg                                                                           <  2:270299/100‑1 (MQ=255)
                          tGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACgg                                                                          <  2:2862037/100‑1 (MQ=255)
                          tGATATGTTTCTGATATAGCGGATTAGCCATCTTTTACCTCCTGACGCCTGTGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGAGCAACgg                                                                          <  2:778177/100‑1 (MQ=255)
                            atatGTTTCTGACATAGCCGATTAGCCAACTTTTATCTCCTGACGATTGGGCAAAAAAAAGCCCCCCGCTTGAGGCGCTGGGAATGGGTGATCAACGGga                                                                        >  2:2161258/1‑100 (MQ=255)
                             tatTTTTCTGATATAGCGGATTAGCCATCTTTTATCTTCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGgan                                                                       <  1:1295098/100‑2 (MQ=255)
                             tatGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGgaa                                                                       >  1:2641444/1‑100 (MQ=255)
                               tGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGgaaga                                                                     >  2:3712849/1‑100 (MQ=255)
                               tGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGgaaga                                                                     <  1:1254794/100‑1 (MQ=255)
                                gTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGgaagaa                                                                    >  2:858096/1‑100 (MQ=255)
                                gTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACACCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGga                                                                        >  2:1792753/1‑96 (MQ=255)
                                 tttCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGaaa                                                                   >  2:1324082/1‑100 (MQ=255)
                                 tttCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGaaa                                                                   >  1:547779/1‑100 (MQ=255)
                                  ttCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGAATGAGGGAAAGGGAATGGGTGATCAACGGGAAGaaaa                                                                  >  1:1149835/1‑100 (MQ=255)
                                     tGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAAc                                                                <  1:3819224/99‑1 (MQ=255)
                                     tGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACg                                                               >  2:3931442/1‑100 (MQ=255)
                                       atatAGCTGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGc                                                             <  1:1689150/100‑1 (MQ=255)
                                       atatAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGc                                                             <  1:3844679/100‑1 (MQ=255)
                                        tataGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCa                                                            <  1:3926476/100‑1 (MQ=255)
                                        tataGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATCGGTGATCAACGGGAAGAAAAACGGCa                                                            <  1:2384847/100‑1 (MQ=255)
                                            gCGGATTAGCTATTTTTTATCTTCTGACGTCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGcc                                                        <  1:1551553/100‑1 (MQ=255)
                                            gCGGATTAGCCATCTTTTATCTCTTGACGCCTTGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGcc                                                        <  1:3222453/100‑1 (MQ=255)
                                              ggATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAg                                                      >  1:596510/1‑100 (MQ=255)
                                                aTTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAg                                                      <  1:1324096/98‑1 (MQ=255)
                                                   aGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCt                                                 >  2:487036/1‑100 (MQ=255)
                                                         cTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTc                                           <  2:2802376/100‑1 (MQ=255)
                                                         cTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTc                                           >  2:70189/1‑100 (MQ=255)
                                                            ttATCTCCTGACGTCTGTGCAAAAAAATGTTTTTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCtgaaaaaaga                                        <  1:179737/100‑10 (MQ=255)
                                                             tATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGAc                                       <  2:107500/100‑1 (MQ=255)
                                                             tATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGAc                                       >  1:353054/1‑100 (MQ=255)
                                                                ctcCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGAcgcg                                    <  1:649558/100‑1 (MQ=255)
                                                                ctcCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGAcgcg                                    <  2:2641420/100‑1 (MQ=255)
                                                                 tcCTGACGCCTGGGCAAAAAAAAGCCACTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCATAcgc                                     >  2:939920/1‑98 (MQ=255)
                                                                  ccTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGt                                  >  2:1803296/1‑100 (MQ=255)
                                                                   cTGACGCCTGGGGAAAGAAAAGCGCTTCGATGGAGGGGCGGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTa                                 <  2:3311208/100‑1 (MQ=255)
                                                                     gACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAg                               >  1:3067289/1‑100 (MQ=255)
                                                                       cGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGAc                             >  2:2876303/1‑100 (MQ=255)
                                                                        gCCTGGGAAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGAc                             <  2:425401/99‑1 (MQ=255)
                                                                          cTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACaaa                          >  1:518725/1‑100 (MQ=255)
                                                                           tGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACaaaa                         <  2:251014/100‑1 (MQ=255)
                                                                            gggCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAAt                        <  2:3199622/100‑1 (MQ=255)
                                                                              tcAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGt                      <  1:1922254/99‑1 (MQ=255)
                                                                              gCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGt                      <  1:3361109/100‑1 (MQ=255)
                                                                               cAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAAt                        >  1:1530665/1‑97 (MQ=255)
                                                                               cAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATg                       <  1:3797668/98‑1 (MQ=255)
                                                                               cAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGta                     >  2:1922238/1‑99 (MQ=255)
                                                                               cAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTc                     >  1:2937546/1‑100 (MQ=255)
                                                                               cAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAAAGGCAGGGCAACGGCTTTTTTCAGAAGCCGTAaaaaaaaa                         >  1:3743407/1‑96 (MQ=255)
                                                                                aaaaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCg                    <  1:1727020/100‑1 (MQ=255)
                                                                                 ttccgatctcccTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGa                   <  2:251566/91‑1 (MQ=255)
                                                                                 aaaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGa                   <  1:2075438/100‑1 (MQ=255)
                                                                                 aaaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGa                   >  2:3733250/1‑100 (MQ=255)
                                                                                  tctaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGaa                  <  2:3701625/97‑1 (MQ=255)
                                                                                   aaaaaGCTCTTTGTTTGAGGGGCTGGGAATTGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAac                 <  1:3246825/100‑1 (MQ=255)
                                                                                   aaaaaGCCCCTCGATTGTGGGGCTGGGAATTTGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAac                 <  1:3713490/100‑1 (MQ=255)
                                                                                   aaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAac                 <  1:2321220/100‑1 (MQ=255)
                                                                                   aaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAac                 >  2:3115167/1‑100 (MQ=255)
                                                                                   aaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAac                 >  1:3629558/1‑100 (MQ=255)
                                                                                    cccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAaca                >  1:2789822/4‑100 (MQ=255)
                                                                                    cccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAaca                >  2:2899148/4‑100 (MQ=255)
                                                                                    cccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAaca                >  2:3461426/4‑100 (MQ=255)
                                                                                    cccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAAACAAAATGTCGAaca                >  1:2723529/4‑100 (MQ=255)
                                                                                    cccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAAACGCGGTAAGACAAAATGTCGAaca                >  1:1374443/4‑100 (MQ=255)
                                                                                    cccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAAAGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTAGAaca                >  2:1264522/4‑100 (MQ=255)
                                                                                    aaaaGCACCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAac                 <  2:2023517/99‑1 (MQ=255)
                                                                                     ccaGCCCCTCAATCGAGGGGGTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacac               >  2:2381555/3‑100 (MQ=255)
                                                                                     ccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacac               >  2:1891670/3‑100 (MQ=255)
                                                                                     ccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCCAacac               >  1:2342999/3‑100 (MQ=255)
                                                                                     ccaGCCCCTCAATCGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACACGGTAAAACAAAATGTCAAaaac               >  1:281028/3‑100 (MQ=255)
                                                                                     ccaGCCCCTCAATAGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCGTTTTTCAGACGCGGTAAGACAAAATGTCGAaca                >  2:3522269/3‑99 (MQ=255)
                                                                                     aaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGaaaga               >  1:3701664/1‑97 (MQ=255)
                                                                                      cgatctcTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacaca              <  2:2993465/94‑1 (MQ=255)
                                                                                       aGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacacac             <  2:2724564/100‑1 (MQ=255)
                                                                                          cccTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGAAAAAATGTCGAACACACTGa          >  1:1284712/1‑100 (MQ=255)
                                                                                           ccTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACTCGGTAAGACAAAATGTCGAACACACTGaa         >  1:1424293/1‑100 (MQ=255)
                                                                                           ccTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGaa         >  2:29371/1‑100 (MQ=255)
                                                                                           ccTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGaa         >  1:2519277/1‑100 (MQ=255)
                                                                                            cTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAacacaagatcg        >  1:2993492/1‑94 (MQ=255)
                                                                                             tCGATTGAGGGGCTGGGAATTGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAAcc       <  2:1845210/100‑1 (MQ=255)
                                                                                              cGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAAcc       <  2:2063596/99‑1 (MQ=255)
                                                                                               gATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCAt     >  2:3817783/1‑100 (MQ=255)
                                                                                                  tGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACa  <  1:168483/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
ATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACA  >  NC_000913/4472339‑4472536

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: