Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 782,800 C→T 5.9% A239V (GCG→GTG)  nadA → quinolinate synthase, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913782,8000CT 5.9% 347.8 / 7.9 102A239V (GCG→GTG) nadAquinolinate synthase, subunit A
Reads supporting (aligned to +/- strand):  ref base C (54/42);  new base T (3/3);  total (57/45)
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.70e-01

Only 100 of 108 total aligned reads displayed.

GGGTGGAGACATTCTATGCTGGCAGGGTGCCTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGT  >  NC_000913/782651‑782948
                                                                                                                                                     |                                                                                                                                                    
gggTGGAGACATTCTATGCTGGCAGGGTGCCTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGACTCACCACAAGCTATTGTCGATATGGCGGATGc                                                                                                                                                      >  1:1471600/1‑150 (MQ=255)
 ggTGGAGACATTCTATGCTGGCAGGGTGCCTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGcg                                                                                                                                                     <  1:1128323/150‑1 (MQ=255)
   tGGAGACATTCTATGCTGGCAGGGTGCCTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGcggt                                                                                                                                                   <  2:814417/150‑1 (MQ=255)
    ggAGACATTCTATGCTGGCAGGGTGCCTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGcggtc                                                                                                                                                  <  2:127052/150‑1 (MQ=255)
     gagaCATTCTATGCTGGCAGGGTGCCTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGcggtcg                                                                                                                                                 <  1:1016154/150‑1 (MQ=255)
     gagaCATTCTATGCTGGCAGGGTGCCTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGcggtcg                                                                                                                                                 <  1:983294/150‑1 (MQ=255)
        aCATTCTATGCTGGCAGGGTGCCTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGtt                                                                                                                                              <  1:333757/150‑1 (MQ=255)
          aTTCTATGCTGGCAGGGTGCCTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTcc                                                                                                                                            >  2:400138/1‑150 (MQ=255)
           ttCTATGCTGGCAGGGTGCCTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTcca                                                                                                                                           >  2:803583/1‑150 (MQ=255)
             cTATGCTGGCAGGGTGCCTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCGCCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTccacc                                                                                                                                         <  2:1129086/150‑1 (MQ=255)
              tATGCTGGCAGGGTGCCTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTccacca                                                                                                                                        >  1:771197/1‑150 (MQ=255)
              tATGCTGGCAGGGTGCCTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTccacca                                                                                                                                        >  1:1161594/1‑150 (MQ=255)
               aTGCTGGCAGGGTGCCTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAg                                                                                                                                       <  1:1498898/150‑1 (MQ=255)
                   tGGCAGGGTGCCTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCaa                                                                                                                                   >  1:1599797/1‑150 (MQ=255)
                   tGGCAGGGTGCCTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCaa                                                                                                                                   >  1:762919/1‑150 (MQ=255)
                     gCAGGGTGCCTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACt                                                                                                                                 >  2:1091823/1‑150 (MQ=255)
                     gCAGGGTGCCTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACt                                                                                                                                 >  2:1537219/1‑150 (MQ=255)
                     gCAGGGTGCCTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACt                                                                                                                                 >  2:1091825/1‑150 (MQ=255)
                     gCAGGGTGCCTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACt                                                                                                                                 >  1:129355/1‑150 (MQ=255)
                             ccTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCgctg                                                                                                                         <  2:191998/150‑1 (MQ=255)
                             ccTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCgctg                                                                                                                         <  2:192020/150‑1 (MQ=255)
                             ccTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCgctg                                                                                                                         <  2:454918/150‑1 (MQ=255)
                             ccTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCgctg                                                                                                                         <  1:1004447/150‑1 (MQ=255)
                              cTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCgctgc                                                                                                                        <  1:1267436/150‑1 (MQ=255)
                                  aTTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCg                                                                                                                    <  2:896149/150‑1 (MQ=255)
                                  aTTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCg                                                                                                                    <  2:148761/150‑1 (MQ=255)
                                  aTTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCg                                                                                                                    <  2:1676680/150‑1 (MQ=255)
                                  aTTGTGCAGGATGAATTTAAGACTCAGGCGTTAAACCGCTTGCAAGTAGTATACCCGGATGCTGACATCATGGTGCATACAGAATCACCACAAGCTATTGTCGATATGGGGTATGCGGTCGTTTCCACCAGTCAACTGATCGCTGCTTCg                                                                                                                    <  2:391023/150‑1 (MQ=255)
                                   ttGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGa                                                                                                                   <  1:927422/150‑1 (MQ=255)
                                     gtgCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGaaa                                                                                                                 >  2:148164/1‑150 (MQ=255)
                                     gtgCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGaaa                                                                                                                 >  2:1021172/1‑150 (MQ=255)
                                      tgCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGaaaa                                                                                                                >  2:1563038/1‑150 (MQ=255)
                                      tgCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGaaaa                                                                                                                >  2:1126612/1‑150 (MQ=255)
                                      tgCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGaaaa                                                                                                                >  1:679073/1‑150 (MQ=255)
                                      tgCATGATGAATCTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGaaaa                                                                                                                >  2:1563037/1‑150 (MQ=255)
                                       gCAGGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATCCTGCCATACTGGTGCATCCAGAATCACCACAAGCAATTGTCGAAAGGGCGGATGCGGTCGGTTCCACCAGACAACAGATCGCTGCTGCGAAAAc                                                                                                               >  1:232049/1‑150 (MQ=255)
                                        cATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACa                                                                                                              <  1:1297374/150‑1 (MQ=255)
                                         atgatgAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACAt                                                                                                             >  1:665855/1‑150 (MQ=255)
                                            atgaATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGc                                                                                                          >  1:617954/1‑150 (MQ=255)
                                             tgaATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGcc                                                                                                         >  2:401483/1‑150 (MQ=255)
                                                aTTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCcaca                                                                                                      <  2:312415/150‑1 (MQ=255)
                                                 tttAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAGCATTGCCACAt                                                                                                     >  2:597270/1‑150 (MQ=255)
                                                 tttAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACAt                                                                                                     >  2:414720/1‑150 (MQ=255)
                                                        cTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGc                                                                                              <  2:406333/150‑1 (MQ=255)
                                                         tCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCt                                                                                             <  1:803583/150‑1 (MQ=255)
                                                           aGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTa                                                                                           >  2:1011380/1‑150 (MQ=255)
                                                           aGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTa                                                                                           >  2:1580291/1‑150 (MQ=255)
                                                           aGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTa                                                                                           >  2:527379/1‑150 (MQ=255)
                                                           aGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTa                                                                                           >  2:578695/1‑150 (MQ=255)
                                                                  aaCCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGc                                                                                    >  1:845517/1‑150 (MQ=255)
                                                                  aaCCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGc                                                                                    >  1:845510/1‑150 (MQ=255)
                                                                  aaCCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGc                                                                                    <  2:472309/150‑1 (MQ=255)
                                                                   aCCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCa                                                                                   >  2:127945/1‑150 (MQ=255)
                                                                     ccGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAAc                                                                                 >  1:311291/1‑150 (MQ=255)
                                                                      cGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGTGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAAcc                                                                                >  1:621970/1‑150 (MQ=255)
                                                                      cGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGTGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAAcc                                                                                >  1:621965/1‑150 (MQ=255)
                                                                      cGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGTGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAAcc                                                                                >  1:621963/1‑150 (MQ=255)
                                                                            cAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCgg                                                                          <  2:1086243/150‑1 (MQ=255)
                                                                            cAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCgg                                                                          >  1:532003/1‑150 (MQ=255)
                                                                                        ccGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTAc                                                              >  2:81456/1‑150 (MQ=255)
                                                                                           gATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACaaa                                                           >  2:1528383/1‑150 (MQ=255)
                                                                                                      aCTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGc                                                >  2:1258687/1‑150 (MQ=255)
                                                                                                       cTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCg                                               <  2:873969/150‑1 (MQ=255)
                                                                                                        tGGTGCCTCCAGAATCACCCCAAGCTCTTGTCGAGACGGCGGCTGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATAGCCACAACAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAACTGCACCAGGCgc                                              <  2:1471600/150‑2 (MQ=255)
                                                                                                            gCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGcc                                          <  2:1161594/150‑1 (MQ=255)
                                                                                                            gCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGcc                                          <  2:771197/150‑1 (MQ=255)
                                                                                                               tCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGa                                       >  1:1060289/1‑150 (MQ=255)
                                                                                                                 cAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATa                                     >  1:635479/1‑150 (MQ=255)
                                                                                                                  aGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATaa                                    <  2:1722413/150‑1 (MQ=255)
                                                                                                                      tCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAagag                                >  2:532596/1‑150 (MQ=255)
                                                                                                                      tCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAagag                                >  1:883182/1‑150 (MQ=255)
                                                                                                                         ccacAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTa                             >  1:329447/1‑150 (MQ=255)
                                                                                                                         ccacAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTa                             >  1:1233534/1‑150 (MQ=255)
                                                                                                                          cacaAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTAc                            >  1:100886/1‑150 (MQ=255)
                                                                                                                          cacaAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTAc                            >  1:526175/1‑150 (MQ=255)
                                                                                                                             aaGCTATTGTAGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTgg                         >  2:279765/1‑150 (MQ=255)
                                                                                                                               gattttGTCGTTTGGGCGGTTGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGaa                       <  1:500647/146‑1 (MQ=255)
                                                                                                                                 tATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGc                     <  2:1599797/150‑1 (MQ=255)
                                                                                                                                  aTTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCa                    >  1:1636853/1‑150 (MQ=255)
                                                                                                                                  aTTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACAGATAGCTGCTGCGAAAACATTGCCACAACAGAGGCTTATTGTGGCAACCGAAAGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCa                    >  1:1624630/1‑150 (MQ=255)
                                                                                                                                   ttGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCAc                   <  1:1532890/150‑1 (MQ=255)
                                                                                                                                   ttGTCGATATGGCGGATGCGGGCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCAc                   <  1:1532908/150‑1 (MQ=255)
                                                                                                                                    tGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCAcc                  <  1:1021172/150‑1 (MQ=255)
                                                                                                                                     gTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCa                 >  2:562948/1‑150 (MQ=255)
                                                                                                                                         atatGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCg             >  2:1590686/1‑150 (MQ=255)
                                                                                                                                         atatGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCg             >  2:1590689/1‑150 (MQ=255)
                                                                                                                                         atatGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCGTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATTCAGCAGGCGGTGCCAGATAAAGAGTTACAGGAAGCACCAAccc             >  1:609154/1‑149 (MQ=255)
                                                                                                                                         atatGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCGGCTGCGAAAACCTTGCCACCTCAGAGGCGGATTGTGGCCACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCCGATAAAGAGTTACTGGACGCACCAAccc             >  1:191212/1‑149 (MQ=255)
                                                                                                                                            tGGCGGATGTGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAg          <  2:621965/150‑1 (MQ=255)
                                                                                                                                            tGGCGGATGTGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAg          <  2:621970/150‑1 (MQ=255)
                                                                                                                                            tGGCGGATGTGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAg          <  2:621963/150‑1 (MQ=255)
                                                                                                                                            tGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAg          <  2:1384352/150‑1 (MQ=255)
                                                                                                                                            tGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAg          <  2:780193/150‑1 (MQ=255)
                                                                                                                                            tGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAg          <  2:780267/150‑1 (MQ=255)
                                                                                                                                            tGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAg          <  2:780291/150‑1 (MQ=255)
                                                                                                                                              gCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGACAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCCCAggg        >  1:1625432/1‑149 (MQ=255)
                                                                                                                                              gCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGt        >  2:1711540/1‑150 (MQ=255)
                                                                                                                                              gCGGATGCGGTCGGTTCCACCAGTCAACAGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGt        >  1:809486/1‑150 (MQ=255)
                                                                                                                                                    gCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGt  <  2:1439017/150‑1 (MQ=255)
                                                                                                                                                    gCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGt  <  1:1132910/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGGTGGAGACATTCTATGCTGGCAGGGTGCCTGTATTGTGCATGATGAATTTAAGACTCAGGCGTTAACCCGCTTGCAAGAAGAATACCCGGATGCTGCCATACTGGTGCATCCAGAATCACCACAAGCTATTGTCGATATGGCGGATGCGGTCGGTTCCACCAGTCAACTGATCGCTGCTGCGAAAACATTGCCACATCAGAGGCTTATTGTGGCAACCGATCGGGGTATTTTCTACAAAATGCAGCAGGCGGTGCCAGATAAAGAGTTACTGGAAGCACCAACCGCAGGTGAGGGT  >  NC_000913/782651‑782948

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