Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,433,665 C→T intergenic (‑105/‑25) smf ← / → def protein Smf/peptide deformylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,433,6650CT89.5% 84.8 / 3.5 38intergenic (‑105/‑25)smf/defprotein Smf/peptide deformylase
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (21/13);  total (23/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.73e-01

CGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAACTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAAC  >  NC_000913/3433421‑3433897
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                        
cgaaaTGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATCCGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGa                                                                                                                                                                                                                                    <  1:20440/248‑1 (MQ=255)
          gCCCTGTTTGCTGAAAACCAAACGCATCACTATACGACTGTCTCGCCTCCCAGTGCGCGATACCGACCATATACTACACGTACAACCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATAGACCATCCTTATCTCCCTGCCATAAGACGCCTTAGCCCTCTTTGCTATTGGTAAGTGATGCTGTCAATAACAGGGGGATGTGTCTAGAGTAGCAGAACTAATCTTTCTAATTCCTGacccaatctct                                                                                                                                                                                                                          >  2:21306/1‑241 (MQ=255)
             cTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTgg                                                                                                                                                                                                                        >  1:159939/1‑250 (MQ=255)
              tGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAACTCATGAACACATCTCTgg                                                                                                                                                                                                                        >  1:100574/1‑249 (MQ=255)
                    ctgcAATACAACCGCATCAATATGCGACTGTTTTGCCACCAAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGCGGGGGATTTGTCTCGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATATCTGGAGATTTAt                                                                                                                                                                                                                >  2:380996/1‑251 (MQ=255)
                        aaTACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGaa                                                                                                                                                                                                                                   >  1:178731/1‑228 (MQ=255)
                         aTACAACCGCAGCAATATGCGTCTGTTTTGCCACCCAGTGAGCGATACGGACCATATCACCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAg                                                                                                                                                                                                           <  1:320486/251‑1 (MQ=255)
                                     cAATATGCGACTGTTTTGACACCCAGTGAGCGATACGCACCATATCTGCGCCGTACAAGCTGCTGATATCCATAAATCGCAGCCAAATTTTAGTATCGACCATCCTTAACTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAgtgt                                                                                                                                                                                               <  2:32665/251‑1 (MQ=255)
                                       atatTCGACTGTTTTGCCACCCAGTGAGCGATGCGGACCATATCATCGCCGTACAAGCTGCTGTTACTCATTAAACGCAGCCAAATTTCTTTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTa                                                                                                                                                                                             <  1:5814/251‑1 (MQ=255)
                                             gACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATAtt                                                                                                                                                                                       >  1:301618/1‑251 (MQ=255)
                                                           ccAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCAAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCt                                                                                                                                                                         <  1:71958/251‑1 (MQ=255)
                                                                       aCGGGCCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCGCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTGATGTCAGTTTTGCAAGTGTTACATATGCCGGACGAGCGGCTTCGCAAAGTTg                                                                                                                                                              <  1:364583/250‑1 (MQ=255)
                                                                        cGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCt                                                                                                                                                            >  1:284287/1‑251 (MQ=255)
                                                                                      gCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTaga                                                                                                                                                 >  1:381932/1‑248 (MQ=255)
                                                                                              aGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATATTTGCGATTGGTCCGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGATAGAAGAAGTGAATg                                                                                                                                      >  1:352923/1‑251 (MQ=255)
                                                                                                               aaaCGCAGCCAAATTTATGTATCGACCATCCTTATCTCCCTGCAATAAGCAGCCTTAGCAATATTTGCGATTGGTCAGTGATGCTGTAAATCAGAGGGGGATTTGTCTCGAATAGAAGAAATATTCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCGAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAATTGCTCAACCGATAGAAGCAGTGAATGACGAAATTCAGCGtctc                                                                                                                     >  2:32752/1‑248 (MQ=255)
                                                                                                                         aaaTTTCTGTATCGACCATCCGTATCCCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGata                                                                                                           <  1:19328/251‑1 (MQ=255)
                                                                                                                          aaTTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAACTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGatat                                                                                                          <  1:30463/251‑1 (MQ=255)
                                                                                                                                               ttttccccTGTCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGTGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAATTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCa                                                                                       <  1:59340/244‑1 (MQ=255)
                                                                                                                                               tATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTAcac                                                                                        >  1:39863/1‑246 (MQ=255)
                                                                                                                                                        gCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGt                                                                              <  1:256954/249‑1 (MQ=255)
                                                                                                                                                                   ccTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTggcgg                                                                 >  1:82076/1‑251 (MQ=255)
                                                                                                                                                                    cTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCACTCAGAGGGGGATTTGTATAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGGATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTggcgg                                                                 >  2:270387/1‑250 (MQ=255)
                                                                                                                                                                                         tCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGCATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAACGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATc                                            >  1:40018/1‑250 (MQ=255)
                                                                                                                                                                                           aGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGCACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGACGTGAATGCAGAAATTCAGCGTATCGTCGCTGATATGTTCGAGACGATGTACGCAGAAGAAGGGATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAAc                                         >  1:171485/1‑251 (MQ=255)
                                                                                                                                                                                           aGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCCGGTTGATATCCATCAAc                                         >  1:199005/1‑251 (MQ=255)
                                                                                                                                                                                                 gCTGTAAATCAGAGGGCGATTTGTCTAGAATAGAAGAAATAATCTTTCTACTTCCTGAACACAACTCTGGAGATTTATTTCAGTTTTGCAAGTGTTACATTGACCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTTGAAGAAGTGAATGCATAAATTCAGCGTATCGTCGAGGATATGTTCGAGACGATGGACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACTTATc                                    <  2:284287/250‑1 (MQ=255)
                                                                                                                                                                                                      cAATCAGAGGGGGATTTGTCTAGAATATAAAAAATAATCTGTCTAATTTCTGAACACATCTCTGGAGATTGATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTGCGCTAAGTTGCTAAACCGGTATAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGTAGAATAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGt                               <  1:313121/250‑1 (MQ=255)
                                                                                                                                                                                                         tCAGAGGGGGATTTGTCTAGAATAGAAGAAATCATCTTTCTAATTCCGGAACACATCACTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGATTCGCAAAGTTGCTAAACCGGTAGAAGAAGAGAATGCAGACATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGGCTGGCGGCAACGCAGGTTGATATCCATCACCGTATCATTGTTAtt                           >  1:129219/1‑251 (MQ=255)
                                                                                                                                                                                                          cAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTg                          >  1:50203/1‑251 (MQ=255)
                                                                                                                                                                                                          cAGAGGGGGATGTGTCTAGAATAGAAGAAATAATCTTGCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCTGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTg                          <  1:66259/251‑1 (MQ=255)
                                                                                                                                                                                                             agGGGGATTTGTCTAGAATAGCAGAAATAATCTTTCTAATTCATGAACACATCTCTGGAGATATATGTAAGATTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAACGTTGCTAAACCGGTAGAAGAAGTGAATGCAGCAATTCCGAGTCTCGTCGATGATATGTTCGTGACGATGTACGCAGAAGACGGTATTGGCCGGGCGGCACCCCAGGGTGTTATCCATGAACGTATCATTGTTAttagg                        >  2:329418/1‑247 (MQ=255)
                                                                                                                                                                                                             agGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATGCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTg                          <  1:329418/248‑1 (MQ=255)
                                                                                                                                                                                                             agGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAATTCATGAACACATCTCTGGAGATTTAAGTAAGTTTTGCAAGTGTTACATATTCCGGCCGAGCGTCTTCGCAAAGTTGCTAAACCAGTAGAAGTAGTGAATGCAGAAATTCAGAGTATCGTCCATGAAATGTTCGAGCCGATGTAAGCAGAAGCAGGAACTGGCCTGACGGCAACCCAGGATTATATCCATCAACGTCTCATTGTTAttactg                       >  2:124458/1‑247 (MQ=255)
                                                                                                                                                                                                                           agaatagaaGAAATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGt         >  1:185415/1‑251 (MQ=255)
                                                                                                                                                                                                                             aatagaaGAAATAATCTTTCTAACTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAAt                                                                                                                                       <  2:125745/123‑1 (MQ=255)
                                                                                                                                                                                                                             aatagaaGAAATAATCTTTCTAACTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAAt                                                                                                                                       >  1:125745/1‑123 (MQ=255)
                                                                                                                                                                                                                                  aagaaATAATCTTTCTAATTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAAc  >  1:188812/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                        
CGCAATGTAAGCCCTGTTTGCTGCAATACAACCGCATCAATATGCGACTGTTTTGCCACCCAGTGAGCGATACGGACCATATCATCGCCGTACAAGCTGCTGATACTCATTAAACGCAGCCAAATTTCTGTATCGACCATCCTTATCTCCCTGCCATAAGCAGCCTTAGCAATCTTTGCGATTGGTCAGTGATGCTGTCAATCAGAGGGGGATTTGTCTAGAATAGAAGAAATAATCTTTCTAACTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAAC  >  NC_000913/3433421‑3433897

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