Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 529,930 G→A 10.1% intergenic (+24/‑160) PP_5SC → / → birA 5S ribosomal RNA/bifunctional biotin‑[acetyl‑CoA‑carboxylase] ligase/biotin operon repressor BirA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947529,9300GA10.1% 170.0 / 10.0 69intergenic (+24/‑160)PP_5SC/birA5S ribosomal RNA/bifunctional biotin‑[acetyl‑CoA‑carboxylase] ligase/biotin operon repressor BirA
Reads supporting (aligned to +/- strand):  ref base G (24/38);  new base A (3/4);  total (27/42)
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.21e-01

TTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTATGCAAGCCAACAGCCCATGGCTATCATGCCTGCC  >  NC_002947/529807‑530068
                                                                                                                           |                                                                                                                                          
ttcgAACCACCTGATTCCATTCCGAACTCAGTAGTTAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTtgt                                                                                                                  <  2:1384778/147‑1 (MQ=2)
ttGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTtgt                                                                                                                  <  1:303051/150‑1 (MQ=11)
ttGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTtgt                                                                                                                  <  1:1344052/150‑1 (MQ=11)
ttGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTtgt                                                                                                                  <  2:403324/150‑1 (MQ=11)
ttGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTtgt                                                                                                                  >  1:1844535/1‑150 (MQ=11)
ttGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTtgt                                                                                                                  >  1:1838589/1‑150 (MQ=11)
ttGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCGATCTGCGCATGCAGGTAGGGGTTTTGTCTTtgt                                                                                                                  <  1:1587813/150‑1 (MQ=11)
  ggAACCACCTGATCCCAACCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCATGTAGCGGTTTTGTCTTtgtcg                                                                                                                <  2:563850/150‑3 (MQ=2)
    aaCCACCTGATCGCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGc                                                                                                              <  1:171088/150‑1 (MQ=17)
    aaCCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTgg                                                                                                                >  1:1282430/1‑148 (MQ=32)
    aaCCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTgg                                                                                                                <  2:1282430/148‑1 (MQ=32)
           tGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCGTaaaaa                                                                                                       >  1:293409/1‑150 (MQ=25)
                 catcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTgg                                                                                                 >  1:584051/1‑150 (MQ=255)
                  atcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATc                                                                                                     >  2:554247/1‑145 (MQ=40)
                  atcccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATc                                                                                                     <  1:554247/145‑1 (MQ=40)
                     ccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCtt                                                                                             <  2:1838589/150‑1 (MQ=255)
                     ccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCt                                                                                              >  2:1872330/1‑149 (MQ=255)
                     ccGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCt                                                                                              <  1:1872330/149‑1 (MQ=255)
                       gAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGTCCTTGc                                                                                           <  2:584051/150‑1 (MQ=255)
                       gAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGCAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGc                                                                                           <  1:1844725/150‑1 (MQ=255)
                        aaCTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAGGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGcc                                                                                          >  2:1255600/1‑150 (MQ=255)
                               tagtGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGctc                                                                                   <  2:583688/150‑1 (MQ=255)
                                  tGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGctcctc                                                                                >  1:926769/1‑150 (MQ=255)
                                    aaaCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTccc                                                                              >  1:1165247/1‑150 (MQ=255)
                                      aCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCtg                                                                            <  2:926769/150‑1 (MQ=255)
                                      aCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTCGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCtg                                                                            <  1:1569578/150‑1 (MQ=255)
                                        gATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCtgtg                                                                          >  2:1100643/1‑150 (MQ=255)
                                         aTGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGa                                                                         >  2:877837/1‑150 (MQ=255)
                                            cATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCtg                                                                            >  2:1152850/1‑144 (MQ=255)
                                            cATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCtg                                                                            <  1:1152850/144‑1 (MQ=255)
                                             aTCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTTGGTTTAAAAATCGTGGCCTTGCCAATGCTCCTCCCTGTGATTGc                                                                     <  2:1191221/150‑1 (MQ=255)
                                              tCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGc                                                                                           <  1:1756162/127‑1 (MQ=255)
                                              tCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGc                                                                                           >  2:1756162/1‑127 (MQ=255)
                                                   gATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGCGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATAcg                                                               <  2:1394463/150‑1 (MQ=255)
                                                    aTGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATAcgc                                                              <  2:1875292/150‑1 (MQ=255)
                                                     tGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATAcgcg                                                             <  2:1895263/150‑1 (MQ=255)
                                                     tGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATAcgcg                                                             <  2:1844535/150‑1 (MQ=255)
                                                           tgtgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCGAATGCGCCTCCCTGTTATTGCGATACGCGGCctct                                                       >  2:804121/1‑150 (MQ=255)
                                                                   ttCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAgg                                               <  2:929934/150‑1 (MQ=255)
                                                                   ttCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAgg                                               <  1:804121/150‑1 (MQ=255)
                                                                    tCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGc                                              <  2:1216094/150‑1 (MQ=255)
                                                                           gtgAGAGTGGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTc                                       <  2:1919154/150‑1 (MQ=255)
                                                                           gtgAGAGTAGGTCATCGTCAGGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTc                                       <  1:1255600/150‑1 (MQ=255)
                                                                              agagTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGa                                                                   <  1:263443/119‑1 (MQ=255)
                                                                              agagTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGa                                                                   >  2:263443/1‑119 (MQ=255)
                                                                               gagTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTgg                                                                                                                >  1:318946/1‑73 (MQ=34)
                                                                               gagTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTgg                                                                                                                <  2:318946/73‑1 (MQ=34)
                                                                               aaggccaggatttttaagcccacaaagaCAAAACCCCTACCTGCATGCGCAGATAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTAt                                   <  1:1751606‑M1/122‑1 (MQ=255)
                                                                                  gccaggatttttaagcccacaaagaCAAAACCCCTACCTGCATGCGCAGATAGGGGTTTTGCCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTATGCa                                <  2:1543469‑M1/125‑1 (MQ=255)
                                                                                   aGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCtg                                                                            >  2:1317999/1‑105 (MQ=255)
                                                                                   aGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCtg                                                                            <  1:1317999/105‑1 (MQ=255)
                                                                                    caggatttttaagcccacaaagaCAAAACCCCTACCTGCATGCGCAGATAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCa                                                 <  2:604638‑M1/108‑1 (MQ=255)
                                                                                    caggatttttaagcccacaaagaCAAAACCCCTACCTGCATGCGCAGATAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCa                                                 >  1:604638‑M1/24‑131 (MQ=255)
                                                                                     aggatttttaagcccacaaagaCAAAACCCCTACCTGCATGCGCAGATAGGGGTTTTGCCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTATGCAAGc                             >  1:1543469‑M1/23‑150 (MQ=255)
                                                                                        aTCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTATGCAAGCCaa                          <  2:101092/150‑1 (MQ=255)
                                                                                            tCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTATGCAAGCCAACAGc                      <  1:877837/150‑1 (MQ=255)
                                                                                                   cacaaagaCAAAACCCCTACCTGCATGCGCAGATAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTAt                                   >  1:792933‑M1/9‑130 (MQ=255)
                                                                                                   cacaaagaCAAAACCCCTACCTGCATGCGCAGATAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTAt                                   <  2:792933‑M1/122‑1 (MQ=255)
                                                                                                     aTTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTg                                                                                                  <  2:1792595/65‑1 (MQ=255)
                                                                                                     aTTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTg                                                                                                  >  1:1792595/1‑65 (MQ=255)
                                                                                                      tttCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCt                                                                                                             >  2:982512/1‑53 (MQ=255)
                                                                                                      tttCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCt                                                                                                             <  1:982512/53‑1 (MQ=255)
                                                                                                      tttCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTATGCAAGCCAACAGCCCATGGCTAt            >  2:1654948/1‑150 (MQ=255)
                                                                                                      tttCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTATGCAAGCCAACAGCCCATGGCTAt            >  2:894298/1‑150 (MQ=255)
                                                                                                      tttCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTATGCAAGCCAACAGCCCATGGCTAt            >  2:1134785/1‑150 (MQ=255)
                                                                                                         cGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTATGCAAGCCAACAGCCCATGGCTATCAt         <  2:1165247/150‑1 (MQ=255)
                                                                                                           CAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTtgtg                                                                                                                 <  2:426582‑M1/44‑1 (MQ=255)
                                                                                                           CAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTtgtg                                                                                                                 >  1:426582‑M1/1‑44 (MQ=255)
                                                                                                                ccacTATCTTCGCATGCAGGTAGGGGTATTGCCTTTGTGCGCTTGAGAATGCTGGCCGTGCCAGTACTGCTCCCTGTGATTGCGGTACGCGGCATGTTGTGCAGGCCCTTTTCCTCTGCAAGCCAACAGCCCATGGCTATCAtgcctgcc  <  2:293409/147‑1 (MQ=255)
                                                                                                                           |                                                                                                                                          
TTGGAACCACCTGATCCCATCCCGAACTCAGTAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCATTTCGCAAAACCCCTATCTGCGCATGCAGGTAGGGGTTTTGTCTTTGTGGGCTTAAAAATCCTGGCCTTGCCAATGCTCCTCCCTGTGATTGCGATACGCGGCCTCTCCTGCAGGCCCTTTTCCTATGCAAGCCAACAGCCCATGGCTATCATGCCTGCC  >  NC_002947/529807‑530068

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: