Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 182,659 C→T 7.1% A66V (GCG→GTG)  cdaR → DNA‑binding transcriptional activator CdaR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913182,6590CT 7.1% 213.3 / 4.1 70A66V (GCG→GTG) cdaRDNA‑binding transcriptional activator CdaR
Reads supporting (aligned to +/- strand):  ref base C (30/35);  new base T (3/2);  total (33/37)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTACCATGCGCATCATCGATACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAACTGGTCTGCATGACGGCTGAAATGATG  >  NC_000913/182511‑182807
                                                                                                                                                    |                                                                                                                                                    
cGTACCATGCGCATCATCGATACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGAcgcg                                                                                                                                                     >  2:45163/1‑150 (MQ=255)
 gTACCATGCGCATCATCGATACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCgg                                                                                                                                                    <  1:1415834/150‑1 (MQ=255)
      attcgcTTCATCGATACCAATATCATCGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTTGTACTTTCACAGGGACTAGTCGTCGATATCGATGACGCGGTAGCa                                                                                                                                               <  2:1270206/147‑1 (MQ=255)
         cgcATCATCGATACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGt                                                                                                                                            >  2:356108/1‑150 (MQ=255)
         cgcATCATCGATACCAATATCAAAGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTAATTTCACATGGACGAGTCGTCGATATCGATGACCCGGTAGCACGt                                                                                                                                            >  2:78407/1‑150 (MQ=255)
            atcatcGATACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGTGCTACCGCGTCAt                                                                                                                                         >  1:1430859/1‑150 (MQ=255)
             tcatcGATACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGTGCTACCGCGTCATc                                                                                                                                        <  2:216899/150‑1 (MQ=255)
               atcGATACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTTCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTg                                                                                                                                      >  1:1170035/1‑150 (MQ=255)
                tcGATACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGTGCTACCGCGTCATcgat                                                                                                                                     <  2:1430859/150‑4 (MQ=255)
                tcGATACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGc                                                                                                                                     >  2:173523/1‑150 (MQ=255)
                  gATACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCAc                                                                                                                                   >  1:1343396/1‑150 (MQ=255)
                    tACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACgg                                                                                                                                 <  2:1343396/150‑1 (MQ=255)
                      ccAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGgtg                                                                                                                               <  1:863569/150‑1 (MQ=255)
                       cAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGgtgt                                                                                                                              >  2:301713/1‑150 (MQ=255)
                        aaTATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGtgtg                                                                                                                             <  2:394299/150‑1 (MQ=255)
                                    aTGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGAtt                                                                                                                 >  1:480819/1‑150 (MQ=255)
                                    aTGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGAtt                                                                                                                 <  2:107481/150‑1 (MQ=255)
                                           cccGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTAc                                                                                                          >  1:9578/1‑150 (MQ=255)
                                           cccGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTAc                                                                                                          >  1:1171101/1‑150 (MQ=255)
                                            ccGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTAcc                                                                                                         >  1:1356346/1‑150 (MQ=255)
                                            ccGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTAcc                                                                                                         >  1:931186/1‑150 (MQ=255)
                                              gTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGt                                                                                                       <  2:196031/150‑1 (MQ=255)
                                                 ggCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTAc                                                                                                    <  1:173523/150‑1 (MQ=255)
                                                  gCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACg                                                                                                   >  2:190407/1‑150 (MQ=255)
                                                  gCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACg                                                                                                   >  1:190316/1‑150 (MQ=255)
                                                   cGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACgg                                                                                                  <  1:1183679/150‑1 (MQ=255)
                                                      aTTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTg                                                                                               <  2:61345/150‑1 (MQ=255)
                                                         aTCGGCAGCGCCGATCGAGAGCGTATTGGAGAATTGCACGAACGTGCATTGCTGGTACTTTCACAGCGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGCGGATTAATCTACCGTTACGGCTGGaa                                                                                            >  2:518114/1‑150 (MQ=255)
                                                           cGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAgg                                                                                          >  1:1243344/1‑150 (MQ=255)
                                                            ggCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGt                                                                                         >  1:764825/1‑150 (MQ=255)
                                                             gCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTg                                                                                        <  2:532201/150‑1 (MQ=255)
                                                             gCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTg                                                                                        >  1:537754/1‑150 (MQ=255)
                                                             gCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTg                                                                                        <  1:657395/150‑1 (MQ=255)
                                                                gcggcgATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGaaa                                                                                     >  2:880930/1‑150 (MQ=255)
                                                                 cggcgATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAAt                                                                                    <  2:1077479/150‑1 (MQ=255)
                                                                     gATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTc                                                                                <  1:356108/150‑1 (MQ=255)
                                                                      aTCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCg                                                                               <  1:351891/150‑1 (MQ=255)
                                                                         gTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCg                                                                            <  1:1335773/150‑1 (MQ=255)
                                                                            aGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTaa                                                                         >  1:658488/1‑150 (MQ=255)
                                                                            aGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTaa                                                                         <  1:880930/150‑1 (MQ=255)
                                                                               gTATTGGTGAATTGCACGAAGGTGGATTGCTGGTACTTTCACTGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTg                                                                      <  2:728091/150‑1 (MQ=255)
                                                                                tATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTgg                                                                     <  1:45163/150‑1 (MQ=255)
                                                                                tATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTgg                                                                     <  2:473820/150‑1 (MQ=255)
                                                                                    ggTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTg                                                                 >  2:1084909/1‑150 (MQ=255)
                                                                                       gAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACa                                                              <  2:190316/150‑1 (MQ=255)
                                                                                        aaTTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAg                                                             >  2:432722/1‑150 (MQ=255)
                                                                                              aCGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAAc                                                       >  1:1381276/1‑150 (MQ=255)
                                                                                               cGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAAcc                                                      >  1:250304/1‑150 (MQ=255)
                                                                                                         ttGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGTGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTg                                            <  2:1201804/150‑1 (MQ=255)
                                                                                                         ttGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTg                                            <  2:1243344/150‑1 (MQ=255)
                                                                                                          tGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGc                                           <  1:1179260/150‑1 (MQ=255)
                                                                                                          tGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATAAATATACCGTTACGCGTGGAAGGTGAAATTGGCGGCGTAATTGGCCGGACAGGTGAACCAGAGAATCTGc                                           >  2:154585/1‑150 (MQ=255)
                                                                                                              ggTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTaa                                       >  2:1084407/1‑150 (MQ=255)
                                                                                                                 aCTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAata                                    <  1:301713/150‑1 (MQ=255)
                                                                                                                  cTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCTTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGTTGAACCAGAGAATCTGCGTAAatat                                   <  1:1059700/150‑1 (MQ=255)
                                                                                                                   tttCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACTGTGTGCGGCAGGGGATTAATCTACCGTTACGTCTGTAAGGTGAAATTGTCGGCGTTATTTGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATg                                  <  1:337132/150‑1 (MQ=255)
                                                                                                                    ttCCCAGGGACGCGTCGGCGATAGCGATGACGCGGTAGCACGTCATCTGCACGGTGGGCGGCAGGGGATTAATCTACCGTTCCGGCTGGAAGGGCAAATTGGCGGCGTAATGGGCCGGACAGGTGACCCAGAGAATCTGCGTAAATAggg                                 <  2:1044381/150‑1 (MQ=255)
                                                                                                                    ttCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATgg                                 >  1:100392/1‑150 (MQ=255)
                                                                                                                     tCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGATAATCTGCGTAAATATGGc                                <  2:626607/150‑1 (MQ=255)
                                                                                                                      cacaGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCg                               >  2:1021957/1‑150 (MQ=255)
                                                                                                                       acaGGGACGAGTCGTCGATATCGATCACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGa                              >  2:520136/1‑150 (MQ=255)
                                                                                                                         aGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAAc                            <  2:1127894/150‑1 (MQ=255)
                                                                                                                         aGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAAc                            >  1:856694/1‑150 (MQ=255)
                                                                                                                                aGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAACTGGTCTg                     >  1:729644/1‑150 (MQ=255)
                                                                                                                                 gtcgtcGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGCGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAACTGGTCTGc                    <  2:430410/150‑1 (MQ=255)
                                                                                                                                 ggcgtcgATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCGTGACTGGTGTACCAGAGAATCTGCGTAAATATGGCGAACTGGTCTGc                    <  1:78407/148‑1 (MQ=255)
                                                                                                                                   cgtcgATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAACTGGTCTGCAt                  <  1:432722/150‑1 (MQ=255)
                                                                                                                                      cgATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAACTGGTCTGCATGAc               <  1:518114/150‑1 (MQ=255)
                                                                                                                                      ccgtgcagATGACGTGCTACCGCGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAACTGGTCTGCATGAc               >  2:136091/8‑150 (MQ=255)
                                                                                                                                        gtgcagATGACGTGCTACCGCGTCATCTGCACGGGGTGCGGCAGGGGATTAATCTACCGGTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAACTGGTCTGCATGACgg             >  1:661914/6‑150 (MQ=255)
                                                                                                                                               tGACGCGGTAGCACGTCATCTGCACGGCGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAACTGGTCTGCATGACGGCTGAAat      <  2:1170035/150‑1 (MQ=255)
                                                                                                                                                   gcgGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAACTGGTCTGCATGACGGCTGAAatgatg  <  2:658488/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGTACCATGCGCATCATCGATACCAATATCAACGTAATGGATGCCCGTGGGCGAATTATCGGCAGCGGCGATCGTGAGCGTATTGGTGAATTGCACGAAGGTGCATTGCTGGTACTTTCACAGGGACGAGTCGTCGATATCGATGACGCGGTAGCACGTCATCTGCACGGTGTGCGGCAGGGGATTAATCTACCGTTACGGCTGGAAGGTGAAATTGTCGGCGTAATTGGCCTGACAGGTGAACCAGAGAATCTGCGTAAATATGGCGAACTGGTCTGCATGACGGCTGAAATGATG  >  NC_000913/182511‑182807

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