Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 126,147 T→G 6.9% A151A (GCT→GCG aceF → pyruvate dehydrogenase, dihydrolipoyltransacetylase component E2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913126,1470TG 6.9% 234.4 / 3.5 79A151A (GCT→GCGaceFpyruvate dehydrogenase, dihydrolipoyltransacetylase component E2
Reads supporting (aligned to +/- strand):  ref base T (24/49);  new base G (4/2);  total (28/51)
Fisher's exact test for biased strand distribution p-value = 1.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTG  >  NC_000913/126000‑126295
                                                                                                                                                   |                                                                                                                                                    
tGCGGCGGCAAAAGACGGTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGGTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTcc                                                                                                                                                    <  1:734550/150‑1 (MQ=255)
 gcggcgGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCg                                                                                                                                                   <  2:1572653/150‑1 (MQ=255)
  cggcggCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGt                                                                                                                                                  >  1:1347089/1‑150 (MQ=255)
       gcAAAAGACGTTACCGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGGGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCt                                                                                                                                             <  2:131219/150‑1 (MQ=255)
       gcAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCt                                                                                                                                             <  1:753692/150‑1 (MQ=255)
       gcAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCt                                                                                                                                             <  2:1319638/150‑1 (MQ=255)
        cAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTg                                                                                                                                            <  2:1659076/150‑1 (MQ=255)
        cAAAAGACGTTAACGTTCCGGAAATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCATCTATGGAAGTTCCGGCTCCGtatacgg                                                                                                                                            >  1:249200/1‑144 (MQ=255)
          aaaGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGc                                                                                                                                          <  1:573324/150‑1 (MQ=255)
            aGACGTTAACGTTCCTGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTTGTGAAAGTTTGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCAc                                                                                                                                        <  2:71517/150‑1 (MQ=255)
                gTTAAGGTTCCGGAAATCTGCAGCGACGAAGTTGATGTGACCGAAATGCTGTTGAAAGTTGGCTATAATGTTGAAGCTGATCAGTCGCTGATCACCGAATAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTg                                                                                                                                    <  2:883425/150‑1 (MQ=18)
                    aCGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAATTTGGCGATAAAGTTGAAGCTGTACAGTCGCTGATCTCCGTAGAAGGCGACAATGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAag                                                                                                                                <  1:1274932/150‑1 (MQ=255)
                    aCGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCTCCGTAGAAGGCGACAAGGCTTCTATGGATGTTCCGGCTCCGTTTGCTGGCACCGTGAAag                                                                                                                                <  1:1119318/150‑1 (MQ=255)
                      gTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCGCCGTTTGCTGGCACCGTGAAagag                                                                                                                              >  1:1124528/1‑150 (MQ=255)
                        tCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGAt                                                                                                                            <  2:814621/150‑1 (MQ=255)
                                   gCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGTGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACg                                                                                                                 <  1:11576/150‑1 (MQ=255)
                                   gCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACg                                                                                                                 <  1:130237/150‑1 (MQ=255)
                                   gCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACg                                                                                                                 <  2:1103521/150‑1 (MQ=255)
                                       cgacgaAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTggg                                                                                                             <  2:428633/150‑1 (MQ=255)
                                         acgaAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTg                                                                                                           >  1:1407975/1‑150 (MQ=255)
                                          cgaAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGa                                                                                                          >  1:704007/1‑150 (MQ=255)
                                           gaAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGAc                                                                                                         >  2:1625801/1‑150 (MQ=255)
                                            aaGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACa                                                                                                        >  1:1573914/1‑150 (MQ=255)
                                             aGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCTCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGaag                                                                                                                                 <  2:455390/124‑2 (MQ=14)
                                             aGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGaag                                                                                                                                 >  1:455390/1‑123 (MQ=14)
                                               ttGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAg                                                                                                     <  2:1581291/150‑1 (MQ=255)
                                               ttGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGGGACAAAg                                                                                                     <  1:842193/150‑1 (MQ=255)
                                               gtGAAGTGAACGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGGCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCAGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAg                                                                                                     <  1:313913/149‑1 (MQ=255)
                                                tGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAgt                                                                                                    <  1:536548/150‑1 (MQ=255)
                                                  aaGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCATAGTGAACGTGGGTGACAAAgtgt                                                                                                  <  1:1271492/150‑1 (MQ=255)
                                                  aaGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAgtgt                                                                                                  >  2:1634892/1‑150 (MQ=255)
                                                   aGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTc                                                                                                 >  1:336108/1‑150 (MQ=255)
                                                   aGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTc                                                                                                 <  1:814688/150‑1 (MQ=255)
                                                      gACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTAc                                                                                              >  2:1408708/1‑150 (MQ=255)
                                                        ccGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCg                                                                                            <  1:1417169/150‑1 (MQ=255)
                                                            aaTCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTc                                                                                        <  1:1565234/150‑1 (MQ=255)
                                                              tCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGc                                                                                      <  2:996120/150‑1 (MQ=255)
                                                                cTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGGCGCGGATCACCGTAGAAGTCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAACGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTg                                                                                    <  2:1059090/150‑1 (MQ=255)
                                                                  ggTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGAt                                                                                  <  2:1475171/150‑1 (MQ=255)
                                                                    tGAAAGTTTGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGATCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTa                                                                                <  2:1670364/150‑1 (MQ=255)
                                                                    tGAAAGTTGTCTATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTTAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTa                                                                                <  1:1196217/150‑1 (MQ=255)
                                                                    tGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTa                                                                                >  2:1413651/1‑150 (MQ=255)
                                                                             gCGATAATGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCg                                                                       <  2:1723730/150‑1 (MQ=255)
                                                                               gATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGaa                                                                     <  1:413878/150‑1 (MQ=255)
                                                                                aTAAAGTTTAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAg                                                                    <  1:1027060/150‑1 (MQ=255)
                                                                                 aaaaGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGGGAACGTGGGGGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGt                                                                   <  2:1573253/149‑1 (MQ=255)
                                                                                    aGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGcag                                                                                                                                            >  1:628167/1‑72 (MQ=11)
                                                                                    aGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGcag                                                                                                                                            <  2:628167/74‑3 (MQ=11)
                                                                                     gTTGAAGCTGAACAGGCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGGCTTCGAAGTcgcg                                                               <  2:1342477/150‑1 (MQ=255)
                                                                                       tGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGAGTATGGTCTTCGAAGTCGCggg                                                             >  1:1711512/1‑150 (MQ=255)
                                                                                       tGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGcgtc                                                                                                                                       >  1:821085/1‑73 (MQ=11)
                                                                                             cgAACAGTCGCTGATCACCGTAGAAGGCGACAACGCCTCTATGGAAGTTCCGGCTCCGGTTGCTGGGAACGTGAAAGAGATCAAAGTGAACGTGGGTGAGAAAGAGTCTACCGGCGCGCTGATTATGGTATTCGATTTCGCGGGTGGAGc                                                       >  1:485803/2‑150 (MQ=21)
                                                                                              gAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCa                                                      <  1:1423970/150‑1 (MQ=255)
                                                                                                   gTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTg                                                                                                                                    >  2:1429994/1‑67 (MQ=255)
                                                                                                   gTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTg                                                                                                                                    <  1:1429994/67‑1 (MQ=255)
                                                                                                      gCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGc                                              <  2:680477/150‑1 (MQ=255)
                                                                                                      gCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGGGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGc                                              <  1:131231/150‑1 (MQ=255)
                                                                                                      gCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTGCCGGCTCCGTTTGCTCGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGGCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGc                                              <  1:111666/150‑1 (MQ=255)
                                                                                                       cTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGcc                                             >  2:1157690/1‑149 (MQ=255)
                                                                                                        tGATCACCGTAGAAGTCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCGGGCACCGTGAAAGAGATCAAAGTCAACGTGGGCGACAAAGAGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAg                                            <  2:1163522/150‑1 (MQ=25)
                                                                                                        tGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAg                                            >  1:827455/1‑150 (MQ=255)
                                                                                                             aCCGTAGAAGTCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTTGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCg                                       <  2:767753/150‑1 (MQ=255)
                                                                                                                gTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGcc                                    >  2:269089/1‑150 (MQ=255)
                                                                                                                 tAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCg                                   >  2:1281224/1‑150 (MQ=255)
                                                                                                                 caGAAGGCGACAAGGCGTCTATGGAAGTTCCGGCTCCGATTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCg                                   <  1:1042670/149‑1 (MQ=255)
                                                                                                                  aGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGc                                  <  2:1347089/150‑1 (MQ=255)
                                                                                                                   gAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCt                                 <  2:1336284/150‑1 (MQ=255)
                                                                                                                      ggCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTaaa                              >  1:1422190/1‑150 (MQ=255)
                                                                                                                       gCGACAAGGCGTCTATGGAAGTTCCGGCTCCGTTTGCTGGCCCCGGGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAAc                             <  2:1124528/150‑1 (MQ=255)
                                                                                                                        cGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACa                            >  2:3066/1‑150 (MQ=255)
                                                                                                                         gACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAg                           >  2:1217957/1‑150 (MQ=255)
                                                                                                                           cAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGa                         >  2:903370/1‑150 (MQ=255)
                                                                                                                               gCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAagca                     <  2:485803/150‑1 (MQ=255)
                                                                                                                                    tATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAAGCAGCTcc                >  2:1338372/1‑150 (MQ=255)
                                                                                                                                       ggAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAAGCAGCTCCGGc             <  2:704007/150‑1 (MQ=255)
                                                                                                                                          aGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTAGCGGGGGCAGCAGGAGCGGCAGCTCCGGCCTCTAAGGACGCAGCAGCGCCGGCAGc          >  2:1188313/1‑150 (MQ=255)
                                                                                                                                                 gCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCt   <  1:1217957/150‑1 (MQ=255)
                                                                                                                                                 gCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCt   <  1:1413651/150‑1 (MQ=255)
                                                                                                                                                  cTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTg  >  2:499979/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGCGGCGGCAAAAGACGTTAACGTTCCGGATATCGGCAGCGACGAAGTTGAAGTGACCGAAATCCTGGTGAAAGTTGGCGATAAAGTTGAAGCTGAACAGTCGCTGATCACCGTAGAAGGCGACAAGGCTTCTATGGAAGTTCCGGCTCCGTTTGCTGGCACCGTGAAAGAGATCAAAGTGAACGTGGGTGACAAAGTGTCTACCGGCTCGCTGATTATGGTCTTCGAAGTCGCGGGTGAAGCAGGCGCGGCAGCTCCGGCCGCTAAACAGGAAGCAGCTCCGGCAGCGGCCCCTG  >  NC_000913/126000‑126295

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