Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,587,225 A→T 7.3% M95L (ATG→TTG)  yhhA → DUF2756 domain‑containing protein YhhA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,587,2250AT 7.3% 167.8 / 4.4 64M95L (ATG→TTG) yhhADUF2756 domain‑containing protein YhhA
Reads supporting (aligned to +/- strand):  ref base A (30/28);  new base T (3/2);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.05e-01

TCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAA  >  NC_000913/3587085‑3587367
                                                                                                                                            |                                                                                                                                              
tCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAAtg                                                                                                                                              >  2:709772/1‑143 (MQ=255)
 cgaGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGCGTTGCAGTAACAGCCGGGGGGGCGAAATCCCGCCCGGAAGCGAAtgc                                                                                                                                             >  2:326796/3‑143 (MQ=255)
  cAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAAtgct                                                                                                                                            >  1:158578/1‑143 (MQ=255)
   aGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCACAtgctg                                                                                                                                           >  2:133280/1‑143 (MQ=255)
     cAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCGGCAAGTGCAGcc                                                                                                                                         >  1:172092/1‑143 (MQ=255)
        aaaGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGCGAGCGACCTAACGCCCGGCAGCACATGCTGCCCaa                                                                                                                                      >  2:167807/1‑143 (MQ=255)
         gaGGGATGTTGAATCAGCTACTGACAACGCAACCAGAATTGCTACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGTAGCGAAATCCCGTCCGGCAGCAAATGCTGCCCAac                                                                                                                                     <  2:468636/142‑1 (MQ=255)
           gggATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCGCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCATTTGCTGccgggcgg                                                                                                                                   <  1:989313/143‑7 (MQ=255)
              aTGCTGAATCCGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCaa                                                                                                                                <  2:569469/143‑1 (MQ=255)
              aTGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCGGCCCGGCAGCAAATGCTGCCCAACACCaa                                                                                                                                >  2:964677/1‑143 (MQ=255)
                  tGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCTCGCCCGGCAGCAAATGCTGCCCAACACCAAcgga                                                                                                                            >  1:273864/1‑142 (MQ=255)
                   gAATCAGCGACTGAAAACGCAAACTCATTTGCAACGGCAGCATTTAGAAAACCAGATACACAATCATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGTCAGCAAATGCTGCCCAACACCAAcggcg                                                                                                                           <  2:898030/143‑1 (MQ=255)
                   gAATCAGCAACTGAAAACGCAAACTCAGTTGCAAGAGCAGCATTTAGAAAACCAGATAAACAATAATACTCAGCGGGTGCTGCAGTAACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCCGCCCAACACCAAcagcg                                                                                                                           >  2:895806/1‑143 (MQ=255)
                    aaTCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGCTGCAGTCACACCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGGTGCCCCACACCAGcggcag                                                                                                                          >  2:553501/1‑141 (MQ=255)
                         gCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGCTCTTGCAGTCACAGCCGGGGGCGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACCCCACCGGCGGGATgc                                                                                                                     >  1:917429/1‑142 (MQ=255)
                              tGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGACAGCCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAAc                                                                                                                <  2:607980/143‑1 (MQ=255)
                              tGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGACGGGATGGTAAAc                                                                                                                >  2:284338/1‑143 (MQ=255)
                                aaaaCGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAg                                                                                                              >  2:906897/1‑143 (MQ=255)
                                   aCGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCaa                                                                                                           >  2:909261/1‑143 (MQ=255)
                                        aaCTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGCTAAACAGCAACCGTa                                                                                                      >  2:87222/1‑143 (MQ=255)
                                          cTCATTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAAt                                                                                                    <  1:770746/143‑1 (MQ=255)
                                           tCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGCGGCGCGAAATCCCGCCCGTCAGCAAATGCTGCCCAACACCAACGACAGAAAAATAAACAGCAACCGTAATc                                                                                                   >  1:50227/1‑143 (MQ=255)
                                             aGTTGCAACAGCGGCAGTTAGAAAACCACCAAATCAAACATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCg                                                                                                 <  2:219716/143‑1 (MQ=255)
                                               ttGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGa                                                                                               >  1:746972/1‑143 (MQ=255)
                                                 gCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCACTCACAGCCGGGGGCGCGAAATCCCGCCCGGCAGCCAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGTTa                                                                                             >  1:557442/1‑143 (MQ=255)
                                                   atgagcagcATTTAGAACACCACATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCCGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGt                                                                                           <  2:533254/140‑1 (MQ=255)
                                                   aaaagcagcATTTAGAAAACCAGATAAACGATAATCCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGt                                                                                           <  1:282618/140‑1 (MQ=255)
                                                       gcagcaTTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGt                                                                                       <  1:102030/143‑1 (MQ=255)
                                                         agcaTTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGTGTATCCAGCCCAACAGCGGCGGCGGGATGTTAAACAGCAACCGTAATCGGGATAGTTCGtag                                                                                     >  1:396041/1‑141 (MQ=255)
                                                              ttAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATca                                                                                >  1:45464/1‑143 (MQ=255)
                                                                aGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATcagc                                                                              >  2:294135/1‑143 (MQ=255)
                                                                 gAACACCAGATAAACAATAATGCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATcagca                                                                             <  1:428962/143‑1 (MQ=255)
                                                                 gAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATcagca                                                                             <  1:838353/143‑1 (MQ=255)
                                                                 gAAAACCAGATAAACAACAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATcagca                                                                             <  2:945983/143‑1 (MQ=255)
                                                                        aGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCCACCGTAATCCGGGTAGTTCGTTGCATCCGCCGCatata                                                                      >  1:376323/1‑142 (MQ=255)
                                                                          agaacccataatTCTCAGCGGGTGCTGCAGTCACAGCCGGGGCAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATAtgct                                                                    <  1:96763/136‑1 (MQ=255)
                                                                            aaaCAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATAtgctgc                                                                  <  1:33271/143‑1 (MQ=255)
                                                                             caCAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGcc                                                                 <  1:379449/142‑1 (MQ=255)
                                                                             aaCAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGCTCGTTGAATCAGCAGCATATgccgcc                                                                 >  2:174737/1‑143 (MQ=255)
                                                                              gcAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCg                                                                <  2:601207/142‑1 (MQ=255)
                                                                              aCAATAATTCTCGGGTGGTGGTGCAGGCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATAtgctgctg                                                                <  1:1054910/143‑3 (MQ=255)
                                                                                    atTCTCAGCGGGGGTGGCAGTCACGGCCGGGGGACCTAAATGTGGCTCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCggagagg                                                          <  2:422553/143‑1 (MQ=255)
                                                                                     tgctcAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCggagagga                                                         <  1:351671/141‑1 (MQ=255)
                                                                                      tctcAGCGGGTGTTGCAGTCACAGCCGGGGGAGCAGAAGCCCTCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCggagaggag                                                        <  1:725069/143‑1 (MQ=255)
                                                                                      gctcAGCGGGTGGTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCggagaggag                                                        <  1:656990/142‑1 (MQ=255)
                                                                                        ccAGCGGGTGTGGCAGTCGCAGCCGGGGGTGCTAAAGTCCGCCCGGCAGCTAATGCTGCCCACCACCAACGGCGGGATGTTAACCCGCAACCGTAAACCGTATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGaa                                                      <  1:538342/142‑1 (MQ=255)
                                                                                                 gtTGCAGTCACAGCCGGGGCAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGCATAGTTCGTTGAATCAACAGCACATGCTCCCGGACAGGCGAAACAGCGACa                                             >  2:537775/1‑143 (MQ=255)
                                                                                                  ttGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACAt                                            >  1:398864/1‑143 (MQ=255)
                                                                                                         cacaGCCGTGGGAGGGGAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAAt                                     <  1:1058641/143‑1 (MQ=255)
                                                                                                         cacaGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCGGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAAt                                     <  2:248522/143‑1 (MQ=255)
                                                                                                                 ggggAGCGAAATCCCGCCCGGCAGAAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGGTAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGGAACGGCGACATGCTGAATCAGCCCAg                             >  1:293926/1‑143 (MQ=255)
                                                                                                                   ggAGCGAACGCCCGCCCGGCAGCCATCGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATCTGGTACCGGATAGGAGAAACGGCGACATGCTGAATCAGCCCAGCa                           <  2:1049554/143‑1 (MQ=255)
                                                                                                                         aaaTCCCGCCCGGCAGCAAAGGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCAcgccgc                     <  2:1067227/143‑1 (MQ=255)
                                                                                                                            tCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGc                  >  2:294357/1‑143 (MQ=255)
                                                                                                                               cgcccgGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGcccg               >  1:115211/1‑141 (MQ=255)
                                                                                                                               cgcccgGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTg               >  2:9229/1‑143 (MQ=255)
                                                                                                                                 ttgggCAGCATTTGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGCTAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCGAGAACGGCGCAGccgcac             >  2:922480/4‑139 (MQ=255)
                                                                                                                                 tgcgGCAGCAAATGCCGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGAc             <  1:29325/141‑1 (MQ=255)
                                                                                                                                  ccgGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACa            >  2:276313/1‑143 (MQ=255)
                                                                                                                                   gggCAGCATTTGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGCTAGTTCGTTGAATCAGCAGCACATGCTGCCGGCGAGGCCAAACGGCCACATGCTCAATCAGCCCAGAACGCCGAAGCCTGAccc           >  2:1015101/2‑141 (MQ=255)
                                                                                                                                   gggCAGCATTTGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGCGGAGAAACCGCGACATGCTGAATCAACCCCGCACGCCGCAGccgcgcat           >  2:1069438/2‑137 (MQ=255)
                                                                                                                                    ggCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACAtt          <  1:744679/143‑1 (MQ=255)
                                                                                                                                     gcagcaAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTc         <  1:745742/143‑1 (MQ=255)
                                                                                                                                         caAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGtt     >  1:1029748/1‑143 (MQ=255)
                                                                                                                                            aTGCTGCCCAACGCCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGaa  <  2:1010443/143‑1 (MQ=255)
                                                                                                                                            aTGCTGCCCAACACCAACGGCGGGAAGCTAAACAGCAACCGTAACCCGGATACTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGCAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGaa  <  1:895806/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                              
TCCAGCAAAAAGGGATGCTGAATCAGCAACTGAAAACGCAAACTCAGTTGCAACAGCAGCATTTAGAAAACCAGATAAACAATAATTCTCAGCGGGTGTTGCAGTCACAGCCGGGGGAGCGAAATCCCGCCCGGCAGCAAATGCTGCCCAACACCAACGGCGGGATGTTAAACAGCAACCGTAATCCGGATAGTTCGTTGAATCAGCAGCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAA  >  NC_000913/3587085‑3587367

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