Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,782,257 T→G 6.1% N51T (AAC→ACC)  iscS‑II ← cysteine desulfurase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,782,2570TG 6.1% 192.0 / 4.3 66N51T (AAC→ACC) iscS‑IIcysteine desulfurase
Reads supporting (aligned to +/- strand):  ref base T (27/35);  new base G (2/2);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

GACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATCGAAGTAGAGTGGGGCGCTAGGCATTGAGGGCAGGTCCGTGTCAGTGGCCGTGAAATCGGCGATCATTGACGTAAGCGTACC  >  NC_002947/2782110‑2782402
                                                                                                                                                   |                                                                                                                                                 
accTCGGGCGCGACCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTa                                                                                                                                                 >  1:30789/3‑150 (MQ=255)
   tCGGGCGCGACCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTAccc                                                                                                                                              <  2:30789/150‑1 (MQ=255)
      gTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCcagg                                                                                                                                           >  1:1419615/1‑150 (MQ=255)
       tCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCcaggc                                                                                                                                          >  1:1332196/1‑150 (MQ=255)
        cgcgCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGGTAACcaggcc                                                                                                                                         <  1:589073/150‑2 (MQ=255)
         gcgcCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCcaggcag                                                                                                                                        <  2:638853/150‑1 (MQ=255)
           gcCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCaggcaggc                                                                                                                                      >  1:875281/1‑150 (MQ=255)
             ccGAGGTCCATTCCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCGGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGAGGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGa                                                                                                                                    <  1:1361401/150‑1 (MQ=255)
             ccGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGa                                                                                                                                    >  2:967404/1‑150 (MQ=255)
             ccGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGa                                                                                                                                    >  1:894781/1‑150 (MQ=255)
             ccGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGa                                                                                                                                    >  2:86853/1‑150 (MQ=255)
             ccGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGGTAACcaggc                                                                                                                                          >  1:388870/1‑144 (MQ=255)
             ccGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGGTAACcaggc                                                                                                                                          <  2:388870/144‑1 (MQ=255)
                     cATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGtt                                                                                                                            >  1:739708/1‑150 (MQ=255)
                      aTACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGttt                                                                                                                           >  1:758633/1‑150 (MQ=255)
                         ccAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGa                                                                                                                        >  1:1379503/1‑150 (MQ=255)
                         ccAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACCTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGa                                                                                                                        >  1:1151002/1‑150 (MQ=255)
                             ctgctCGGCGTGCGCGCCGACGCACTCGGCGACCTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAAc                                                                                                                    <  2:1151002/150‑1 (MQ=255)
                                                aCGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTa                                                                                                 <  2:739708/150‑1 (MQ=255)
                                                   cACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTAc                                                                                                                                                >  1:447147/1‑100 (MQ=255)
                                                   cACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTAc                                                                                                                                                <  2:447147/100‑1 (MQ=255)
                                                   cACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTg                                                                                              <  1:302881/150‑1 (MQ=255)
                                                    aCTCGGCGACTTGCTGGCGCGCCTGCTCCACGGAGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTc                                                                                                                          <  2:1040619/121‑1 (MQ=255)
                                                    aCTCGGCGACTTGCTGGCGCGCCTGCTCCACGGAGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTc                                                                                                                          >  1:1040619/1‑121 (MQ=255)
                                                      tCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTggcg                                                                                           <  2:198186/150‑1 (MQ=255)
                                                            aCTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTaa                                                                                     <  2:1419615/150‑1 (MQ=255)
                                                             cTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAAt                                                                                    <  2:740609/150‑1 (MQ=255)
                                                             cTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAAt                                                                                    >  2:106231/1‑150 (MQ=255)
                                                              ttGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATc                                                                                   >  2:819633/1‑150 (MQ=255)
                                                                 cTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATCGaa                                                                                <  1:819633/150‑1 (MQ=255)
                                                                 cTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATCGaa                                                                                <  2:894781/150‑1 (MQ=255)
                                                                  tGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATCGAAg                                                                               <  2:1332196/150‑1 (MQ=255)
                                                                   ggCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATCGAAGt                                                                              <  1:86853/150‑1 (MQ=255)
                                                                    gcgcgcCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGTGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATCGAAGTa                                                                             <  1:726334/150‑1 (MQ=255)
                                                                         ccTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATCGAAGTAGAGTg                                                                        <  1:106231/150‑1 (MQ=255)
                                                                             cTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATCGAAGTAGAGTGGGgc                                                                    <  2:602119/150‑1 (MQ=255)
                                                                                 aCGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATCGAAGTAGAGTGGGGCGCTa                                                                <  1:1195083/150‑1 (MQ=255)
                                                                                 aCGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATCGAAGTAGAGTGGGGCGCTa                                                                >  2:1195083/1‑150 (MQ=255)
                                                                                    gCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATCGAAGTAGAGTGGGGCGCTAGGc                                                             <  1:40985/150‑1 (MQ=255)
                                                                                           gTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATCGAAGTAGAGTGGGGCGCTAGGAATTGAgg                                                      <  2:240325/150‑1 (MQ=255)
                                                                                            tGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTa                                                                                                 <  1:150977/106‑1 (MQ=255)
                                                                                            tGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTa                                                                                                 >  2:150977/1‑106 (MQ=255)
                                                                                                     gCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATCGAAGTAGATTGGGGCGCTAGGCATTGAGGGCAGGTCCgt                                            <  1:469224/150‑1 (MQ=255)
                                                                                                           ggCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATCGAAGTAGAGTGGGGCGCTAGGCATTGAGGGCAGGTCCGTGTCAGt                                      <  2:982351/150‑1 (MQ=255)
                                                                                                            gCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATCGAAGTAGAGTGGGGCGCTAGGCATTGAGGGCAGGTCCGTGTCAGTg                                     <  2:226037/150‑1 (MQ=255)
                                                                                                             cGTGGCCGCTGGACGCCGGGTTGCCGAAGTCGGCCTGGTTCAACGGGAAGGGGGCCATCGTGTCGGGAACGCGGTCGTCGACGGGGGTAGTGGCGGCGTCATCGACGTAGAGTGGGGCGCTAGGCCTTGTGTGCCTGTCCGTGTCattcg                                    <  2:1379503/150‑5 (MQ=255)
                                                                                                                    gCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATCGAAGTAGAGTGGGGCGCTAGGCATTGAGGGCAGGTCCGTGTCAGTGGCCGTGaa                             >  1:769309/1‑150 (MQ=255)
                                                                                                                           gCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAAcc                                                                                                                   <  1:625179/57‑1 (MQ=255)
                                                                                                                           gCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAAcc                                                                                                                   >  2:625179/1‑57 (MQ=255)
                                                                                                                           gCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACgg                                                                                                      <  1:377947/70‑1 (MQ=255)
                                                                                                                           gCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACgg                                                                                                      >  2:377947/1‑70 (MQ=255)
                                                                                                                           gCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATCGAAGTAGAGTGGGGCGCTAGGCATTGAGGGCAGGTCCGTGTCAGTGGCCGTGAAATCGGCg                      >  2:84877/1‑150 (MQ=255)
                                                                                                                                    ccGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATCGAAGTAGAGTGGGGCGCTAGGCATTGAGGGCAGGTCCGTGTCAGTGGCCGTGAAATCGGCGATc                   >  2:459662/1‑144 (MQ=255)
                                                                                                                                    ccGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATCGAAGTAGAGTGGGGCGCTAGGCATTGAGGGCAGGTCCGTGTCAGTGGCCGTGAAATCGGCGATc                   <  1:459662/144‑1 (MQ=255)
                                                                                                                                     cGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATCGAAGTAGAGTGGGGCGCTAGGCATTGAGGGCAGGTCCGTGTCAGTGGCCGTGAAATCGGCGATCATTGACg            >  1:846991/1‑150 (MQ=255)
                                                                                                                                         gTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATCGAAGTAGAGTGGGGCGCTAGGCATTGAGGGCAGGTCCGTGTCAGTGGCCGTGAAATCGGCGATCATTGACGTAAg        <  2:875281/150‑1 (MQ=255)
                                                                                                                                           ctgCCTGGGTAACCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATCGAAGTAGAGTGGGGCGCTAGGCATTGAGGGCAGGTCCGTGTCAGTGGCCGTGAAATCGGCGATCATTGACGTAAGCg      >  2:1018816/3‑150 (MQ=255)
                                                                                                                                             gCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATCGAAGTAGAGTGGGGCGCTAGGCATTGAGGGCAGGTCCGTGTCAGTGGCCGTGAAATCGGCGATCATTGACGTAAGCGTa    <  2:846991/150‑1 (MQ=255)
                                                                                                                                              ccTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATCGAAGTAGAGTGGGGCGCTAGGCATTGAGGGCAGGTCCGTGTCAGTGGCCGTGAAATCGGCGATCATTGACGTAAGCGTAc   <  1:1203401/150‑1 (MQ=255)
                                                                                                                                               cTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATCGAAGTAGAGTGGGGCGCTAGGCATTGAGGGCAGGTCCGTGTCAGTGGCCGTGAAATCGGCGATCATTGa              <  1:150926/138‑1 (MQ=255)
                                                                                                                                               cTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATCGAAGTAGAGTGGGGCGCTAGGCATTGAGGGCAGGTCCGTGTCAGTGGCCGTGAAATCGGCGATCATTGa              >  2:150926/1‑138 (MQ=255)
                                                                                                                                               cTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATCGAAGTAGAGTGGGGCGCTAGGCATTGAGGGCAGGTCCGTGTCAGTGGCCGTGAAATCGGCGATCATTGACGTAAGCGTAcc  >  1:1147183/1‑150 (MQ=255)
                                                                                                                                                tGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATCGAAGTAGAGTGGGGCGCTAGGCATTGa                                                        <  1:319148/95‑1 (MQ=255)
                                                                                                                                                tGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATCGAAGTAGAGTGGGGCGCTAGGCATTGa                                                        >  2:319148/1‑95 (MQ=255)
                                                                                                                                                tGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATCGAAGTAGAGTGGGGCGCTAGGCATTGAGGGCAGGTCCgt                                            >  1:471606/1‑107 (MQ=255)
                                                                                                                                                tGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATCGAAGTAGAGTGGGGCGCTAGGCATTGAGGGCAGGTCCgt                                            <  2:471606/107‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GACTCGGTCGCGCCCGAGGTCCATACCAGCTGCTCGGCGTGCGCGCCGACGCACTCGGCGACTTGCTGGCGCGCCTGCTCCACGGCGTGACGTGCAGCCTGGCCGAAGGCGTGGCCGCTGGACGCCGGGTTGCCGAAGTTGGCCTGGTTACCCAGGCAGGCGAGCATGGTTTCGATAACCCGGTCGTCGACGGGGGTAGTGGCGGCGTAATCGAAGTAGAGTGGGGCGCTAGGCATTGAGGGCAGGTCCGTGTCAGTGGCCGTGAAATCGGCGATCATTGACGTAAGCGTACC  >  NC_002947/2782110‑2782402

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: