Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 378,000 G→A 5.8% N123N (AAC→AAT)  frmB ← S‑formylglutathione hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913378,0000GA 5.8% 230.4 / 3.6 69N123N (AAC→AAT) frmBS‑formylglutathione hydrolase
Reads supporting (aligned to +/- strand):  ref base G (38/27);  new base A (2/2);  total (40/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-01

AGACAATGGGCGAAAACGCCGAGACGCTGACATATTCATCTGGGTTACGTAACGCCAGCACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGCCCGAGATCGTAACGGTCAGCATCTGCGACATGACTGCCTCGCGGACTGGTGTCCGGCGCAACAACAATAA  >  NC_000913/377852‑378148
                                                                                                                                                    |                                                                                                                                                    
agaCAATGGGCGAAAACGCCGAGACGCTGACATATTCATCTGGGTTACGTAACGCCAGCACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCg                                                                                                                                                      >  2:596107/1‑149 (MQ=255)
    aaTGGGCGAAAACGCCGAGACGCTGACATATTCATCTGGGTTACGTAACGCCAGCACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGc                                                                                                                                                  >  3:6380/1‑149 (MQ=255)
          cGAAAACGCCGAGACGCTGACATATTCATCTGGGTTACGTAACGCCAGCACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGatat                                                                                                                                            >  3:610564/1‑149 (MQ=255)
             aaaCGCCGAGCCGCTGACATATTCATCTGGGTTACGTAACGCCAGCACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGt                                                                                                                                         >  2:21491/1‑149 (MQ=255)
               aCGCCGAGACGCTGACATATTCATCTGGGTTACGTAACGCCAGCACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCa                                                                                                                                       >  1:227813/1‑149 (MQ=255)
                    gagaCGCTGACATATTCATCTGGGTTACGTAACGCCAGCACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACAt                                                                                                                                  <  2:743754/149‑1 (MQ=255)
                       aCGCTGACATATTCATCTGGGTTACGTAACGCCATCACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACAtttt                                                                                                                               <  3:655582/149‑1 (MQ=255)
                       aCGCTGACATATTCATCTGGGTTACGTAACGCCAGCACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGCGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACAtttt                                                                                                                               <  4:576331/149‑1 (MQ=255)
                               atatTCATCTGGGTTACGTAACGCCAGCACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGtt                                                                                                                       <  1:596107/149‑1 (MQ=255)
                                tatTCATCTGGGTTACGTAACGCCAGCACCAGCGCGCCCAGCCCGCCCATAGAATGACAAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTc                                                                                                                      >  3:147599/1‑149 (MQ=255)
                                 atTCATCTGGGTTACGTAACGCCAGCACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGt                                                                                                                        >  4:175061/1‑146 (MQ=255)
                                 atTCATCTGGGTTACGTAACGCCAGCACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGt                                                                                                                        <  3:175061/146‑1 (MQ=255)
                                  ttCATCTGGGTTACGTAACGCCAGCACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCAt                                                                                                                    >  1:148355/1‑149 (MQ=255)
                                     aTCTGGGTTACGTAACGCCAGCACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCa                                                                                                                     <  4:144003/145‑1 (MQ=255)
                                     aTCTGGGTTACGTAACGCCAGCACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCa                                                                                                                     >  3:144003/1‑145 (MQ=255)
                                      tCTGGGTTACGTAACGCCAGCACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCa                                                                                                                >  2:655353/1‑149 (MQ=255)
                                         gggTTACGTAACGCCAGCACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACgg                                                                                                             >  4:193933/1‑149 (MQ=255)
                                          ggTTACGTAACGCCAGCACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGgc                                                                                                            >  3:103716/1‑149 (MQ=255)
                                          ggTTACGTAACGCCAGCACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGgc                                                                                                            <  1:497823/149‑1 (MQ=255)
                                            ttACGTAACGCCAGCACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGt                                                                                                                                         <  3:256000/118‑1 (MQ=255)
                                            ttACGTAACGCCAGCACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGt                                                                                                                                         >  4:256000/1‑118 (MQ=255)
                                            acctgcaaCGCCAGCACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGgcgc                                                                                                          >  3:313341/7‑149 (MQ=255)
                                              aCGTAACGCCAGCACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCtt                                                                                                        <  1:166668/149‑1 (MQ=255)
                                                gTAACGCCAGCACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGc                                                                                                      >  1:743672/1‑149 (MQ=255)
                                                  aaCGCCAGCACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGt                                                                                                    >  1:697004/1‑149 (MQ=255)
                                                     gCCAGCACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTcgc                                                                                                 <  1:722907/149‑1 (MQ=255)
                                                     gCCAGCACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCATTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTcgc                                                                                                 >  2:578379/1‑149 (MQ=255)
                                                     gCCAGCACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCATTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCACTTGCGTcgc                                                                                                 >  4:597575/1‑149 (MQ=255)
                                                      ccAGCACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTcgcg                                                                                                <  1:272266/149‑1 (MQ=255)
                                                       cAGCACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGt                                                                                               <  2:148355/149‑1 (MQ=255)
                                                         gCACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTc                                                                                             <  3:303063/149‑1 (MQ=255)
                                                          cACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCa                                                                                            <  4:610564/149‑1 (MQ=255)
                                                           aCCAGCGCGCCCAGCCCGCCCATAGAATGACAAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAg                                                                                           <  4:147599/149‑1 (MQ=255)
                                                               gcgcgcCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTaa                                                                                       >  1:537889/1‑149 (MQ=255)
                                                                    cccAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAAccc                                                                                  <  2:295049/149‑1 (MQ=255)
                                                                     ccAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCg                                                                                 >  2:406775/1‑149 (MQ=255)
                                                                              ccATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGc                                                                        >  4:361786/1‑149 (MQ=255)
                                                                              ccATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGc                                                                        <  4:258610/149‑1 (MQ=255)
                                                                              ccATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGc                                                                        >  2:251883/1‑149 (MQ=255)
                                                                              ccATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGc                                                                        <  1:21491/149‑1 (MQ=255)
                                                                                aTAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGccc                                                                      >  3:218632/1‑149 (MQ=255)
                                                                                aTAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGccc                                                                      <  4:218632/149‑1 (MQ=255)
                                                                                aTAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGccc                                                                      >  3:727039/1‑149 (MQ=255)
                                                                                aTAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGccc                                                                      >  4:560637/1‑149 (MQ=255)
                                                                                  aGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGCCCga                                                                    >  4:336139/1‑149 (MQ=255)
                                                                                     aTGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGCCCGAGAt                                                                 <  4:103716/149‑1 (MQ=255)
                                                                                         ccAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGCCCGAGATCGTa                                                             >  2:648540/1‑149 (MQ=255)
                                                                                          cAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGCCCGAGATCGTaa                                                            >  2:746661/1‑149 (MQ=255)
                                                                                           agagATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAACCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGCCCGAGATCGTAAc                                                           >  3:487335/1‑149 (MQ=255)
                                                                                             agaTAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGCTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGCCCGAGATCGTAACgg                                                         >  1:654280/1‑149 (MQ=255)
                                                                                             agaTAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCATTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGCCCGAGATCGTAACgg                                                         <  1:578379/149‑1 (MQ=255)
                                                                                             agaTAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCATTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCACTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGCCCGAGATCGTAACgg                                                         <  3:597575/149‑1 (MQ=255)
                                                                                               aTAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGCCCGAGATCGTAACGGTc                                                       >  3:575721/1‑149 (MQ=255)
                                                                                               aTAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGCCCGAGATCGTAACGGTc                                                       >  1:14277/1‑149 (MQ=255)
                                                                                                tcgACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAgt                                                                                          <  3:317231/111‑2 (MQ=255)
                                                                                                tAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAgg                                                                                          >  4:317231/1‑113 (MQ=255)
                                                                                                tAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGCCCGAGATCGTAACGGTCa                                                      <  3:84014/149‑1 (MQ=255)
                                                                                                 aGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGCCCGAGATCGTAACGGTCAg                                                     >  1:650742/1‑149 (MQ=255)
                                                                                                 aGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGCCCGAGATCGTAACGGTCAg                                                     >  3:371730/1‑149 (MQ=255)
                                                                                                          ggCCGTTGCCGGAAAATGATGCATCACTGAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGCCCGAGATCGTAACGGTCAGCATCTGCGa                                            <  3:686645/149‑1 (MQ=255)
                                                                                                          ggCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTTGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGCCCGAGAGCGTAACGGTCAGCATCTGCGa                                            >  4:518981/1‑149 (MQ=255)
                                                                                                                            atgCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGCCCGAGATCGTAACGGTCAGCATCTGCGACATGACTGCCTCGCGGAc                          <  4:371730/149‑1 (MQ=255)
                                                                                                                              gCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGCCCGAGATCGTAACGGTCAGCATCTGCGACATGACTGCCTCGCGGACTg                        <  3:336139/149‑1 (MQ=255)
                                                                                                                                aTCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGCCCGAGATCGTAACGGTCAGCATCTGCGACATGACTGCCTCGCGGACTGgt                      <  3:90478/149‑1 (MQ=255)
                                                                                                                                aTCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGCCCGAGATCGTAACGGTCAGCATCTGCGACATGACTGCCTCGCGGACTGgt                      >  4:112876/1‑149 (MQ=255)
                                                                                                                                aTCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGCCCGAGATCGTAACGGTCAGCATCTGCGACATGACTGCCTCGCGGACTGgt                      >  2:222983/1‑149 (MQ=255)
                                                                                                                                aTCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGCCCGAGATCGTAACGGTCAGCATCTGCGACATGACTGCCTCGCGGACTGgt                      >  2:544954/1‑149 (MQ=255)
                                                                                                                                  cACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGCCCGAGATCGTAACGGTCAGCATCTGCGACATGACTGCCTCGCGGACTGgtgt                    >  1:324133/1‑149 (MQ=255)
                                                                                                                                                  tCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGCCCGAGATCGTAACGGTCAGc                                                    >  3:88610/1‑101 (MQ=255)
                                                                                                                                                  tCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGCCCGAGATCGTAACGGTCAGc                                                    <  4:88610/101‑1 (MQ=255)
                                                                                                                                                    gTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGCCCGAGATCGTAACGGTCAGCATCTGCGACATGACTGCCTCGCGGACTGGTGTCCGGCGCAACAACaataa  <  2:14277/149‑1 (MQ=255)
                                                                                                                                                    gTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGCCCGAGATCGTAACGGTCAGCATCTGCGACATGACTGCCTCGCGGACTGGTGTCCGGCGCAACAACaataa  <  4:575721/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGACAATGGGCGAAAACGCCGAGACGCTGACATATTCATCTGGGTTACGTAACGCCAGCACCAGCGCGCCCAGCCCGCCCATAGAATGACCAGAGATAGACTTTTTGGCCGTTGCCGGAAAATGATGCATCACTAAATCCGGCAGCTCGTTGCGGATATAGTCATACATTTTGTAATGTTCATTCCACGGCGCTTGCGTCGCGTTCAGGTAAAACCCGGCACCTTGCCCGAGATCGTAACGGTCAGCATCTGCGACATGACTGCCTCGCGGACTGGTGTCCGGCGCAACAACAATAA  >  NC_000913/377852‑378148

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑