Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 818,694 C→T 5.8% A42V (GCG→GTG)  moaC → molybdopterin biosynthesis, protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913818,6940CT 5.8% 228.0 / 3.3 69A42V (GCG→GTG) moaCmolybdopterin biosynthesis, protein C
Reads supporting (aligned to +/- strand):  ref base C (28/37);  new base T (2/2);  total (30/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-02

CCATCCACATTTGAAGAAATAAGTATGTCGCAACTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGATCCCGCTCTGTCATCCGCTGATGCTCAGCAAAGTTGAAGTCAATTTACAGGCC  >  NC_000913/818546‑818839
                                                                                                                                                    |                                                                                                                                                 
ccatccaCATTTGAAGAAATAAGTATGTCGCAACTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGc                                                                                                                                                   <  3:453389/149‑1 (MQ=255)
     cacaTTTGAAGAAATAAGTATGTCGCAACTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCgatga                                                                                                                                              <  3:34859/149‑1 (MQ=255)
        aTTTGAAGAAATAAGTATGTCGCAACTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGatta                                                                                                                                           <  3:96393/149‑1 (MQ=255)
        aTTTGAAGAAATAAGTATGTCGCAACTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGatta                                                                                                                                           <  1:397599/149‑1 (MQ=255)
         tttGAAGAAATAAGTATGTCGCAACTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGattat                                                                                                                                          <  2:321254/149‑1 (MQ=255)
          ttGAAGAAATAAGTATGTCGCAACTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGattatt                                                                                                                                         <  2:274411/149‑1 (MQ=255)
            gatgaaATCAGTATGTCGCAACTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGCTGTCTCCGCCAAAGCGGAAACCGTGCGTCAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGa                                                                                                                                       <  4:349127/146‑1 (MQ=255)
             aagaaATAAGTATGTCGCAACTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGAt                                                                                                                                      <  2:22259/149‑1 (MQ=255)
               gaaATAAGTATGTCGCAACTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATgg                                                                                                                                    <  3:409555/149‑1 (MQ=255)
                 aaTAAGTATGTCGCAACTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTc                                                                                                                                  <  3:69695/149‑1 (MQ=255)
                 aaTAAGTATGTCGCAACTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTc                                                                                                                                  <  2:211858/149‑1 (MQ=255)
                        aTGTCGCAACTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGccacca                                                                                                                           <  1:33868/149‑1 (MQ=255)
                               aaCTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGc                                                                                                                    >  1:51634/1‑149 (MQ=255)
                                     cccATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTa                                                                                                              <  1:549937/149‑1 (MQ=255)
                                      ccATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTAt                                                                                                             >  2:437086/1‑149 (MQ=255)
                                              aCGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACt                                                                                                     >  4:262778/1‑149 (MQ=255)
                                                       gCGAAGCTCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTAtt                                                                                            >  2:314582/1‑149 (MQ=255)
                                                       gCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTAtt                                                                                            >  2:490099/1‑149 (MQ=255)
                                                       gCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTAtt                                                                                            >  1:513516/1‑149 (MQ=255)
                                                           aGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCg                                                                                        >  4:460013/1‑149 (MQ=255)
                                                           aGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCg                                                                                        >  2:51329/1‑149 (MQ=255)
                                                            gcgcACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCgg                                                                                       <  3:52620/149‑1 (MQ=255)
                                                                  aTGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCa                                                                                 >  2:452798/1‑149 (MQ=255)
                                                                   tggtggATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCAACACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAg                                                                                <  1:115881/149‑1 (MQ=255)
                                                                     gtggATGGCTCCGCCCAAGCGGAAACCGTGCGTGAAGGGCGGGCCGAAGCCTTTGTCACCATGCGCAGCGAGACGGTTGCGATTATTATTTATGGTCGCCAACCCAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAggc                                                                              <  1:187502/149‑1 (MQ=255)
                                                                            tctcCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAg                                                                                                                      <  1:303533/102‑1 (MQ=255)
                                                                            tctcCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAg                                                                                                                      >  2:303533/1‑102 (MQ=255)
                                                                            tctcCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCaaaa                                                                       >  2:215412/1‑149 (MQ=255)
                                                                             ctcCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAAc                                                                      <  1:443673/149‑1 (MQ=255)
                                                                              tcCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACGAAGGCGa                                                                                                                  >  4:519629/1‑104 (MQ=255)
                                                                              tcCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACGAAGGCGa                                                                                                                  <  3:519629/104‑1 (MQ=255)
                                                                              tcCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACg                                                                     >  2:267564/1‑149 (MQ=255)
                                                                                cGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGTGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGTCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCa                                                                   >  2:565086/1‑149 (MQ=255)
                                                                                cGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGTGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCa                                                                   >  4:216572/1‑149 (MQ=255)
                                                                                  ccAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCAcc                                                                 >  3:268415/1‑149 (MQ=255)
                                                                                   cAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCt                                                                >  2:440772/1‑149 (MQ=255)
                                                                                     aaGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTCTTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTgg                                                              >  3:466677/1‑149 (MQ=255)
                                                                                        cGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGAt                                                           <  4:356849/149‑1 (MQ=255)
                                                                                         ggAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATc                                                          >  4:271638/1‑149 (MQ=255)
                                                                                          gAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCt                                                         >  2:87029/1‑149 (MQ=255)
                                                                                             aCCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGAt                                                      >  1:243082/1‑149 (MQ=255)
                                                                                                gtgcgtgAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGATccc                                                   >  4:562498/1‑149 (MQ=255)
                                                                                                      gAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGATCCCGCTCTg                                             >  3:415918/1‑149 (MQ=255)
                                                                                                       aaGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGATCCCGCTCTGt                                            <  2:286084/149‑1 (MQ=255)
                                                                                                         gcgcGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGATCCCGCTCTGTCa                                          <  4:193590/149‑1 (MQ=255)
                                                                                                         gcgcGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGATCCCGCTCTGTCa                                          >  3:119985/1‑149 (MQ=255)
                                                                                                            cgggcggAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGATCCCGCTCTGTCATcc                                       <  3:311776/149‑1 (MQ=255)
                                                                                                                 ggAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGATCCCGCTCTGTCATCCGCTGa                                  <  1:539862/149‑1 (MQ=255)
                                                                                                                       cTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGATCCCGCTCTGTCATCCGCTGATGCTCa                            <  2:275551/149‑1 (MQ=255)
                                                                                                                              aCCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGATCCCGCTCTGTCATCCGCTGATGCTCAGCAAAGt                     >  2:45520/1‑149 (MQ=255)
                                                                                                                              aCCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGATCCCGCTCTGTCATCCGCTGATGCTCAGCAAAGt                     >  3:522158/1‑149 (MQ=255)
                                                                                                                              aCCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGATCCCGCTCTGTCATCCGCTGATGCTCAGCAAAGt                     >  3:462866/1‑149 (MQ=255)
                                                                                                                              aCCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGATCCCGCTCTGTCATCCGCTGATGCTCAGCAAAGt                     <  3:87089/149‑1 (MQ=255)
                                                                                                                               ccATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGAt                                                           >  3:504191/1‑110 (MQ=255)
                                                                                                                               ccATGCGCACCGAGACGCTGGCGATGATGATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGAt                                                           <  4:504191/110‑1 (MQ=255)
                                                                                                                                  tGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGATCCCGCTCTGTCATCCGCTGATGCTCTGCAAAGTTGaa                 <  3:453726/149‑1 (MQ=255)
                                                                                                                                  tGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGATCCCGCTCTGTCATCCGCTGATGCTCAGCAAAGTTGaa                 <  2:86473/149‑1 (MQ=255)
                                                                                                                                  tGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGATCCCGCTCTGTCATCCGCTGATGCTCAGCAAAGTTGaa                 <  3:94488/149‑1 (MQ=255)
                                                                                                                                  tGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGATCCCGCTCTGTCATCCGCTGATGCTCAGCAAAGTTGaa                 <  2:243082/149‑1 (MQ=255)
                                                                                                                                  tGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGATCCCGCTCTGTCATCCGCTGATGCTCAGCAAAGTTGaa                 <  3:347276/149‑1 (MQ=255)
                                                                                                                                  tGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGATCCCGCTCTGTCATCCGCTGATGCTCAGCAAAGTTGaa                 <  1:408592/149‑1 (MQ=255)
                                                                                                                                  tGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGATCCCGCTCTGTCATCCGCTGATGCTCAGCAAAGTTGaa                 <  3:208202/149‑1 (MQ=255)
                                                                                                                                       aGCGAGACGCTGGTGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGTCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGATCCCGCTCTGTCATCCGCTGATGCTCAGCAAAGTTGAAGTCaa            <  1:565086/149‑1 (MQ=255)
                                                                                                                                       aGCGAGACGCTGGTGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGATCCCGCTCTGTCATCCGCTGATGCTCAGCAAAGTTGAAGTCaa            <  3:216572/149‑1 (MQ=255)
                                                                                                                                        gCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGATCCCGCTCTGTCATCCGCTGATGCTCAGCAAAGTTGAAGTCAAt           <  1:215412/149‑1 (MQ=255)
                                                                                                                                         cGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGATCCCGCTCTGTCATCCGCTGATGCTCAGCAAAGTTGAAGTCAAtt          <  3:262778/149‑1 (MQ=255)
                                                                                                                                           agaCGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGATCCCGCTCTGTCATCCGCTGATGCTCAGCAAAGTTGAAGTCAATTTa        >  3:174059/1‑149 (MQ=255)
                                                                                                                                                cTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGATCCCGCTCTGTCATCCGCTGATGCTCAGCAAAGTTGAAGTCAATTTACAGGc   <  4:407767/149‑1 (MQ=255)
                                                                                                                                                cTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGATCCCGCTCTGTCATCCGCTGATGCTCAGCAAAGTTGAAGTCAATTTACAGGc   <  2:150134/149‑1 (MQ=255)
                                                                                                                                                 tGGCGATGATTATCGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGATCCCGCTCTGTCATCCGCTGATGCTCAGCAAAGTTGAAGTCAATTTACAGGcc  >  2:118903/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CCATCCACATTTGAAGAAATAAGTATGTCGCAACTGACCCATATCAACGCCGCTGGCGAAGCGCACATGGTGGATGTCTCCGCCAAAGCGGAAACCGTGCGTGAAGCGCGGGCGGAAGCCTTTGTCACCATGCGCAGCGAGACGCTGGCGATGATTATTGATGGTCGCCACCACAAAGGCGACGTATTTGCCACTGCGCGTATTGCCGGTATTCAGGCGGCAAAACGCACCTGGGATCTGATCCCGCTCTGTCATCCGCTGATGCTCAGCAAAGTTGAAGTCAATTTACAGGCC  >  NC_000913/818546‑818839

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: