Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,018,454 G→A 5.9% V799M (GTG→ATG)  ygfK → putative oxidoreductase, Fe‑S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,018,4540GA 5.9% 216.3 / 2.3 68V799M (GTG→ATG) ygfKputative oxidoreductase, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base G (30/34);  new base A (2/2);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

TAACCTTGCGCGTTATGTCGCTTGGCGAACCGGATGAGAAAGGTCGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTATGATCGGCGACGTACAGCGCGGACCATCCTCCATTGTCGCTGCTGTCGGAACCGCGCGTCGGGCGACCGATGCCATCCT  >  NC_000913/3018305‑3018602
                                                                                                                                                     |                                                                                                                                                    
tAACCTTGCGCGTTATGTCGCTTGGCGAACCGGATGAGAAAGGTCGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCg                                                                                                                                                      <  2:243000/150‑1 (MQ=255)
   ccTTGCGCGTTATGTCGCTTGGCGAACCGGATGAGAAAGGTCGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACTGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGc                                                                                                                                                   >  1:118519/1‑150 (MQ=255)
      tGCGCGTTATGTCGCTTGGCGAACCGGATGAGAAAGGTCGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGc                                                                                                                                                >  2:872528/1‑150 (MQ=255)
                     ttGGCGAACCGGATGAGAAAGGTCGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCt                                                                                                                                 <  2:186630/150‑1 (MQ=255)
                        gCGAACCGGATGAGAAAGGTCGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAACGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGc                                                                                                                              <  2:940233/150‑1 (MQ=255)
                          gAACCGGATGAGAAAGGTCGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCa                                                                                                                            <  2:18949/150‑1 (MQ=255)
                           aaCCGGATGAGAAAGGTCGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAg                                                                                                                           >  1:628810/1‑150 (MQ=255)
                            aCCGGATGAGAAAGGTCGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGTTGGCCAGa                                                                                                                          >  1:433531/1‑150 (MQ=255)
                             ccGGATGAGAAAGGTCGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAATACGGCTGGCCAGAc                                                                                                                         >  1:735894/1‑150 (MQ=255)
                                gATGAGAAAGGTCGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTc                                                                                                                      <  1:346351/150‑1 (MQ=255)
                                  tGAGAAAGGTCGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGTATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGa                                                                                                                    <  2:1389036/150‑1 (MQ=255)
                                  tGAGAAAGGTCGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGa                                                                                                                    >  2:322873/1‑150 (MQ=255)
                                  tGAGAAAGGTCGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGa                                                                                                                    <  2:1000247/150‑1 (MQ=255)
                                    agaAAGGTCGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGAcc                                                                                                                  >  2:1077279/1‑150 (MQ=255)
                                        aGGTCGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATaa                                                                                                              <  1:760364/150‑1 (MQ=255)
                                        aGGTCGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATaa                                                                                                              <  1:755299/150‑1 (MQ=255)
                                          gtcgtcgtcgtcCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCt                                                                                                                                               <  1:1316487/115‑1 (MQ=255)
                                          gtcgtcgtcgtcCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCt                                                                                                                                               >  2:1316487/1‑115 (MQ=255)
                                          gtcgtcgtcgtcCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATg                                                                                                            >  1:437897/1‑150 (MQ=255)
                                                 tcgtcCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGGGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAAc                                                                                                     <  1:1077279/150‑1 (MQ=255)
                                                 tcgtcCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAAc                                                                                                     <  2:529924/150‑1 (MQ=255)
                                                 tcgtcCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAAc                                                                                                     <  2:886609/150‑1 (MQ=255)
                                                  cgtcCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACt                                                                                                    <  2:1174548/150‑1 (MQ=255)
                                                    tcCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCg                                                                                                  >  2:501768/1‑150 (MQ=255)
                                                     ccGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGt                                                                                                 >  2:19794/1‑150 (MQ=255)
                                                       ggTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTct                                                                                               <  1:459554/150‑1 (MQ=255)
                                                        gTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTctg                                                                                              >  2:194031/1‑150 (MQ=255)
                                                        gTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTctg                                                                                              <  1:691961/150‑1 (MQ=255)
                                                          tGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTctgac                                                                                            <  1:501768/150‑1 (MQ=255)
                                                            aaaCCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTctgactg                                                                                          >  2:999176/1‑150 (MQ=255)
                                                             aaCCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTctgactga                                                                                         <  2:437897/150‑1 (MQ=255)
                                                              aCCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCtgactgac                                                                                        >  1:429125/1‑150 (MQ=255)
                                                              aCCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAAAAGCAGGATACTGACGCCATGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGCCCAGACGTCGACCATAATGCCGAAACTCGTCtgactgac                                                                                        >  1:670113/1‑150 (MQ=255)
                                                                     aaaCAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTa                                                                                 <  2:646823/150‑1 (MQ=255)
                                                                      aaCAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTAt                                                                                <  2:1135923/150‑1 (MQ=255)
                                                                      aaCAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTAt                                                                                <  2:628810/150‑1 (MQ=255)
                                                                            aaCACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCATGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTATGATCg                                                                           >  1:557587/1‑149 (MQ=255)
                                                                            aaCACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCATGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTATGATCg                                                                           >  1:562267/1‑149 (MQ=255)
                                                                            aaCACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCATGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTATGATCg                                                                           <  2:562267/149‑1 (MQ=255)
                                                                            aaCACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCATGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTATGATCg                                                                           <  2:557587/149‑1 (MQ=255)
                                                                                  gCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTATGATCGGCGACGt                                                                    >  1:1248819/1‑150 (MQ=255)
                                                                                  gCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTATGATCGGCGACGt                                                                    >  2:1300864/1‑150 (MQ=255)
                                                                                     tGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTATGATCGGCGACGTACa                                                                 >  1:677681/1‑150 (MQ=255)
                                                                                      gTAGACAGCCTGATCACAGCCATTTGTGAACAGCAGTATACTGAATCCCTGAATGCGATGGGCGTGCCTCTGTACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTATGATCTGCTACGTACAg                                                                <  2:666454/150‑1 (MQ=255)
                                                                                         gACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTATGATCGGCGACGTACAgcgc                                                             >  2:224440/1‑150 (MQ=255)
                                                                                         gACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTATGATCGGCGACGTACAgcgc                                                             >  2:953577/1‑150 (MQ=255)
                                                                                           cAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTATGATCGGCGACGTACAGCGCgg                                                           <  1:657362/150‑1 (MQ=255)
                                                                                           cAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTATGATCGGCGACGTACAGCGCgg                                                           <  1:645126/150‑1 (MQ=255)
                                                                                             gCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTATGATCGGCGACGTACAGCGCGGAc                                                         >  1:1051514/1‑150 (MQ=255)
                                                                                             gCCTGATCACCGCCATTGGTGAACAACAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTATGATCGGCGACGTACAGCGCGGAc                                                         <  1:912874/150‑1 (MQ=255)
                                                                                                    cACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTATGATCGGCGACGTACAGCGCGGa                                                          <  1:800445/142‑1 (MQ=255)
                                                                                                    cACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTATGATCGGCGACGTACAGCGCGGACCAtcctc                                                  >  2:1171811/1‑150 (MQ=255)
                                                                                                    cACCGCCATTGGTGAACAGCAGGATACTGAAGCACTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTATGATCGGCGACGTACAGcgcgca                                                          >  2:800445/1‑140 (MQ=255)
                                                                                                      ccGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTATGATCGGCGACGTACAGCGCGGACCATCCTCCa                                                <  1:872528/150‑1 (MQ=255)
                                                                                                            ttGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGTTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTATGATCGGCGACGTACAGCGCGGACCATCCTCCATTGTCg                                          <  2:433531/150‑1 (MQ=255)
                                                                                                             tGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTATGATCGGCGACGTACAGCGCGGACCATCCTCCATTGTCgc                                         <  1:311960/150‑1 (MQ=255)
                                                                                                                   aCAGCAGGATACTGACGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCGTTATGATCGGCGACGTACAGCGCGGACCATCCTCCATTGTCGCTGCTGt                                   <  1:158042/150‑1 (MQ=255)
                                                                                                                    cagcagGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTATGATCGGCGACGTACAGCGCGGACCATCCTCCATTGTCGCTGCTGTc                                  >  2:563701/1‑150 (MQ=255)
                                                                                                                         ggATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTATGATCGGCGACGTACAGCGCGGACCATCCTCCATTGTCGCTGCTGTCGGAAc                             <  1:1195631/150‑1 (MQ=255)
                                                                                                                           aTACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAGACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTATGATCGGCGACGTACAGCGCGGACCATCCTCCATTGTCGCTGCTGTCGGAACcg                           >  1:1089894/1‑150 (MQ=255)
                                                                                                                             aCTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTATGATCGGCGACGTACAGCGCGGACCATCCTCCATTGTCGCTGCTGTCGGAACcgcg                         <  1:953577/150‑1 (MQ=255)
                                                                                                                               tGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTATGATCGGCGACGTACAGCGCGGACCATCCTCCATTGTCGCTGCTGTCGGAACcgcgcg                       >  1:1379212/1‑150 (MQ=255)
                                                                                                                                gAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGa                                                                                                        <  1:97937/68‑1 (MQ=255)
                                                                                                                                gAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGa                                                                                                        >  2:97937/1‑68 (MQ=255)
                                                                                                                                gAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTATGATCGGCGACGTACAGCGCGGACCATCCTCCATTGTCGCTGCTGTCGGAACCGCGCGt                      >  2:421210/1‑150 (MQ=255)
                                                                                                                                 aaGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTATGATCGGCGACGTACAGCGCGGACCATCCTCCATTGTCGCTGCTGTCGGAACCGCGCGTc                     <  2:118519/150‑1 (MQ=255)
                                                                                                                                        gAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTATGATCGGCGACGTACAGCGCGGACCATCCTCCATTGTCGCTGCTGTCGGAACCGCGCGTCGGGCGAc              >  1:260572/1‑150 (MQ=255)
                                                                                                                                           tGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTATGATCGGCGACGTACAGCGCGGACCATCCTCCATTGTCGCTGCTGTCGGAACCGCGCGTCGGGCGACCGa           >  2:115093/1‑150 (MQ=255)
                                                                                                                                                    cGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTATGATCGGCGACGTACAGCGCGGACCATCCTCCATTGTCGCTGCTGTCGGAACCGCGCGTCGGGCGACCGATGCCATCCt  <  2:1051514/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TAACCTTGCGCGTTATGTCGCTTGGCGAACCGGATGAGAAAGGTCGTCGTCGTCCGGTTGAAACCAATGAAACAGTAACACTGCTTGTAGACAGCCTGATCACCGCCATTGGTGAACAGCAGGATACTGAAGCCCTGAATGCGATGGGCGTGCCGCTGGACAAAAACGGCTGGCCAGACGTCGACCATAATGGCGAAACTCGTCTGACTGACGTCTTTATGATCGGCGACGTACAGCGCGGACCATCCTCCATTGTCGCTGCTGTCGGAACCGCGCGTCGGGCGACCGATGCCATCCT  >  NC_000913/3018305‑3018602

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: