Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,409,134 T→G 6.2% intergenic (+72/‑203) PP_1230 → / → nadA hypothetical protein/quinolinate synthase iron‑sulfur cluster subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,409,1340TG 6.2% 182.7 / 4.1 65intergenic (+72/‑203)PP_1230/nadAhypothetical protein/quinolinate synthase iron‑sulfur cluster subunit
Reads supporting (aligned to +/- strand):  ref base T (29/31);  new base G (2/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

AATGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTGGCGGTCAACATGCCGTCCGGCACCCTTAGCCGAT  >  NC_002947/1408985‑1409275
                                                                                                                                                     |                                                                                                                                             
aaTGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTCCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGtt                                                                                                                                               >  2:196430/1‑150 (MQ=255)
  tGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCa                                                                                                                                             >  2:58641/1‑150 (MQ=255)
  tGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCa                                                                                                                                             <  2:884756/150‑1 (MQ=255)
    gCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGACCCCGCCCCTACCAGGCTTCCACAGTGTAGGGGCGGGTCCAcc                                                                                                                                           >  2:451452/1‑150 (MQ=255)
    gCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCAcc                                                                                                                                           >  2:945683/1‑150 (MQ=255)
     ccGAAAAATGACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTGGGGGCGGGTTCAccc                                                                                                                                          >  2:258706/1‑150 (MQ=255)
        gaaaaTGACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCcgcg                                                                                                                                       >  1:996653/2‑150 (MQ=255)
          aaaTGACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTGAACCcgcccc                                                                                                                                     <  1:1239294/150‑4 (MQ=255)
             tGACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACCAACTAGGGGCGGGTTCACCCGCGAAggg                                                                                                                                  <  1:630916/150‑1 (MQ=255)
              gACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCAGCCCCTTCGCGGGTGAAACCGCCCATACGCGCATATTACAGTGTATGGGCGGGTTCACCCGGGAAgggg                                                                                                                                 <  2:431231/150‑1 (MQ=255)
               aCCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGc                                                                                                                                >  1:227369/1‑150 (MQ=255)
                ccGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGcc                                                                                                                               <  1:945683/150‑1 (MQ=255)
                    tGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGca                                                                                                                           <  2:227369/150‑1 (MQ=255)
                      aagaaTTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGaa                                                                                                                                     <  1:517151/138‑1 (MQ=255)
                      aagaaTTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGaa                                                                                                                                     >  2:517151/1‑138 (MQ=255)
                       agaaTTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACaa                                                                                                                        <  2:695544/150‑1 (MQ=255)
                           tttGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCAc                                                                                                                    >  1:75802/1‑150 (MQ=255)
                                     tCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAAcc                                                                                                          >  2:1588490/1‑150 (MQ=255)
                                          gcggcgCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTGAACCcgccccta                                                                                                                                   <  2:754511/120‑6 (MQ=38)
                                          gcggcgCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTGAACCcgccccta                                                                                                                                   >  1:754511/1‑115 (MQ=38)
                                             gcgcATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACaaaaa                                                                                                            <  2:976315/140‑1 (MQ=255)
                                             gcgcATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACaaaaa                                                                                                            >  1:976315/1‑140 (MQ=255)
                                               gcATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCAccc                                                                                                                                          <  2:1535071/108‑1 (MQ=255)
                                               gcATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCAccc                                                                                                                                          >  1:1535071/1‑108 (MQ=255)
                                                cATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAAt                                                                                               >  1:307639/1‑150 (MQ=255)
                                                cATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAAt                                                                                               <  1:227224/150‑1 (MQ=255)
                                                 aTCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCAcc                                                                                                                                           >  2:246617/1‑105 (MQ=255)
                                                 aTCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCAcc                                                                                                                                           <  1:246617/105‑1 (MQ=255)
                                                 aTCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCGTAGGGGGGGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCGCTACAAAAACGCCTTTAGCAATc                                                                                              >  2:227224/1‑150 (MQ=255)
                                                     gggggCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTg                                                                                          <  2:75802/150‑1 (MQ=255)
                                                      ggggCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTCCGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTgg                                                                                         <  1:488977/150‑1 (MQ=255)
                                                       gggCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGa                                                                                        <  2:1097042/150‑1 (MQ=255)
                                                            cTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTAc                                                                                                                 >  2:64612/1‑120 (MQ=255)
                                                            cTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTAc                                                                                                                 <  1:64612/120‑1 (MQ=255)
                                                            cTGTCGCCGCAGGGGTGAGGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCTCCCACGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTc                                                                                   <  1:451452/150‑1 (MQ=255)
                                                                        gggTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCgg                                                                       >  2:1388097/1‑150 (MQ=255)
                                                                             aTGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCGGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTcc                                                                  >  1:1332568/1‑150 (MQ=255)
                                                                              tGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTg                                                                                 <  2:1409784/134‑1 (MQ=255)
                                                                              tGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTg                                                                                 >  1:1409784/1‑134 (MQ=255)
                                                                              tGGTCTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATc                                                                                              >  1:976705/1‑121 (MQ=255)
                                                                              tGGTCTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATc                                                                                              <  2:976705/121‑1 (MQ=255)
                                                                                       cTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATc                                                        >  2:206032/1‑150 (MQ=255)
                                                                                        tGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATcc                                                       <  1:1303076/150‑1 (MQ=255)
                                                                                        tGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATcc                                                       <  2:996653/150‑1 (MQ=255)
                                                                                                 ctctTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCa                                              <  2:307639/150‑1 (MQ=255)
                                                                                                  tctTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCag                                             >  1:1210389/1‑150 (MQ=255)
                                                                                                  tctTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCag                                             <  1:58641/150‑1 (MQ=255)
                                                                                                  tctTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGAGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCag                                             >  2:1215697/1‑150 (MQ=255)
                                                                                                           gTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAg                                                                                                                                    <  2:744596/54‑1 (MQ=255)
                                                                                                           gTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAg                                                                                                                                    >  1:744596/1‑54 (MQ=255)
                                                                                                           gTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTa                                                                               >  2:894286/1‑107 (MQ=255)
                                                                                                           gTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTa                                                                               <  1:894286/107‑1 (MQ=255)
                                                                                                            tGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCAATACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTgg                                   <  2:627823/150‑1 (MQ=255)
                                                                                                              aaCCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTGGCg                                 <  1:196430/150‑1 (MQ=255)
                                                                                                                 ccGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTGGCGGTc                              >  1:1076135/1‑150 (MQ=255)
                                                                                                                 ccGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTGGCGGTc                              >  1:1241672/1‑150 (MQ=255)
                                                                                                                   gCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTGGCGGTCaa                            <  2:1241672/150‑1 (MQ=255)
                                                                                                                      ccTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTGGCGGTCAACAt                         <  1:1215697/150‑1 (MQ=255)
                                                                                                                          cGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTGGCGGTCAACATGCCg                     <  1:1588490/150‑1 (MQ=255)
                                                                                                                          cGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTGGCGGTCAACATGCCg                     >  2:932424/1‑150 (MQ=255)
                                                                                                                                tttCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTGGCGGTCAACATGCCGTCCGGc               >  2:202683/1‑150 (MQ=255)
                                                                                                                                       gtgtAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTGGCGGTCAACATGCCGTCCGGCACCCTTa        <  1:932424/150‑1 (MQ=255)
                                                                                                                                           cttcgCGGGTGAACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTGGCGGTCAACATGCCGTCCGGCACCCTTAGCCg    >  2:1353823/5‑150 (MQ=255)
                                                                                                                                             gggCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTGGCGGTCAACATGCCGTCCGGCACCCTTAGCCGAt  >  1:1257739/1‑150 (MQ=255)
                                                                                                                                               gCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTGGCg                                 >  1:327590/1‑117 (MQ=255)
                                                                                                                                               gCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTGGCg                                 <  2:327590/117‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
AATGGCCGAAAAATGACCGTTGAAGAATTTGCCATGTTCATGGCGGCGCATCTGGGGGCCCTGTCGCCGCAGGGGTGATGGTTTCTCCTGTTGCGGCCTCTTCGCGGGTGAACCCGCTCCTACGAGGATTTCACAGTGTAGGGGCGGGTTCACCCGCGAAGGGGCCCGCACAAGCACTACAAAAACCCCTTTAGTAATCGCTGGATTGTCTGTAGCCTTCGGCCTCCGATAATAATCCGCGCCCGCAGAGGGCCGTGGCGGTCAACATGCCGTCCGGCACCCTTAGCCGAT  >  NC_002947/1408985‑1409275

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: ‑