Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,410,289 T→A 5.5% I320N (ATT→AAT)  nadA → quinolinate synthase iron‑sulfur cluster subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,410,2890TA 5.5% 227.2 / 2.5 73I320N (ATT→AAT) nadAquinolinate synthase iron‑sulfur cluster subunit
Reads supporting (aligned to +/- strand):  ref base T (29/40);  new base A (2/2);  total (31/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

GCAGCAGTTGTGCCCGGACAAGGAATTCGTCGAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGGCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGACGCCAGACCGCGTCATACTCTTCGCGGGCTCGCCCACTCCCACAGAT  >  NC_002947/1410143‑1410436
                                                                                                                                                  |                                                                                                                                                   
gcagcaGTTGTGCCCGGACAAGGAATTCGTCGAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAAttt                                                                                                                                                   >  2:50224/1‑149 (MQ=255)
  agcagTTGTGCCCGGACAAGGAATTCGTCGAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTc                                                                                                                                                 <  3:242495/149‑1 (MQ=255)
  agcagTTGTGCCCGGACAAGGAATTCGTCGAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTc                                                                                                                                                 <  2:314608/149‑1 (MQ=255)
  agcagTTGTGCCCGGACAAGGAATTCGTCGAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTc                                                                                                                                                 <  1:89080/149‑1 (MQ=255)
       ttGTGCCCGGACAAGGAATTCGTCGAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGa                                                                                                                                            >  4:95850/1‑149 (MQ=255)
         gtgCCCGGACAAGGAATTCGTCGAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAAATTTCGTTggt                                                                                                                                           <  2:256635/148‑3 (MQ=255)
            cccGGACAAGGAATTCGTCGAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCgg                                                                                                                                       >  8:237362/1‑149 (MQ=255)
                   aaGGAATTCGTCGAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGt                                                                                                                                <  3:177017/149‑1 (MQ=255)
                     ggAATTCGTCGAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGc                                                                                                                              <  4:303938/149‑1 (MQ=255)
                         ttCGTCGAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCaa                                                                                                                          <  5:237924/149‑1 (MQ=255)
                           cgtcgAAGCCCCCACCGCCGGCAACGGCGCGGCATGCTGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAg                                                                                                                         <  6:180850/148‑1 (MQ=255)
                            gtcgAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAgg                                                                                                                        >  5:160530/1‑148 (MQ=255)
                            gtcgAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAgg                                                                                                                        <  7:237362/148‑1 (MQ=255)
                            gtcgAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGc                                                                                                                       >  5:201755/1‑149 (MQ=255)
                            gtcgAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGc                                                                                                                       >  1:89364/1‑149 (MQ=255)
                            gtcgAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGc                                                                                                                       >  7:208366/1‑149 (MQ=255)
                              cgAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGa                                                                                                                     <  4:237279/149‑1 (MQ=255)
                               gAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGAt                                                                                                                    <  3:26894/149‑1 (MQ=255)
                               gAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGAt                                                                                                                    <  6:154868/149‑1 (MQ=255)
                                  gCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATc                                                                                                                   >  3:21722/1‑147 (MQ=255)
                                    ccccACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGc                                                                                                               >  5:329442/1‑149 (MQ=255)
                                    ccccACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGc                                                                                                               <  3:95850/149‑1 (MQ=255)
                                     cccACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGcc                                                                                                              <  2:36480/149‑1 (MQ=255)
                                        aCCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCg                                                                                                                                                >  2:272364/1‑112 (MQ=255)
                                        aCCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCg                                                                                                                                                <  1:272364/112‑1 (MQ=255)
                                             cgGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATc                                                                                                      >  4:211447/1‑149 (MQ=255)
                                             cgGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGAGAAAGCCGCTTaaac                                                                                                      >  4:192760/1‑147 (MQ=255)
                                                     gcgcGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGcc                                                                                                              <  6:214880/133‑1 (MQ=255)
                                                     gcgcGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGcc                                                                                                              >  5:214880/1‑133 (MQ=255)
                                                     gcgcGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCt                                                                                                           >  4:321935/1‑136 (MQ=255)
                                                     gcgcGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTg                                                                                              >  3:334329/1‑149 (MQ=255)
                                                     gcgcGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTg                                                                                              >  1:237876/1‑149 (MQ=255)
                                                         ggCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACt                                                                                          <  6:36392/149‑1 (MQ=255)
                                                          gccccccGCACCGGGGCGCATGGCCCGGGGAACGCGAGGAACACCCCGGCGCGGGGGCTGGAATGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGCACtt                                                                                         <  4:210163/144‑1 (MQ=255)
                                                               ccGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTAccc                                                                                    <  3:211447/149‑1 (MQ=255)
                                                               ccGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTAcc                                                                                     <  3:2220/148‑1 (MQ=255)
                                                               ccGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAAATTTCGTTggtgc                                                                                                                                         >  5:22579/1‑92 (MQ=255)
                                                               ccGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAAATTTCGTTggtgc                                                                                                                                         <  6:22579/96‑5 (MQ=255)
                                                                   aGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAggc                                                                                <  2:89364/149‑1 (MQ=255)
                                                                        tgCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGgcggcgc                                                                           <  6:201755/149‑1 (MQ=255)
                                                                           gcATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGgcggcgcggc                                                                        >  8:23400/1‑149 (MQ=255)
                                                                           gcATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGgcggcgcggc                                                                        <  2:237876/149‑1 (MQ=255)
                                                                             aTTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGGCGGCGCGGCTg                                                                      <  6:317946/149‑1 (MQ=255)
                                                                                gCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGGCGGCGCGGCTGaa                                                                    <  1:50224/148‑1 (MQ=255)
                                                                                 cccGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGGCGGCGCGGCTGAAGc                                                                  >  2:102750/1‑149 (MQ=255)
                                                                                            cGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTCCCCCGGCGGCGCGGCTGAAGCTGTCTGGTAAt                                                       >  4:183145/1‑149 (MQ=255)
                                                                                               tGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGc                                                                                                            >  2:251989/1‑93 (MQ=255)
                                                                                               tGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGc                                                                                                            <  1:251989/93‑1 (MQ=255)
                                                                                                gAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGGCGGCGCGGCTGAAGCTGTCTGGTAATGCCt                                                   >  3:55554/1‑149 (MQ=255)
                                                                                                gAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGGCGGCGCGGCTGAAGCTGTCTGGTAATGCCt                                                   >  2:80153/1‑149 (MQ=255)
                                                                                                 aaCACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGGCGGCGCGGCTGAAGCTGTCTGGTAATGCCt                                                   <  8:204518/148‑1 (MQ=255)
                                                                                                 aaCACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGGCGGCGCGGCTGAAGCTGTCTGGTAATGCCt                                                   >  4:205158/1‑148 (MQ=255)
                                                                                                 aaCACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGGCGGCGCGGCTGAAGCTGTCTGGTAATGCCTg                                                  <  4:21722/149‑1 (MQ=255)
                                                                                                   cacCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGGCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGAc                                                <  8:208366/149‑1 (MQ=255)
                                                                                                    acCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGGCGGCGCGGCTGAAg                                                                   <  7:291043/129‑1 (MQ=255)
                                                                                                    acCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGGCGGCGCGGCTGAAg                                                                   >  8:291043/1‑129 (MQ=255)
                                                                                                      ccTGGAGCGGTTGATGGATTGCATGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGATGGTGCCCAAGTCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGGCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGACGc                                              >  8:34296/1‑148 (MQ=255)
                                                                                                      ccTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGGCGGCGCGGCTGAAGCt                                                                 >  7:200583/1‑129 (MQ=255)
                                                                                                      ccTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGGCGGCGCGGCTGAAGCt                                                                 <  8:200583/129‑1 (MQ=255)
                                                                                                       cTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGGCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGACGCCa                                            <  7:23400/149‑1 (MQ=255)
                                                                                                          gAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGGCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGACGCCAGAc                                         <  7:7768/149‑1 (MQ=255)
                                                                                                                      gagtGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGGCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGACGCCAGACCGCGTCATActc                             <  3:192760/146‑1 (MQ=255)
                                                                                                                           cccgCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGGCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGACGCCAGACCGCGTCATACTCTTcgc                        >  1:141250/4‑149 (MQ=255)
                                                                                                                           ccTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGGCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGACGCCAGACCGCGTCATACTCTTcgc                        >  7:264945/1‑149 (MQ=255)
                                                                                                                            cTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGgcggcgc                                                                           >  5:221759/1‑97 (MQ=255)
                                                                                                                            cTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGgcggcgc                                                                           <  6:221759/97‑1 (MQ=255)
                                                                                                                            cTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGGCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGACGCCAGACCGCGTCATACTCTTcgcg                       <  4:334329/149‑1 (MQ=255)
                                                                                                                             tGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGGCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGACGCCAGACCGCGTCATACTCTTCGCgg                      >  6:120140/1‑149 (MQ=255)
                                                                                                                               cGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGGCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGACGCCAGACCGCGTCATACTCTTCGCgg                      <  5:120140/147‑1 (MQ=255)
                                                                                                                                   aCGGCACCAACTAAATTTTCGTTGACCCGGCGGTGCTGCCCAAGCCGATCCAGCCGCTTAACCGCATGCTGGACTTTACCCATGCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGACGCCAGACCGCGTCATACTCTTCGCGGGCTCGc                <  1:102750/149‑1 (MQ=255)
                                                                                                                                       cACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGGCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGACGCCAGACCGCGTCATACTCTTCGCGGGCTCGCCCAc            <  4:55554/149‑1 (MQ=255)
                                                                                                                                        gtcAACGAAAATTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGGCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGACGCCAGACCGCGTCATACTCTTCGCGGGCTCGCCCACt           >  7:203764/3‑149 (MQ=255)
                                                                                                                                              taaaaTTTCGTTGACCCGGCGCTGGTGCACAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGGCGGCGCGGCAGAAGCTGTCTGGTAATGc                                                     >  2:293057/2‑101 (MQ=255)
                                                                                                                                              gtaaTTTTCGTTGACCCGTCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATGGCATGCTGGACTTTACCCAGGCGGCGCGGCTTAAGCTGTCTGGTAATGc                                                     <  1:293057/99‑1 (MQ=255)
                                                                                                                                                 aTTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGGCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGACGCCAGACCGCGTCATACTCTTCGCGGGCTCGCCCACTCCCACAGat  <  6:329442/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCAGCAGTTGTGCCCGGACAAGGAATTCGTCGAAGCCCCCACCGCCGGCAACGGCGCGGCATGCCGCAGCTGTGCGCATTGCCCGTGGATGGCGATGAACACCCTGGAGCGGGTGCTGGATTGCCTGCGTAACGGCACCAACGAAATTTTCGTTGACCCGGCGCTGGTGCCCAAGGCGATCAAGCCGCTTAATCGCATGCTGGACTTTACCCAGGCGGCGCGGCTGAAGCTGTCTGGTAATGCCTGACGCCAGACCGCGTCATACTCTTCGCGGGCTCGCCCACTCCCACAGAT  >  NC_002947/1410143‑1410436

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: