Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,872,578 A→G 6.5% N180D (AAC→GAC)  dgcJ → putative diguanylate cyclase DgcJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,872,5780AG 6.5% 182.6 / 3.0 61N180D (AAC→GAC) dgcJputative diguanylate cyclase DgcJ
Reads supporting (aligned to +/- strand):  ref base A (29/28);  new base G (2/2);  total (31/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 = 7.94e-01

GATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTATTTACCGGGCGTACCGTATTGTCGAGCATTTACCAGGATCAGCGTACTAAACAGAACGTGATGAGTTTGCTGACGCCGGTATATGTCGC  >  NC_000913/1872431‑1872726
                                                                                                                                                   |                                                                                                                                                    
gATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAAc                                                                                                                                                    >  2:252453/1‑150 (MQ=255)
 aTCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACt                                                                                                                                                   <  2:794424/150‑1 (MQ=255)
    aGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGt                                                                                                                                                >  1:210550/1‑150 (MQ=255)
     ggCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGtc                                                                                                                                               >  1:227888/1‑150 (MQ=255)
         gTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGtccttc                                                                                                                                           >  2:504585/1‑150 (MQ=255)
          tGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGtccttcc                                                                                                                                          <  1:61266/150‑1 (MQ=255)
           ggAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGtccttcct                                                                                                                                         >  2:435993/1‑150 (MQ=255)
           ggAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGtccttcct                                                                                                                                         >  2:672699/1‑150 (MQ=255)
            gAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTg                                                                                                                                        >  1:282939/1‑150 (MQ=255)
             aaGGAAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGc                                                                                                                                       <  2:1419543/150‑1 (MQ=255)
               ggCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAg                                                                                                                                     <  1:1051471/150‑1 (MQ=255)
               ggCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAg                                                                                                                                     <  1:504585/150‑1 (MQ=255)
                       ccaccaGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAg                                                                                                                             >  1:1088613/1‑150 (MQ=255)
                        cacaaGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCATGCAGTCAGGAAGt                                                                                                                            >  2:1002753/1‑150 (MQ=255)
                              gatgatTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGc                                                                                                                      <  2:227888/150‑1 (MQ=255)
                                 gatTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCatc                                                                                                                   >  1:631648/1‑150 (MQ=255)
                                    tACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCatcgat                                                                                                                >  2:911469/1‑150 (MQ=255)
                                     aCGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCatcgatc                                                                                                               >  1:1349475/1‑150 (MQ=255)
                                     aCGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCatcgatc                                                                                                               <  1:790957/150‑1 (MQ=255)
                                     aCGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCatcgatc                                                                                                               <  2:1367647/150‑1 (MQ=255)
                                           cAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAg                                                                                                         <  2:631648/150‑1 (MQ=255)
                                           cAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAg                                                                                                         <  1:197506/150‑1 (MQ=255)
                                             ggTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCatcgatc                                                                                                               >  1:438417/1‑142 (MQ=255)
                                             ggTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCatcgatc                                                                                                               <  2:438417/142‑1 (MQ=255)
                                              gTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGata                                                                                                      <  2:1349475/150‑1 (MQ=255)
                                                cTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATAtt                                                                                                    >  2:976938/1‑150 (MQ=255)
                                                 ttGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTg                                                                                                   <  2:282939/150‑1 (MQ=255)
                                                 ttGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTg                                                                                                   <  1:435993/150‑1 (MQ=255)
                                                  tGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGa                                                                                                  <  2:137515/150‑1 (MQ=255)
                                                       tGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAgg                                                                                             <  1:976938/150‑1 (MQ=255)
                                                         gCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTa                                                                                           <  2:180285/150‑1 (MQ=255)
                                                               gatcaaTTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTATTTAcc                                                                                     <  2:1389791/146‑1 (MQ=255)
                                                                   aaTTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTATTTACCGGGc                                                                                 >  1:307522/1‑150 (MQ=255)
                                                                            attatCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTATTTACCGGGCGTACCGTAt                                                                        >  1:1035930/1‑150 (MQ=255)
                                                                            attatCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTATTTACCGGGCGTACCGTAt                                                                        >  2:1076525/1‑150 (MQ=255)
                                                                             ttatCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTATTTACCGGGCGTACCGTAtt                                                                       <  2:1035930/150‑1 (MQ=255)
                                                                              tatCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTATTTACCGGGCGTACCGTATTg                                                                      >  1:972699/1‑150 (MQ=255)
                                                                                tctcGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATGACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTATTTACCGGGCGTACCGTATTGTc                                                                    >  1:218752/1‑150 (MQ=255)
                                                                                tctcGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATGACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTATTTACCGGGCGTACCGTATTGTc                                                                    >  1:223957/1‑150 (MQ=255)
                                                                                tctcGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTATTTACCGGGCGTACCGTATTGTc                                                                    <  1:1002753/150‑1 (MQ=255)
                                                                                 ctcGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTATTTACCGGGCGTACCGTATTGTCg                                                                   <  2:1082400/150‑1 (MQ=255)
                                                                                   cGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTATTTACCGGGCGTACCGTATTGTCGAg                                                                 >  1:465656/1‑150 (MQ=255)
                                                                                        tcAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTATTTACCGGGCGTACCGTATTGTCGAGCAttt                                                            >  1:678020/1‑150 (MQ=255)
                                                                                           aTGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTATTTACCGGGCGTACCGTATTGTCGAGCATTTAcc                                                         >  2:1315560/1‑150 (MQ=255)
                                                                                                   tATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTATTTACCGGGCGTACCGTATTGTCGAGCATTTACCAGGATCAg                                                 <  2:210550/150‑1 (MQ=255)
                                                                                                             tttATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTATTTACCGGGCGTACCGTATTGTCGAGCATTTACCAGGATCAGCGTACTAAAc                                       >  2:566572/1‑150 (MQ=255)
                                                                                                                  gAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTATTTACCGGGCGTACCGTATTGTCGAGCATTTACCAGGATCAGCGTACTAAACAGAAc                                  >  2:956730/1‑150 (MQ=255)
                                                                                                                  gAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTATTTACCGGGCGTACCGTATTGTCGAGCATTTACCAGGATCAGCGTACTAAACAGAAc                                  <  1:126235/150‑1 (MQ=255)
                                                                                                                    gCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTATTTACCGGGCGTACCGTATTGTCGAGCATTTACCAGGATCAGCGTACTAAACAGAACGt                                >  2:1331678/1‑150 (MQ=255)
                                                                                                                     ccGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTATTTACCGGGCGTACCGTATTGTCGAGCATTTACCAGGATCAGCGTACTAAACAGAACGtg                               >  1:181437/1‑150 (MQ=255)
                                                                                                                        gTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTATTTACCGGGCGTACCGTATTGTCGAGCATTTACCAGGATCAGCGTACTAAACAGAACGtgatg                            <  2:209817/150‑1 (MQ=255)
                                                                                                                            aTGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTATTTACCGGGCGTACCGTATTGTCGAGCATTTACCAGGATCAGCGTACTAAACAGAACGTGATGAGtt                        <  1:672699/150‑1 (MQ=255)
                                                                                                                            aTGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTATTTACCGGGCGTACCGTATTGTCGAGCATTTACCAGGATCAGCGTACTAAACAGAACGTGATGAGtt                        <  2:1088613/150‑1 (MQ=255)
                                                                                                                             tGTTGAATACTTTGCGATGAATGACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTATTTACCGGGCGTACCGTATTGTCGAGCATTTACCAGGATCAGCGTACTAAACAGAACGTGATGAGttt                       <  2:223957/150‑1 (MQ=255)
                                                                                                                             tGTTGAATACTTTGCGATGAATGACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTATTTACCGGGCGTACCGTATTGTCGAGCATTTACCAGGATCAGCGTACTAAACAGAACGTGATGAGttt                       <  2:218752/150‑1 (MQ=255)
                                                                                                                                     aCTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTATTTACCGGGCGTACCGTATTGTCGAGCATTTACCAGGATCAGCGTACt                                           >  1:785745/1‑122 (MQ=255)
                                                                                                                                     aCTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTATTTACCGGGCGTACCGTATTGTCGAGCATTTACCAGGATCAGCGTACt                                           <  2:785745/122‑1 (MQ=255)
                                                                                                                                        ttGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTATTTACCGGGCGTACCGTATTGTCGAGCATTTACCAGGATCAGCGTACTAAACAGAACGTGATGAGTTTGCTGACGCCgg            <  1:129733/150‑1 (MQ=255)
                                                                                                                                          gCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTATTTACCGGGCGTAAAGAATTTAAGAGAATATAAAAAGATAATGTTAGTAAACAGAACTTCATGAGTTTGCGGACGCCGGta          >  2:31296/1‑150 (MQ=255)
                                                                                                                                           cGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTATTTACCGGGCGTACCGTATTGTCGAGCATTTACCAGGATCAGCGTACTAAACAGAACGTGATGAGTTTGCTGACGCCGGtat         >  1:769320/1‑150 (MQ=255)
                                                                                                                                                  tAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTATTTACCGGGCGTACCGTATTGTCGAGCATTTACCAGGATCAGCGTACTAAACAGAACGTGATGAGTTTGCTGACGCCGGTATATGTCGc  >  2:1035968/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GATCAGGCAGTGGAAGGCAACCGCCACCAGGATGATTACGGTTCAGGTCTTGGGATGGCCGAAGAGAAATTTCACTATTATCTCGATCTCAATGACCGCTATGTCTATTTTTATGAGCCGGTTAATGTTGAATACTTTGCGATGAATAACTGGTCCTTCCTGCAGTCAGGAAGTATTGGCATCGATCGCAAAGATATTGAAAAGGTATTTACCGGGCGTACCGTATTGTCGAGCATTTACCAGGATCAGCGTACTAAACAGAACGTGATGAGTTTGCTGACGCCGGTATATGTCGC  >  NC_000913/1872431‑1872726

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