Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,852,038 A→C 6.2% Q51P (CAG→CCG) ‡ pncA → nicotinamidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,852,0380AC 6.2% 295.8 / 8.5 97Q51P (CAG→CCG) ‡pncAnicotinamidase
Reads supporting (aligned to +/- strand):  ref base A (45/46);  new base C (3/3);  total (48/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

GCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACATCCGTTACTGC  >  NC_000913/1851889‑1852187
                                                                                                                                                     |                                                                                                                                                     
gCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCGCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCa                                                                                                                                                       <  2:877041/150‑1 (MQ=255)
  ccccTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAgg                                                                                                                                                     >  2:540553/1‑150 (MQ=255)
  ccccTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAgg                                                                                                                                                     <  1:809790/150‑1 (MQ=255)
   cccTCGCGACCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGACGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGa                                                                                                                                                    >  2:564347/1‑150 (MQ=255)
    ccTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGAc                                                                                                                                                   <  1:1518891/150‑1 (MQ=255)
         cgcCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGTATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCCTGACTGGcg                                                                                                                                              <  2:723035/150‑2 (MQ=255)
            ccTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCCTGACTGGcgata                                                                                                                                           <  2:439306/150‑5 (MQ=255)
              tGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCgg                                                                                                                                         <  2:801662/150‑1 (MQ=255)
              tGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCgg                                                                                                                                         >  2:1208305/1‑150 (MQ=255)
              tGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGACTGATTGAGTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCCTGACTGGcgataac                                                                                                                                         <  2:93413/150‑7 (MQ=255)
               gTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGc                                                                                                                                        <  2:240165/150‑1 (MQ=255)
                  aCTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGaa                                                                                                                                     <  1:1297693/150‑1 (MQ=255)
                  aCTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGaa                                                                                                                                     >  1:813588/1‑150 (MQ=255)
                     ggTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCa                                                                                                                                  <  1:724605/150‑1 (MQ=255)
                     ggTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCa                                                                                                                                  >  2:641826/1‑150 (MQ=255)
                     ggTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCa                                                                                                                                  >  2:481532/1‑150 (MQ=255)
                     ggTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCa                                                                                                                                  <  1:194200/150‑1 (MQ=255)
                      gTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCAc                                                                                                                                 <  1:834040/150‑1 (MQ=255)
                        cGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACgg                                                                                                                               <  1:959420/150‑1 (MQ=255)
                        cGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACgg                                                                                                                               >  2:1066818/1‑150 (MQ=255)
                           tttACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAg                                                                                                                            <  1:540553/150‑1 (MQ=255)
                               cAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGtttt                                                                                                                        >  1:1017556/1‑150 (MQ=255)
                                 aaaTGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGc                                                                                                                      <  1:1066818/150‑1 (MQ=255)
                                  aaTGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGcc                                                                                                                     <  2:1338264/150‑1 (MQ=255)
                                    tGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAg                                                                                                                   >  1:706085/1‑150 (MQ=255)
                                     gATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGt                                                                                                                  >  2:995522/1‑150 (MQ=255)
                                       tttCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCa                                                                                                                <  1:1249361/150‑1 (MQ=255)
                                         tCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGc                                                                                                              >  1:672420/1‑150 (MQ=255)
                                          cTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCa                                                                                                             <  1:977879/150‑1 (MQ=255)
                                          cTGTACTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCa                                                                                                             >  2:769761/1‑150 (MQ=255)
                                           tgtgCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCAc                                                                                                            >  1:1113142/1‑150 (MQ=255)
                                                tggtggCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGgtgt                                                                                                       <  2:690018/150‑1 (MQ=255)
                                                      cgcgcTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCa                                                                                                                >  2:186150/1‑135 (MQ=255)
                                                      cgcgcTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCa                                                                                                                <  1:186150/135‑1 (MQ=255)
                                                       gcgcTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCt                                                                                                <  1:329491/150‑1 (MQ=255)
                                                          cTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTtat                                                                                             <  1:1273363/150‑1 (MQ=255)
                                                               cGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGcc                                                                                        <  2:1113142/150‑1 (MQ=255)
                                                                     ggAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCa                                                                                  >  2:580853/1‑150 (MQ=255)
                                                                     ggAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCa                                                                                  >  2:153121/1‑150 (MQ=255)
                                                                     ggAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCa                                                                                  >  2:321749/1‑150 (MQ=255)
                                                                     ggAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGAATGGCACCCGGCGAAGCCCGGCAGTTTTGCCAGTCAGCACGGGGTAGAGCCTTATACGCCAGGCCa                                                                                  >  2:466484/1‑150 (MQ=255)
                                                                           tGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGa                                                                            >  2:996664/1‑150 (MQ=255)
                                                                           tGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGa                                                                            >  1:1309353/1‑150 (MQ=255)
                                                                                gTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGtt                                                                       >  1:1056255/1‑150 (MQ=255)
                                                                                  aCGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTg                                                                     >  2:901092/1‑150 (MQ=255)
                                                                                  aCGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTg                                                                     >  1:672987/1‑150 (MQ=255)
                                                                                   cGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGc                                                                    <  1:481532/150‑1 (MQ=255)
                                                                                    ggTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGcc                                                                   <  2:813588/150‑1 (MQ=255)
                                                                                       ggATGTCGCTAACCGCCTGATTGACTGGTGCCAGTGGCGCGGTGAAGCGGGTATCGCCAGTCAGGAGTGGCACCCGGCGACTCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATCCGACAGGCCAACTCTACGGTTTGCcaga                                                                <  2:1397649/150‑3 (MQ=255)
                                                                                       ggATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCcaca                                                                <  2:846771/150‑1 (MQ=255)
                                                                                         aTGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACaaa                                                              >  2:466272/1‑150 (MQ=255)
                                                                                          tGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAAc                                                             >  1:1124197/1‑150 (MQ=255)
                                                                                          tGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAAc                                                             >  1:1404730/1‑150 (MQ=255)
                                                                                           gTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAAcc                                                            >  1:569029/1‑150 (MQ=255)
                                                                                           gTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAAcc                                                            >  2:783012/1‑150 (MQ=255)
                                                                                            tCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCt                                                           >  2:354079/1‑150 (MQ=255)
                                                                                             cGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCtt                                                          >  2:1183331/1‑150 (MQ=255)
                                                                                                 aaCCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGCGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAACCCTTCTgg                                                      >  1:685010/1‑150 (MQ=255)
                                                                                                 aaCCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTgg                                                      >  1:349203/1‑150 (MQ=255)
                                                                                                 aaCCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTgg                                                      <  2:586003/150‑1 (MQ=255)
                                                                                                   ccGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGcc                                                    <  2:1156336/150‑1 (MQ=255)
                                                                                                   ccGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGcc                                                    <  2:1309353/150‑1 (MQ=255)
                                                                                                   ccGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGACAACTCGACGGTTTGCCACAAACCTTCTGGcc                                                    >  2:970125/1‑150 (MQ=255)
                                                                                                    cGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCa                                                   <  2:1124197/150‑1 (MQ=255)
                                                                                                    cGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCa                                                   >  2:524211/1‑150 (MQ=255)
                                                                                                     gCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGGACCCGGCGAATCACGGCAGTTTTGCCAGGCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAg                                                  <  2:1141225/150‑1 (MQ=255)
                                                                                                     gCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAg                                                  <  1:564347/150‑1 (MQ=255)
                                                                                                        tGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATc                                               <  2:1099151/150‑1 (MQ=255)
                                                                                                            tGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACtg                                           <  2:777559/150‑1 (MQ=255)
                                                                                                            tGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACtg                                           <  2:196592/150‑1 (MQ=255)
                                                                                                            tGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACtg                                           <  2:381093/150‑1 (MQ=255)
                                                                                                              aCTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACtgtg                                         >  1:713158/1‑150 (MQ=255)
                                                                                                                gggggCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGGCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACtgtgtg                                       <  1:1155059/146‑1 (MQ=255)
                                                                                                                 ggTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACTGTTTGCCACAAACCTACTGGCCACATCACTGTGTGc                                      >  2:682962/1‑150 (MQ=255)
                                                                                                                 ggTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGc                                      <  1:321749/150‑1 (MQ=255)
                                                                                                                  gTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACTGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCa                                     >  2:840422/1‑150 (MQ=255)
                                                                                                                  gTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCa                                     <  1:354079/150‑1 (MQ=255)
                                                                                                                      cAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAAc                                 <  2:1014233/150‑1 (MQ=255)
                                                                                                                        gTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAg                               <  1:1179441/150‑1 (MQ=255)
                                                                                                                          cgcgcgGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATATGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTg                             <  1:641826/150‑1 (MQ=255)
                                                                                                                          cgcgcgGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTg                             >  2:968162/1‑150 (MQ=255)
                                                                                                                              cgGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAgg                         >  2:1218013/1‑150 (MQ=255)
                                                                                                                                gTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGgcg                       <  1:1564113/150‑1 (MQ=255)
                                                                                                                                       ggTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAAtt                >  1:1204718/1‑150 (MQ=255)
                                                                                                                                       gccgggtgCCAGTCCTGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAAtt                >  2:709387/8‑150 (MQ=255)
                                                                                                                                        gTTATCTCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATTCGCGAGGCCAACTCGACGGTTTGCCTCAATGCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTa               <  1:231550/150‑1 (MQ=255)
                                                                                                                                        ccgggtgCCAGTCCTGACCGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTa               >  1:558506/7‑150 (MQ=255)
                                                                                                                                          tATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACa             <  1:901092/150‑1 (MQ=255)
                                                                                                                                          tATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACa             <  2:672987/150‑1 (MQ=255)
                                                                                                                                           ggtgCCAGTCCTGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACAt            >  2:801264/4‑150 (MQ=255)
                                                                                                                                            tCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACATc           >  1:929542/1‑150 (MQ=255)
                                                                                                                                             cGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACATcc          <  2:233796/150‑1 (MQ=255)
                                                                                                                                                cAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCGAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACATCCGtt       <  2:259668/150‑1 (MQ=255)
                                                                                                                                                cAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACATCCGtt       <  2:1017556/150‑1 (MQ=255)
                                                                                                                                                 aGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACATCCGTTa      >  2:1546228/1‑150 (MQ=255)
                                                                                                                                                 aGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACATCCGTTa      >  1:1157975/1‑150 (MQ=255)
                                                                                                                                                    cAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCGGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACATCCGTTACTg   >  2:106516/1‑150 (MQ=255)
                                                                                                                                                    cAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACATCCGTTACTg   >  1:680679/1‑150 (MQ=255)
                                                                                                                                                     aGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACATCCGTTACTGc  <  1:682962/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCCCCTCGCGCCCTGTTACTGGTCGATTTACAAAATGATTTCTGTGCTGGTGGCGCGCTCGCCGTGCCGGAAGGTGACAGTACGGTGGATGTCGCTAACCGCCTGATTGACTGGTGCCAGTCGCGCGGTGAAGCGGTTATCGCCAGTCAGGACTGGCACCCGGCGAATCACGGCAGTTTTGCCAGTCAGCACGGTGTAGAGCCTTATACGCCAGGCCAACTCGACGGTTTGCCACAAACCTTCTGGCCAGATCACTGTGTGCAGAACAGTGAAGGCGCACAATTACATCCGTTACTGC  >  NC_000913/1851889‑1852187

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