Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,613,9040AC40.0% 102.5 / 14.0 60V267G (GTA→GGA) sadsuccinate semialdehyde dehydrogenase, NAD(P)+‑dependent
Reads supporting (aligned to +/- strand):  ref base A (27/9);  new base C (0/24);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 1.09e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGTAAATCAAAACGAGCCATTGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGC  >  NC_000913/1613755‑1614052
                                                                                                                                                       |                                                                                                                                                    
caCGTAAATCAAAACGAGCCATTGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATa                                                                                                                                                      <  2:211034/150‑1 (MQ=255)
 aCGTAAATCAAAACGAGCCATTGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATAc                                                                                                                                                     <  1:607761/150‑1 (MQ=255)
  cGTAAATCAAAACGAGCCATTGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTTTCGCTGCACATACCt                                                                                                                                                   >  1:324509/1‑151 (MQ=255)
  cGTAAATCAAAACGAGCCATTGGTCCGAGAGCGTTCGCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAGCGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATa                                                                                                                                                      <  1:214380/148‑1 (MQ=255)
     aaaTCAAAACGAGCCATTGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTc                                                                                                                                                >  1:389847/1‑151 (MQ=255)
      aaTCAAAACGAGCCATGGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTcc                                                                                                                                               <  1:400361/151‑1 (MQ=255)
       aTCAAAACGAGCCATTGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCg                                                                                                                                              <  2:981103/151‑1 (MQ=255)
       aTCAAAACGAGCCATTGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCg                                                                                                                                              <  2:302163/151‑1 (MQ=255)
       aTCAAAACGAGCCATTGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCg                                                                                                                                              >  2:432171/1‑151 (MQ=255)
       aTCAAAACGAGCCATTGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCg                                                                                                                                              <  1:274743/151‑1 (MQ=255)
           aaaCGAGCCATTGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTAt                                                                                                                                          >  2:648109/1‑151 (MQ=255)
              cGAGCCATTGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCt                                                                                                                                       >  1:640601/1‑151 (MQ=255)
               gAGCCATTGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTg                                                                                                                                      >  1:420/1‑151 (MQ=255)
                      tcccccGAGCGCGTTCTCTTCGTCACGAGGATCGCCCATCTCCAAGGCCGCCGCAGCTGCCACAAACCGCTCGGTAAATGCCGCAGCAATTCCCTCTCCGATAATAAAGCGTTTTGCCGCTGCACATACCTGCCCGGTATTCTCATAAc                                                                                                                                 >  1:955908/5‑149 (MQ=255)
                        gTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTcc                                                                                                                             <  2:1009664/151‑1 (MQ=255)
                          ccGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCg                                                                                                                            <  2:823925/150‑1 (MQ=255)
                                gCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTAccg                                                                                                                     >  1:764704/1‑151 (MQ=255)
                                  gTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTAccgcc                                                                                                                   >  1:457375/1‑151 (MQ=255)
                                  gTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTAccgcc                                                                                                                   <  2:974748/151‑1 (MQ=255)
                                     ctctTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCt                                                                                                                >  1:292599/1‑151 (MQ=255)
                                      tctTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCtt                                                                                                               >  2:102199/1‑151 (MQ=255)
                                        ttCGCCACGGGGATCGACCATTTTCAAGGCTGCCGCAGCTGCCGCAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGCATTCTGATAACGTCCGGCTACCGCCGCtct                                                                                                              >  1:409236/1‑148 (MQ=255)
                                              aCGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGcc                                                                                                       >  2:123931/1‑151 (MQ=255)
                                               cGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGcc                                                                                                       >  1:803661/1‑150 (MQ=255)
                                                ggggATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATCCCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCa                                                                                                      <  2:289123/150‑1 (MQ=255)
                                                ggggATCGCCCATTTCCAAGGTTGCCGCACCGGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCCCATACCTGTCCGGTATTCTGATACCGTCCGGCTACCGCCGCTTTCACCg                                                                                                         <  2:226314/147‑1 (MQ=255)
                                                ggggACCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATCCCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCa                                                                                                      <  2:923260/150‑1 (MQ=255)
                                                     tCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTcc                                                                                                 >  1:340307/1‑150 (MQ=255)
                                                                        ccccggCTCCCACAAACCGTTGGGTAAATGCCAAAGAAATTCCCTCTTGGTAAAAAAAGCGTTTTGCCGCTCCACATCCCTGCCCGGTATTCTGATAACGTCGGGCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGcac                                                                             <  2:646193/146‑1 (MQ=255)
                                                                        ccGCGGCTGCCACAAACCGTTGGGTAAAGGCCGAAGAAATCCCCTCTTCGAAAATAAAGCTTTTTGCCGCTGCACATCCCTGCCCGGTATTCTGAAAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAg                                                                                <  2:961272/148‑1 (MQ=255)
                                                                         cGCAGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGcac                                                                             >  2:59494/1‑150 (MQ=255)
                                                                           tactcgCCACAAACGTCTGGTAAATTCCCAAAGAAATCCCCTTTTCAAAAAAAAAGCTTTTTGCCTCCGCTCTTCCCTGCTGGGTATTCTTATACCTTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCTTTAAg                                                                                <  1:740726/140‑1 (MQ=255)
                                                                            aGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAAt                                                                          >  1:763567/1‑150 (MQ=255)
                                                                                 ccACAAAACGTTCGGTAAATGCCGAAGAAATTCCCTCTTCGATAATAAAGCTTTTTGCCGCTGCACATCCCTTCCCGGTATTCTGATACCGTCCGGTTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATaaa                                                                       <  2:329644/148‑1 (MQ=255)
                                                                                 ccAAAAAACGTTCGATAATTCCCGAACCATTTCCCTCTTCAAAAAAAACGCTTTTTGCCGCTGCCCACCCCTCTCCGATATTCTAAAACCCACCGACACCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATa                                                                         <  1:201988/146‑1 (MQ=255)
                                                                                     aaactgtgcGGTAAATGCCGAGACAATTCCCTCTTGAATAAAAAAGTGTTTGGCCTCTCCACATCCCTGCCCGTTATTCTAATGAGTGCGGGCAACCGCCGGTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGAt                                                                  <  1:610540/141‑1 (MQ=255)
                                                                                      aaaCCTTCGGTAAATGCCCGAAGAAATCCCCTCTCCAAAAAAAAAGATTTTTGCCCCTCCACAACCCTTCCCGTTTTTATAAGAACGTCGGGCACCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGa                                                                   <  1:714394/148‑1 (MQ=255)
                                                                                         cGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATCCCTGTCCGGTATTCTGATACCGTCCGGCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGaa                                                             <  2:700541/150‑1 (MQ=255)
                                                                                         cGTTCGGAAAATGCCGAAGAATTTCCCTCTTATAAAAAAAAGCTTTTTGCCGCTGCACATCCCTGCCCGTTTTTGTGAAAACGTCCGCATACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGa                                                              <  1:458063/149‑1 (MQ=255)
                                                                                                 aaaTGCCAAAGCAATCCCCTCTCCAATAATAAAGCGTTTTGCCGCTCCACATCCCTGCCCGGTATTCTGAAAACCTCCGGCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCa                                                      <  2:760511/149‑1 (MQ=255)
                                                                                                 aaaTGCCAAAGAAATCCCCTCTTCAATAAAAAAGCTTTTTGCCGCTCCACATCCCTGCCCGGTATTCTGATAACTCCCGGCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCa                                                      <  2:194238/149‑1 (MQ=255)
                                                                                                     tGCCGAAGCAAGTCCCTCTTCGATAATAAAGCGTTTGGCCGCTGCACAGACCTGTACGGTATTCTGATACAGTCTGGCTACCGCCGCTTTCACAGCCAGATCCAGACCGTCATCGTTAAGCACAATAAACGTGGCCGAAGCGCCCAGTTcc                                                 >  1:399748/1‑151 (MQ=255)
                                                                                                     atgccaaggaatttCCCTCTTATTTAAAAAGGCTTTTTGCCGCTCCCCTTCCCTGCCCGGAATTCTGAAAACGTCCGGCTCCCGCCTCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTc                                                  <  2:562349/139‑1 (MQ=255)
                                                                                                              cAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGc                                          >  2:764621/1‑150 (MQ=255)
                                                                                                               aaTCCCCTCTTCGAAAAAAAAGCTTTTTGCCGCTGCACATCCCGTTCCGGTATTCTAAAACCGCCCGTCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCAt                                        <  2:69542/151‑1 (MQ=255)
                                                                                                                   cccTCTTGGAAAATAAGGCGTTTGGCCGCTCCCCATCCCTGCCCGTTATTCTGATAACGCCCGCCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCAt                                        <  2:86899/147‑1 (MQ=255)
                                                                                                                    ccTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCAtttttt                                   >  1:1136423/1‑151 (MQ=255)
                                                                                                                      tctTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCAtttttt                                   >  2:154924/1‑149 (MQ=255)
                                                                                                                        tccGATAAAAAAGCTTTTTGCCGCTGCACATCCCTGTCCGGTTTTCTGAAAACGCCCGGCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTc                                  <  2:930088/146‑1 (MQ=255)
                                                                                                                          caaaaaaaaaGCTTTTTCCCGCTCCACATCCCTGCCGGGTATTCTAAAAACGCCCGGCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTg                              <  2:1062071/143‑1 (MQ=255)
                                                                                                                            ataataAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCg                           >  1:28550/1‑151 (MQ=255)
                                                                                                                              aataaGGCTTTTTGCCGCTGCACATCCCTGTCGGTTATTCTAATACCGCCCGCCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCt                         <  2:435400/151‑1 (MQ=255)
                                                                                                                              aataaAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCt                         >  2:62622/1‑151 (MQ=255)
                                                                                                                              aaaaaaGCGTTTTGCCGCTGCACATCCCTGCCGGTTTTTCTAAAACCGCCCGGCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCg                           <  2:1039528/146‑1 (MQ=255)
                                                                                                                                       ttttGCCGCTCCACATCCTTTTCGATTTTTCCAATAACGTCGGCCTCCCGCCGTTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAg                     <  1:585347/146‑1 (MQ=255)
                                                                                                                                       ttttCCCGCTCCACATCCCTGCCGGGTATTCTAATAACGCCCGGCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGcc                   <  2:324699/148‑1 (MQ=255)
                                                                                                                                        tttGCCGCTGCACATCCCTGCCGGGTTTTCTAATACCGCCCGGCTCCCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCtgt                <  2:406583/150‑1 (MQ=255)
                                                                                                                                         ttGCCGCTGCACATCCCTGTCCGGTATTCTGATAACGCCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCtgtg               <  2:328209/150‑1 (MQ=255)
                                                                                                                                               cTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAAt        >  1:304951/1‑151 (MQ=255)
                                                                                                                                                tGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAAt        >  1:473657/1‑150 (MQ=255)
                                                                                                                                                     aTACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGc  >  1:808642/1‑151 (MQ=255)
                                                                                                                                                     aTACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGc  >  2:262392/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
CACGTAAATCAAAACGAGCCATTGGTCCGAGAGCGTTCTCTTCGTCACGGGGATCGCCCATTTTCAAGGCTGCCGCAGCTGCCACAAAACGTTCGGTAAATGCCGAAGCAATTCCCTCTTCGATAATAAAGCGTTTTGCCGCTGCACATACCTGTCCGGTATTCTGATAACGTCCGGCTACCGCCGCTTTCACCGCCAGTTCCAGATCGGCATCGTTAAGCACAATAAACGGATCCGAACCGCCCAGTTCCAGTACGCATTTTTTCAGTGCCGCTCCAGCCTGTGCGCCAATAGCCGC  >  NC_000913/1613755‑1614052

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