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.

TGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAG  >  NC_000913/1950869‑1951063
                                                                                                 |                                                                                                  
tgATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGCCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCCGCGAGATATCgcgc                                                                                                 >  2:3667567/1‑100 (MQ=255)
 gATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGcc                                                                                                >  1:545161/1‑100 (MQ=255)
 gATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGcc                                                                                                <  1:3349048/100‑1 (MQ=255)
    ngATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGc                                                                                             >  1:1695236/2‑100 (MQ=255)
    cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTCGCCGCCAGc                                                                                             <  2:283856/100‑1 (MQ=255)
    cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTCGCCGCCAGc                                                                                             <  1:496529/100‑1 (MQ=255)
    cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTCGCCGCCAGc                                                                                             <  1:446794/100‑1 (MQ=255)
    cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTCGCCGCCAGc                                                                                             <  1:3585111/100‑1 (MQ=255)
    cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTCGCCGCCAGc                                                                                             <  2:735125/100‑1 (MQ=255)
    cGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTCGCCGCCAGc                                                                                             <  1:2090255/100‑1 (MQ=255)
     gATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCg                                                                                            >  2:1233594/1‑100 (MQ=255)
        ttACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTc                                                                                         <  2:3722358/100‑1 (MQ=255)
         tACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCa                                                                                        >  1:3835369/1‑100 (MQ=255)
           cAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGTGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGc                                                                                      >  1:190163/1‑100 (MQ=255)
           cAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGc                                                                                      >  1:1877481/1‑100 (MQ=255)
           cAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGc                                                                                      >  1:4189761/1‑100 (MQ=255)
           cAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGc                                                                                      >  1:1019702/1‑100 (MQ=255)
           cAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGCGAAATTTCGCCGCCAGcttcccgg                                                                                      <  1:1590239/100‑8 (MQ=37)
           cAAGAAGGCAGCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGc                                                                                      <  1:178285/100‑1 (MQ=255)
            aagaagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGcc                                                                                     >  2:3406579/1‑100 (MQ=255)
            aagaagGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGcc                                                                                     <  1:1370575/100‑1 (MQ=255)
             agaagGCATCTTACCTTTTGGCGGTGGTCCACATACTGCCGATTAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGcc                                                                                     <  2:3963047/99‑1 (MQ=255)
              taacgCATCGTACCTTTTGCCGGTTGTCCACATACTTCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGccct                                                                                   <  2:4034852/96‑2 (MQ=255)
                 ggCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCgtgt                                                                                <  1:2206023/100‑1 (MQ=255)
                 ggCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCgtgt                                                                                <  2:2786764/100‑1 (MQ=255)
                  gCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGtt                                                                               <  1:318752/100‑1 (MQ=255)
                    aTCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTc                                                                             >  2:2252810/1‑100 (MQ=255)
                    aTCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTc                                                                             <  2:3196221/100‑1 (MQ=255)
                        tACCTTTTGCCGGTGGTCCACATACCGCCATCGAGGTAGTTAATCGCCCCGCCTAGCTGGCCGCGAAATTTCGCGCCAGCGGTCAGCCCgtcgtcctg                                                                           >  2:3294103/1‑91 (MQ=255)
                          ccTTTTGCCGGTGGTCCACGTACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTgc                                                                       <  2:1624834/100‑1 (MQ=255)
                            ttttGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTgcgc                                                                     >  1:674534/1‑100 (MQ=255)
                               tGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGtt                                                                  >  2:4114009/1‑100 (MQ=255)
                               tGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGtt                                                                  <  1:714349/100‑1 (MQ=255)
                                   ggtggtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCt                                                              <  2:285671/100‑1 (MQ=255)
                                   ggtggtCCACATACTGCCGATGAGGGAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCt                                                              <  1:998351/100‑1 (MQ=255)
                                     tggtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTgg                                                            <  1:1165798/100‑1 (MQ=255)
                                       gtCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTc                                                          <  2:547501/100‑1 (MQ=255)
                                        tCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCt                                                         <  2:3029980/100‑1 (MQ=255)
                                        tCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCt                                                         >  1:1283673/1‑100 (MQ=255)
                                         ccACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTgg                                                            >  2:437307/1‑96 (MQ=255)
                                         ccACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTg                                                        >  1:1338178/1‑100 (MQ=255)
                                         ccACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTg                                                        <  1:1266738/100‑1 (MQ=255)
                                         ccACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTg                                                        >  1:2367803/1‑100 (MQ=255)
                                          cacaTACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGc                                                       >  2:787102/1‑100 (MQ=255)
                                            caTACTGGCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCg                                                     <  2:200844/100‑1 (MQ=255)
                                            caTACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCg                                                     <  1:2268980/100‑1 (MQ=255)
                                            caTACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCg                                                     <  1:3408620/100‑1 (MQ=255)
                                              tACTGCCGATGAGGTGGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGAt                                                   >  2:3517402/1‑100 (MQ=255)
                                              tACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGAt                                                   >  1:4011457/1‑100 (MQ=255)
                                              tACTGCCGATGAGGTAGTTAATCGCCCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGa                                                    >  2:1245872/1‑99 (MQ=255)
                                                cTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTa                                                 <  1:1090811/100‑1 (MQ=255)
                                                 tGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTAc                                                >  2:923804/1‑100 (MQ=255)
                                                  gCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACg                                               >  2:764489/1‑100 (MQ=255)
                                                   ccGATGAGGTAGTTAATCGCCCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGc                                              <  2:710835/100‑1 (MQ=255)
                                                     gATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTg                                            >  2:2448209/1‑100 (MQ=255)
                                                       ggAGGTAGTTAAGCGCGCCGGGAGGCTGGCGGCGAAATTTCGCGCCAGCTGTCAGCCCGTGGCTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGc                                              >  2:1659055/2‑96 (MQ=255)
                                                         aGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGGCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGc                                        <  2:823054/100‑1 (MQ=255)
                                                         aGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGc                                        <  1:2448214/100‑1 (MQ=255)
                                                         aGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGc                                        <  2:2256344/100‑1 (MQ=255)
                                                          ggTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCa                                       <  1:1453929/100‑1 (MQ=255)
                                                            tAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCAtt                                     <  2:3940864/100‑1 (MQ=255)
                                                             aGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTa                                    >  1:786012/1‑100 (MQ=255)
                                                             aGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTa                                    >  2:3008392/1‑100 (MQ=255)
                                                             aGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTa                                    <  2:1193749/100‑1 (MQ=255)
                                                              gTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTATGCTGAAGCATTaa                                   <  1:2842052/100‑1 (MQ=255)
                                                                  aTCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACa                               <  1:102068/100‑1 (MQ=255)
                                                                   tCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAg                              <  2:568267/100‑1 (MQ=255)
                                                                       gcCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCgg                          <  2:1058288/100‑1 (MQ=255)
                                                                       gcCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCgg                          <  2:1060165/100‑1 (MQ=255)
                                                                         ctgaccgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTc                        >  1:3010217/7‑100 (MQ=38)
                                                                          gggAATCTGGCGGCGAAATTTCGTGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCg                       <  2:3831597/100‑1 (MQ=255)
                                                                          gggAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCg                       <  2:3435160/100‑1 (MQ=255)
                                                                           ggAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCg                       >  2:1013845/1‑99 (MQ=255)
                                                                           gaccgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGa                      >  2:2433379/5‑100 (MQ=38)
                                                                            gAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATGACGCTGAAGCATTAAAACAGCCGGTCGAt                     >  2:2523413/1‑100 (MQ=255)
                                                                            accgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGAt                     >  2:4148691/4‑100 (MQ=38)
                                                                            accgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGAt                     >  2:273705/4‑100 (MQ=38)
                                                                            accgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGAt                     <  2:1684699/97‑1 (MQ=38)
                                                                             ccgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg                    >  1:2341857/3‑100 (MQ=39)
                                                                             ccgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg                    >  2:3276916/3‑100 (MQ=39)
                                                                             ccgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg                    >  1:3855001/3‑100 (MQ=39)
                                                                             aaGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATg                    <  2:190163/100‑1 (MQ=255)
                                                                              cgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgc                   >  1:217539/2‑100 (MQ=39)
                                                                              cgCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCAGACGATgc                   >  2:771966/2‑100 (MQ=38)
                                                                               gCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgcg                  >  1:354541/1‑100 (MQ=255)
                                                                               gCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgcg                  >  1:408297/1‑100 (MQ=255)
                                                                               gCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgcg                  >  2:2729198/1‑100 (MQ=255)
                                                                               gCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATgcg                  >  2:2686175/1‑100 (MQ=255)
                                                                               gCTGGCGCGAAATTTCGCCGCCAGCGGTCAGCCCCTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGACGATgcg                  >  2:2500569/1‑100 (MQ=255)
                                                                                 tGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGcc                >  1:246145/1‑100 (MQ=255)
                                                                                 tGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGcc                <  1:565211/100‑1 (MQ=255)
                                                                                  ggcggcGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGccc               <  2:1338177/100‑1 (MQ=255)
                                                                                    cggcgAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTc             >  1:179599/1‑100 (MQ=255)
                                                                                     ggcgAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCa            >  1:3629483/1‑100 (MQ=255)
                                                                                       cgAAATTTCGCGGCAGCGGTCAGCCCGTGTTGCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCAcc          <  2:173883/100‑1 (MQ=255)
                                                                                         aaaTTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCgg        <  1:3832412/100‑1 (MQ=255)
                                                                                          aaTTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGc       >  2:1505654/1‑100 (MQ=255)
                                                                                            tttCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGc       <  2:1300553/98‑1 (MQ=255)
                                                                                               cgcgCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAg  <  2:623456/100‑1 (MQ=255)
                                                                                               cgcgCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAg  <  1:2369205/100‑1 (MQ=255)
                                                                                                 |                                                                                                  
TGATCGATTTACAAGAAGGCATCTTACCTTTTGCCGGTGGTCCACATACTGCCGATGAGGTAGTTAATCGCGCCGGGAAGCTGGCGGCGAAATTTCGCGCCAGCGGTCAGCCCGTGTTTCTGGTGCGCGTTGGCTGGTCTGCCGATTACGCTGAAGCATTAAAACAGCCGGTCGATGCGCCCTCACCGGCTAAAG  >  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: