Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 983,584 T→A 9.4% intergenic (+39/‑50) iscX → / → ndk cysteine desulfurase activity regulator/nucleoside diphosphate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947983,5840TA 9.4% 197.7 / 10.8 64intergenic (+39/‑50)iscX/ndkcysteine desulfurase activity regulator/nucleoside diphosphate kinase
Reads supporting (aligned to +/- strand):  ref base T (34/24);  new base A (3/3);  total (37/27)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.28e-01

ACCGCTGGGTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGAAGAAGCCG  >  NC_002947/983442‑983724
                                                                                                                                              |                                                                                                                                            
aCCGCTGGGTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGGCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTc                                                                                                                                       >  1:1594809/1‑150 (MQ=255)
    cTGGGTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCACTCCAGGCCGTCTCGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCgg                                                                                                                                   <  2:220587/150‑1 (MQ=255)
    cTGGGTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCgg                                                                                                                                   <  2:1381170/150‑1 (MQ=255)
       atctcTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAGCTGCGCGAGGCAATCCCGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGAATTATACGCGGGtt                                                                                                                                <  2:1207389/146‑1 (MQ=255)
            tGGCATTGCCAGAGTTCATCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGAATTATACGCGGGttcagat                                                                                                                           >  1:1207389/1‑145 (MQ=255)
                aTTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTc                                                                                                                       >  1:1132981/1‑150 (MQ=255)
                aTTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTc                                                                                                                       <  1:1490347/150‑1 (MQ=255)
                    ccAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTaa                                                                                                                             >  1:6534/1‑140 (MQ=255)
                    ccAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTaa                                                                                                                             <  2:6534/140‑1 (MQ=255)
                     cAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAac                                                                                                                  >  1:1522181/1‑150 (MQ=255)
                      agagTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTTGtttct                                                                                                                 >  2:712925/1‑146 (MQ=255)
                      agagTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGAATTATACGCGGGtt                                                                                                                                <  2:1633519/135‑1 (MQ=255)
                      agagTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGAATTATACGCGGGtt                                                                                                                                >  1:1633519/1‑135 (MQ=255)
                       gagTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAacac                                                                                                                <  2:1132981/150‑1 (MQ=255)
                                  cgaTCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTgg                                                                                                     >  1:109788/1‑150 (MQ=255)
                                   gaTCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGa                                                                                                    <  2:1522181/150‑1 (MQ=255)
                                                    tGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTgg                                                                                                     >  1:1437716/1‑132 (MQ=255)
                                                    tGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTgg                                                                                                     <  2:1437716/132‑1 (MQ=255)
                                                          aGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCggg                                                                                                                                  <  2:1337635/97‑1 (MQ=255)
                                                          aGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCggg                                                                                                                                  >  1:1337635/1‑97 (MQ=255)
                                                              cTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACttt                                                                         >  2:1531399/1‑150 (MQ=255)
                                                               tCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTtc                                                                        >  1:1437106/1‑150 (MQ=255)
                                                               tCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTtc                                                                        >  2:1068486/1‑150 (MQ=255)
                                                               tCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTtc                                                                        >  1:507708/1‑150 (MQ=255)
                                                                 gAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCAt                                                                                           >  2:182708/1‑129 (MQ=255)
                                                                 gAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCAt                                                                                           <  1:182708/129‑1 (MQ=255)
                                                                   gtcgaTCCAGGCGGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGCTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCca                                                                    >  2:714635/5‑150 (MQ=255)
                                                                    gCAATCCAGGCAGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCcat                                                                   >  2:1269092/1‑150 (MQ=255)
                                                                       aTCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTATGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTc                                                                                                                       <  2:1006157/95‑1 (MQ=255)
                                                                       aTCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTc                                                                                                                       >  1:1006157/1‑95 (MQ=255)
                                                                       aTCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCcatcat                                                                <  2:1566187/150‑1 (MQ=255)
                                                                            ggCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACCTACTTTCTCCATCATCAAGc                                                           >  2:605773/1‑150 (MQ=255)
                                                                                       tcaacgaAGCCGACTAAGCGTCACTTGCAGGTTAGTCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGtt                                                >  2:337642/4‑150 (MQ=255)
                                                                                       tCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGtt                                                <  2:1594809/150‑1 (MQ=255)
                                                                                       tCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGtt                                                >  2:88055/1‑150 (MQ=255)
                                                                                          acgaAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGcc                                             >  2:1398616/1‑150 (MQ=255)
                                                                                           cgaAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCa                                            >  2:529922/1‑150 (MQ=255)
                                                                                                ccGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAAc                                       <  2:109788/150‑1 (MQ=255)
                                                                                                        gCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCgg                               <  1:1269092/150‑1 (MQ=255)
                                                                                                           tCACTTGCAGGTTAGTCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCaa                            <  1:337642/150‑1 (MQ=255)
                                                                                                               ttGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAGCGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATc                        >  1:864319/1‑150 (MQ=255)
                                                                                                               ttGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATc                        >  2:1541792/1‑150 (MQ=255)
                                                                                                               ttGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATc                        >  1:600794/1‑150 (MQ=255)
                                                                                                                  cAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATcacc                     <  1:1531399/150‑1 (MQ=255)
                                                                                                                   aGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATcacca                    >  1:1035269/1‑150 (MQ=255)
                                                                                                                     gTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGc                                                                                        <  1:688371/80‑1 (MQ=255)
                                                                                                                     gTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGc                                                                                        >  2:688371/1‑80 (MQ=255)
                                                                                                                     gTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATcacccat                  >  1:274634/1‑147 (MQ=255)
                                                                                                                          gCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCtt             >  1:1562900/1‑150 (MQ=255)
                                                                                                                           cAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTc            <  2:1035269/150‑1 (MQ=255)
                                                                                                                           cAATCCATCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATCTTCTGAATGGTGTCACGGCTAATCAACGTACTTTCTCCATAATAAACACTGACGCTGTTGCTAAAAACGTCATCGGCAAGATCACCACTCGCTTc            <  1:529922/150‑1 (MQ=255)
                                                                                                                                  ccGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCAt                                  >  2:19231/1‑121 (MQ=255)
                                                                                                                                  ccGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCAt                                  <  1:19231/121‑1 (MQ=255)
                                                                                                                                  ccGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATcaccac                   >  1:637652/1‑136 (MQ=255)
                                                                                                                                  ccGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATcaccac                   <  2:637652/136‑1 (MQ=255)
                                                                                                                                  ccGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGaagaag     >  2:896156/1‑150 (MQ=255)
                                                                                                                                    ggAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGAAGAAGcc   >  1:94636/1‑150 (MQ=255)
                                                                                                                                     gAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGAAGAAGCCg  <  1:408272/150‑1 (MQ=255)
                                                                                                                                       aCCCGCGAATTATACGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTtct                                                                       >  2:297090/1‑79 (MQ=38)
                                                                                                                                       aCCCGCGAATTATACGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTtct                                                                       <  1:297090/79‑1 (MQ=38)
                                                                                                                                         ccGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATc                                 <  2:194475/115‑1 (MQ=255)
                                                                                                                                         ccGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACATACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATc                                 >  1:194475/1‑115 (MQ=255)
                                                                                                                                           gcgTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCaa                            <  2:1573517/118‑1 (MQ=255)
                                                                                                                                           gcgTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCaa                            >  1:1573517/1‑118 (MQ=255)
                                                                                                                                              |                                                                                                                                            
ACCGCTGGGTGCTGGCATTGCCAGAGTTCAGCGACGATCCGTCACGCGGCGGTGAGAAAGTGCTCGAGGCAATCCAGGCCGCCTGGATCGACGAAGCCGACTAAGCGTCACTTGCAGGTTAGGCAATCCCCCGGAACCCGCGTATAATTCGCGGGTTTAATTTTTCGTAACACTCATATTCTGGAGTTTTCCATGGCTGTTCAACGTACTTTCTCCATCATCAAGCCTGACGCTGTTGCCAAAAACGTCATCGGCAAGATCACCACTCGCTTCGAAGAAGCCG  >  NC_002947/983442‑983724

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