Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 980,758 G→C 6.9% S66T (AGT→ACT)  hscB → DnaJ‑like molecular chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947980,7580GC 6.9% 181.4 / 4.8 58S66T (AGT→ACT) hscBDnaJ‑like molecular chaperone
Reads supporting (aligned to +/- strand):  ref base G (34/20);  new base C (2/2);  total (36/22)
Fisher's exact test for biased strand distribution p-value = 6.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

GTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGT  >  NC_002947/980611‑980906
                                                                                                                                                   |                                                                                                                                                    
gTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGAGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTg                                                                                                                                                    >  1:475775/1‑150 (MQ=255)
     gATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGcgccg                                                                                                                                               >  2:574364/1‑150 (MQ=255)
       tctcGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCcgcg                                                                                                                                             >  2:580601/1‑150 (MQ=255)
         tcGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTc                                                                                                                                           >  1:818991/1‑150 (MQ=255)
         tcGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTc                                                                                                                                           <  1:645679/150‑1 (MQ=255)
         tcGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCGGCGCACTGCGCAGCGTc                                                                                                                                           >  1:589621/1‑150 (MQ=255)
         tcGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCGGCGCACTGCGCAGCGTc                                                                                                                                           <  2:589621/150‑1 (MQ=255)
           gACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCGGCGCACTGCGCAGCGTctg                                                                                                                                         <  2:292995/150‑3 (MQ=255)
            aCAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTcgcg                                                                                                                                        >  1:66113/1‑150 (MQ=255)
             cAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTcgcgc                                                                                                                                       >  2:336176/1‑150 (MQ=255)
              aaGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTcgcgcg                                                                                                                                      >  1:383583/1‑150 (MQ=255)
                 cTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCgcgcgcgc                                                                                                                                   >  1:456078/1‑150 (MQ=255)
                       aCTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACctg                                                                                                                             >  2:68613/1‑150 (MQ=255)
                           gCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGCCGCTTACCCGACGCTGCGCCGTGCGCCGCGTCGCGCGCGCTACCTGCTgg                                                                                                                         <  1:574364/150‑1 (MQ=255)
                           gCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTgg                                                                                                                         <  1:402353/150‑1 (MQ=255)
                             tATCGCGAGCTGGCCCTCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGcc                                                                                                                       >  2:336591/1‑150 (MQ=255)
                             tATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGcc                                                                                                                       >  2:959435/1‑150 (MQ=255)
                             tATCGCGAGCTGGCCCCCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGTGCGCTACCAGCTGGcc                                                                                                                       >  2:569554/1‑150 (MQ=255)
                              aTCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCa                                                                                                                      <  2:818991/150‑1 (MQ=255)
                                 gcgAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCa                                                                                                                   <  2:456078/150‑1 (MQ=255)
                                         gCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCAc                                                                                                           <  1:580601/150‑1 (MQ=255)
                                            cgcgAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGACCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGaa                                                                                                        <  1:336176/150‑1 (MQ=255)
                                              cgAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGt                                                                                                      <  1:56451/150‑1 (MQ=255)
                                                   tCCATCCTGATCGCTTTGCCTACGCCTCCGAGCGTGAGCAGCGTATTGCGATGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGcccc                                                                                                 >  2:80443/1‑150 (MQ=255)
                                                   tCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGcccc                                                                                                 >  1:666483/1‑150 (MQ=255)
                                                      aTCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCagg                                                                                              >  1:364064/1‑150 (MQ=255)
                                                      aTCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCagg                                                                                              >  1:105717/1‑150 (MQ=255)
                                                          tGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGt                                                                                          >  2:856781/1‑150 (MQ=255)
                                                            aTCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCa                                                                                        >  1:811961/1‑150 (MQ=255)
                                                             tCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCAc                                                                                       >  1:592383/1‑150 (MQ=255)
                                                              cGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACg                                                                                      <  1:517299/150‑1 (MQ=255)
                                                                 tttGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTc                                                                                   >  1:917911/1‑150 (MQ=255)
                                                                    gCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCAc                                                                                <  2:917911/150‑1 (MQ=255)
                                                                          gCCTCCGTGCGTGAGCTGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCt                                                                          <  2:570062/150‑1 (MQ=255)
                                                                                gAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTc                                                                    <  2:666483/150‑1 (MQ=255)
                                                                                    gTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCctgc                                                                <  2:364064/150‑1 (MQ=255)
                                                                                           gcGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTgcagca                                                         >  2:394239/1‑150 (MQ=255)
                                                                                            cGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGcancag                                                        <  2:696718/150‑1 (MQ=255)
                                                                                                   cgcTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGt                                                 <  2:475775/150‑1 (MQ=255)
                                                                                                      tGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGc                                              >  1:906783/1‑150 (MQ=255)
                                                                                                        gAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGTTCCGTACTATAGATCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGt                                            <  1:569554/150‑1 (MQ=255)
                                                                                                         aaaaaTCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTg                                           >  1:174021/1‑150 (MQ=255)
                                                                                                           aaaTCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGACACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGaa                                         >  1:516144/1‑150 (MQ=255)
                                                                                                                    gCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCgaa                                <  2:2472/150‑1 (MQ=255)
                                                                                                                          aaCGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTg                          >  2:833456/1‑150 (MQ=255)
                                                                                                                          aaCGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCGCGAAGACCTg                          >  1:116538/1‑150 (MQ=255)
                                                                                                                            cgacgCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCa                        >  1:456018/1‑150 (MQ=255)
                                                                                                                             gacgCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAg                       >  2:662180/1‑150 (MQ=255)
                                                                                                                              acgCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAgg                      <  2:66113/150‑1 (MQ=255)
                                                                                                                              acgCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAgg                      <  2:620474/150‑1 (MQ=255)
                                                                                                                                gCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGa                     <  2:139240/149‑1 (MQ=255)
                                                                                                                                gCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGa                     >  1:139240/1‑149 (MQ=255)
                                                                                                                                  ttACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGa                  >  1:556747/1‑150 (MQ=255)
                                                                                                                                      cAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGcc              >  2:6232/1‑150 (MQ=255)
                                                                                                                                      cAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGcc              >  1:813931/1‑150 (MQ=255)
                                                                                                                                      cAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGcc              >  2:903693/1‑150 (MQ=255)
                                                                                                                                        gACGCGGCGCACTGCGCAGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGa            >  2:897991/1‑150 (MQ=255)
                                                                                                                                                  aGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCCGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGgt  >  1:1064704/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GTCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGAAGTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTTGACGGT  >  NC_002947/980611‑980906

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