Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,950,966:1 +C 17.3% coding (135/567 nt) yecD → isochorismatase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,950,9651.C17.3% 470.7 / 113.3 150coding (134/567 nt)yecDisochorismatase family protein
Reads supporting (aligned to +/- strand):  ref base . (59/65);  new base C (18/8);  total (77/73)
Fisher's exact test for biased strand distribution p-value = 5.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.54e-02

Only 100 of 156 total aligned reads displayed.

TGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAG  >  NC_000913/1950869‑1951063
                                                                                                 |                                                                                                  
tgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGCCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCCGCGAGATATCg‑cgc                                                                                                 >  2:3667567/1‑100 (MQ=255)
 gATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGcc                                                                                                >  1:545161/1‑100 (MQ=255)
 gATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGcc                                                                                                <  1:3349048/100‑1 (MQ=255)
    ngATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGc                                                                                             >  1:1695236/2‑100 (MQ=255)
    cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGC‑GCGAAATTTCGCCGCCAGc                                                                                             <  2:283856/100‑1 (MQ=255)
    cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGC‑GCGAAATTTCGCCGCCAGc                                                                                             <  1:496529/100‑1 (MQ=255)
    cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGC‑GCGAAATTTCGCCGCCAGc                                                                                             <  1:446794/100‑1 (MQ=255)
    cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGC‑GCGAAATTTCGCCGCCAGc                                                                                             <  1:3585111/100‑1 (MQ=255)
    cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGC‑GCGAAATTTCGCCGCCAGc                                                                                             <  2:735125/100‑1 (MQ=255)
    cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGC‑GCGAAATTTCGCCGCCAGc                                                                                             <  1:2090255/100‑1 (MQ=255)
     gATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCg                                                                                            >  2:1233594/1‑100 (MQ=255)
        ttACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTc                                                                                         <  2:3722358/100‑1 (MQ=255)
         tACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCa                                                                                        >  1:3835369/1‑100 (MQ=255)
           cAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGTGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGc                                                                                      >  1:190163/1‑100 (MQ=255)
           cAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGc                                                                                      >  1:1877481/1‑100 (MQ=255)
           cAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGc                                                                                      >  1:4189761/1‑100 (MQ=255)
           cAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGc                                                                                      >  1:1019702/1‑100 (MQ=255)
           cAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGC‑GCGAAATTTCGCCGCCAGcttcccgg                                                                                      <  1:1590239/100‑8 (MQ=37)
           cAAGAAGGCAGCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGc                                                                                      <  1:178285/100‑1 (MQ=255)
            aagaagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGcc                                                                                     >  2:3406579/1‑100 (MQ=255)
            aagaagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGcc                                                                                     <  1:1370575/100‑1 (MQ=255)
             agaagGCATCTTACCTTTTGGCGGTGGTCCACATACTGCCGATTAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGcc                                                                                     <  2:3963047/99‑1 (MQ=255)
              taacgCATCGTACCTTTTGCCGGTTGTCCACATACTTCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGccct                                                                                   <  2:4034852/96‑2 (MQ=255)
                 ggCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCgtgt                                                                                <  1:2206023/100‑1 (MQ=255)
                 ggCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCgtgt                                                                                <  2:2786764/100‑1 (MQ=255)
                  gCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGtt                                                                               <  1:318752/100‑1 (MQ=255)
                    aTCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTc                                                                             >  2:2252810/1‑100 (MQ=255)
                    aTCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTc                                                                             <  2:3196221/100‑1 (MQ=255)
                        tACCTTTTGCCGGTGGTCCACATACCGCCATCGAGGTAGTTAATCGCCCCGCCTAGCTGGCCGCGAAATTTCG‑CGCCAGCGGTCAGCCCgtcgtcctg                                                                           >  2:3294103/1‑91 (MQ=255)
                          ccTTTTGCCGGTGGTCCACGTACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTgc                                                                       <  2:1624834/100‑1 (MQ=255)
                            ttttGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTgcgc                                                                     >  1:674534/1‑100 (MQ=255)
                               tGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGtt                                                                  >  2:4114009/1‑100 (MQ=255)
                               tGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGtt                                                                  <  1:714349/100‑1 (MQ=255)
                                   ggtggtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCt                                                              <  2:285671/100‑1 (MQ=255)
                                   ggtggtCCACATACTGCCGATGAGGGAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCt                                                              <  1:998351/100‑1 (MQ=255)
                                     tggtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTgg                                                            <  1:1165798/100‑1 (MQ=255)
                                       gtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTc                                                          <  2:547501/100‑1 (MQ=255)
                                        tCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCt                                                         <  2:3029980/100‑1 (MQ=255)
                                        tCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCt                                                         >  1:1283673/1‑100 (MQ=255)
                                         ccACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTgg                                                            >  2:437307/1‑96 (MQ=255)
                                         ccACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTg                                                        >  1:1338178/1‑100 (MQ=255)
                                         ccACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTg                                                        <  1:1266738/100‑1 (MQ=255)
                                         ccACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTg                                                        >  1:2367803/1‑100 (MQ=255)
                                          cacaTACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGc                                                       >  2:787102/1‑100 (MQ=255)
                                            caTACTGGCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCg                                                     <  2:200844/100‑1 (MQ=255)
                                            caTACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCg                                                     <  1:2268980/100‑1 (MQ=255)
                                            caTACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCg                                                     <  1:3408620/100‑1 (MQ=255)
                                              tACTGCCGATGAGGTGGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGAt                                                   >  2:3517402/1‑100 (MQ=255)
                                              tACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGAt                                                   >  1:4011457/1‑100 (MQ=255)
                                              tACTGCCGATGAGGTAGTTAATCGCCCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGa                                                    >  2:1245872/1‑99 (MQ=255)
                                                cTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTa                                                 <  1:1090811/100‑1 (MQ=255)
                                                 tGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTAc                                                >  2:923804/1‑100 (MQ=255)
                                                  gCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACg                                               >  2:764489/1‑100 (MQ=255)
                                                   ccGATGAGGTAGTTAATCGCCCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGc                                              <  2:710835/100‑1 (MQ=255)
                                                     gATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTg                                            >  2:2448209/1‑100 (MQ=255)
                                                       ggAGGTAGTTAAGCGCGCCGGGAGGCTGGCGGCGAAATTTCG‑CGCCAGCTGTCAGCCCGTGGCTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGc                                              >  2:1659055/2‑96 (MQ=255)
                                                         aGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGGCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGc                                        <  2:823054/100‑1 (MQ=255)
                                                         aGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGc                                        <  1:2448214/100‑1 (MQ=255)
                                                         aGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGc                                        <  2:2256344/100‑1 (MQ=255)
                                                          ggTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCa                                       <  1:1453929/100‑1 (MQ=255)
                                                            tAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCAtt                                     <  2:3940864/100‑1 (MQ=255)
                                                             aGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTa                                    >  1:786012/1‑100 (MQ=255)
                                                             aGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTa                                    >  2:3008392/1‑100 (MQ=255)
                                                             aGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTa                                    <  2:1193749/100‑1 (MQ=255)
                                                              gTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTATGCTGAAGCATTaa                                   <  1:2842052/100‑1 (MQ=255)
                                                                  aTCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACa                               <  1:102068/100‑1 (MQ=255)
                                                                   tCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAg                              <  2:568267/100‑1 (MQ=255)
                                                                       gcCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCgg                          <  2:1058288/100‑1 (MQ=255)
                                                                       gcCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCgg                          <  2:1060165/100‑1 (MQ=255)
                                                                         ctgaccgCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTc                        >  1:3010217/7‑100 (MQ=38)
                                                                          gggAATCTGGCGGCGAAATTTCG‑TGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCg                       <  2:3831597/100‑1 (MQ=255)
                                                                          gggAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCg                       <  2:3435160/100‑1 (MQ=255)
                                                                           ggAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCg                       >  2:1013845/1‑99 (MQ=255)
                                                                           gaccgCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGa                      >  2:2433379/5‑100 (MQ=38)
                                                                            gAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATGACGCTGAAGCATTAAAACAGCCGGTCGAt                     >  2:2523413/1‑100 (MQ=255)
                                                                            accgCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGAt                     >  2:4148691/4‑100 (MQ=38)
                                                                            accgCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGAt                     >  2:273705/4‑100 (MQ=38)
                                                                            accgCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGAt                     <  2:1684699/97‑1 (MQ=38)
                                                                             ccgCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg                    >  1:2341857/3‑100 (MQ=39)
                                                                             ccgCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg                    >  2:3276916/3‑100 (MQ=39)
                                                                             ccgCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg                    >  1:3855001/3‑100 (MQ=39)
                                                                             aaGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg                    <  2:190163/100‑1 (MQ=255)
                                                                              cgCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgc                   >  1:217539/2‑100 (MQ=39)
                                                                              cgCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCAGACGATgc                   >  2:771966/2‑100 (MQ=38)
                                                                               gCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgcg                  >  1:354541/1‑100 (MQ=255)
                                                                               gCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgcg                  >  1:408297/1‑100 (MQ=255)
                                                                               gCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgcg                  >  2:2729198/1‑100 (MQ=255)
                                                                               gCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgcg                  >  2:2686175/1‑100 (MQ=255)
                                                                               gCTGGC‑GCGAAATTTCGCCGCCAGCGGTCAGCCCCTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGACGATgcg                  >  2:2500569/1‑100 (MQ=255)
                                                                                 tGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGcc                >  1:246145/1‑100 (MQ=255)
                                                                                 tGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGcc                <  1:565211/100‑1 (MQ=255)
                                                                                  ggcggcGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGccc               <  2:1338177/100‑1 (MQ=255)
                                                                                    cggcgAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTc             >  1:179599/1‑100 (MQ=255)
                                                                                     ggcgAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCa            >  1:3629483/1‑100 (MQ=255)
                                                                                       cgAAATTTCG‑CGGCAGCGGTCAGCCCGTGTTGCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCAcc          <  2:173883/100‑1 (MQ=255)
                                                                                         aaaTTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCgg        <  1:3832412/100‑1 (MQ=255)
                                                                                          aaTTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGc       >  2:1505654/1‑100 (MQ=255)
                                                                                            tttCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGc       <  2:1300553/98‑1 (MQ=255)
                                                                                               cg‑cgCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAg  <  2:623456/100‑1 (MQ=255)
                                                                                               cg‑cgCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAg  <  1:2369205/100‑1 (MQ=255)
                                                                                                 |                                                                                                  
TGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCG‑CGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAG  >  NC_000913/1950869‑1951063

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: ‑