Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,900,343 A→C 5.6% intergenic (+259/‑239) PP_4308 → / → PP_4309 AsnC family transcriptional regulator/NCS1 family transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,900,3430AC 5.6% 195.3 / 2.2 71intergenic (+259/‑239)PP_4308/PP_4309AsnC family transcriptional regulator/NCS1 family transporter
Reads supporting (aligned to +/- strand):  ref base A (35/32);  new base C (2/2);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

TCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGATCAACTAATCAATCTTTCAAAAATTTAAACTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATACAACTTAAACAAAACATAGATACACAAGAATCCGCGACGCCGACCCCAATGTCAGCCGCCGCGCCCACA  >  NC_002947/4900200‑4900491
                                                                                                                                                |                                                                                                                                                    
tcttcAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGATCAACTAATCAATCTTTCAAAAATTTAAACTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTATCTTGTTAAAAGCTAAGtt                                                                                                                                                <  2:665716/150‑1 (MQ=255)
 cttcAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGATCAACTAATCAATCTTTCAAAAATTTAAACTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTc                                                                                                                                               <  2:1579133/150‑1 (MQ=255)
      gTGGCGCCCCAATTCGTTACACACTTTCAACCCTGTATTTTGACCAAATGATGAACTAATAAATCATTAAAAAATTTAAACATTTTCATACGGTGAACTTATATTGACTATAGTTCATTAACTTAAGTTTTAACAAGATAAGTTCGAAat                                                                                                                                          >  2:62309/1‑150 (MQ=255)
           gcCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGATCAACTAATCAATCTTTCAAAAATTTAAACTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTc                                                                                                                                               >  1:1033927/1‑140 (MQ=255)
           gcCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGATCAACTAATCAATCTTTCAAAAATTTAAACTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTc                                                                                                                                               <  2:1033927/140‑1 (MQ=255)
           gcCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGATCAACTAATCAATCTTTCAAAAATTTAAACTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCg                                                                                                                                     >  1:103375/1‑150 (MQ=255)
            ccacAAGTCGTTACACACTTTTCACCCTCTATTTTCAGCAAATGATCAACTAATCAATCTTTCAAAAATTTAAACTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGc                                                                                                                                    <  2:346479/147‑1 (MQ=255)
                aaTTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGATCAACTAATCAATCTTTCAAAAATTTAAACTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGt                                                                                                                                >  2:1545079/1‑150 (MQ=255)
                 aTTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGATCAACTAATCAATCTTTCAAAAATTTAAACTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGtt                                                                                                                               <  2:467761/150‑1 (MQ=255)
                   tCGTTACACACTTTCCACCCTCTATTTTGACCAAATGATCAACTAATCAATCTTTCAAAAATTTAAACTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGa                                                                                                                             >  2:1238857/1‑150 (MQ=255)
                        acacacTTTCCACCCTCTATTTTGACCAAATGATCAACTAATCAATCTTTCAAAAATTTAAACTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCt                                                                                                                        >  2:789426/1‑150 (MQ=255)
                         cacacTTTCCACCCTCTATTTTGACCAAATGATCAACTAATCAATCTTTCAAAAATTTAAACTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATAATCGCCTGTTGAAACCTg                                                                                                                       >  1:326601/1‑150 (MQ=255)
                            acTTTCCACCCTCTATTTTGACCAAATGATCAACTAATCAATCTTTCAAAAATTTAAACTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGctc                                                                                                                    >  2:506/1‑150 (MQ=255)
                                   aCCCTCTATTTTGACCAAATGATCAACTAATCAATCTTTCAAAAATTTAAACTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTATCTTGTTAAAAGCTAAGt                                                                                                                                                 <  1:202437/114‑1 (MQ=255)
                                   aCCCTCTATTTTGACCAAATGATCAACTAATCAATCTTTCAAAAATTTAAACTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTATCTTGTTAAAAGCTAAGt                                                                                                                                                 >  2:202437/1‑114 (MQ=255)
                                   aCCCTCTATTTTGACCAAATGATCAACTAATCAATCTTTCAAAAATTTAAACTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTAAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCtt                                                                                                             <  2:148269/150‑1 (MQ=255)
                                      ctctATTTTGACCAAATGATCAACTAATCAATCTTTCAAAAATTTAAACTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAGTCATCGCCTGTTGAAACCTGCTCTTCCCTTGGc                                                                                                          >  1:1257787/1‑150 (MQ=255)
                                            tttGACCAAATGATCAACTAATCAATCTTTCAAAAATTTAAACTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGaa                                                                                                    <  1:789426/150‑1 (MQ=255)
                                                aCCAAATGATCAACTAATCAATCTTTCAAAAATTTAAACTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCg                                                                                                >  2:369679/1‑150 (MQ=255)
                                                   aaaTGATCAACTAATCAATCTTTCAAAAATTTAAACTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTCGAACAGTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTATTCCCTTGGCATGGAACTCGCtt                                                                                             >  1:120666/1‑150 (MQ=255)
                                                     aTGATCAACTAATCAATCTTTCAAAAATTTAAACTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCtttt                                                                                           <  1:506/150‑1 (MQ=255)
                                                       gATCAACTTATCAATCTTTCAAAAATTTAATCTTGTTTATACGGACAACTTATATTGTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGCTCGCAGCGGTCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTtgt                                                                                         <  2:271589/150‑1 (MQ=255)
                                                       gATCAACTAATCAATCTTTCAAAAATTTAAACTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTtgt                                                                                         >  1:1502893/1‑150 (MQ=255)
                                                        aTCAACTAATCAATCTTTCAAAAATTTAAACTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTc                                                                                                                 >  2:223266/1‑125 (MQ=255)
                                                        aTCAACTAATCAATCTTTCAAAAATTTAAACTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTc                                                                                                                 <  1:223266/125‑1 (MQ=255)
                                                         tCAACTAATCAATCTTTCAAAAATTTAAACTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTCGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTtgtgt                                                                                       >  1:1476120/1‑150 (MQ=255)
                                                             ctcatcaataTTCCAGAAATTTAAACTTGGTCATACGGTCAACTTATATTCTCTATATTTCATTAGCTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTc                                                                                   >  2:1392217/4‑150 (MQ=255)
                                                             cTAATCAATCTTTCAAAAATTTAAACTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTc                                                                                   <  1:369679/150‑1 (MQ=255)
                                                             cTAATCAATCTTTCAAAAATTTAAACTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTc                                                                                   <  2:1242768/150‑1 (MQ=255)
                                                             cTAATCAATCTTTCAAAAATTTAAACTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTc                                                                                   <  1:1545079/150‑1 (MQ=255)
                                                                 tcaatcTTTCAAAAATTTAAACTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCgtgt                                                                               <  2:521420/150‑1 (MQ=255)
                                                                        ttCAAAAATTTAAACTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGaaa                                                                                                                           <  1:360829/99‑1 (MQ=255)
                                                                        ttCAAAAATTTAAACTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGaaa                                                                                                                           >  2:360829/1‑99 (MQ=255)
                                                                                   aaaCTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATACAACTTaaac                                                             <  2:1858387/150‑1 (MQ=255)
                                                                                       ttGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTGTACAACTTaaacaaaa                                                         <  2:1502893/150‑1 (MQ=255)
                                                                                       ttGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATACAACTTaaacaaaa                                                         >  1:1695641/1‑150 (MQ=255)
                                                                                          ttCATACGGTCAACTTATATTCTCTATATTTAATTAACTTAACTTTTAACAAGATAAGTACGAAATCATCCCCTGCTATCAGTTGATAGGCCTATAGCGTGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATACAGCCAAAACGAAACAt                                                      >  1:442955/1‑150 (MQ=255)
                                                                                             aTACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATACAACTTAAACAAAACATAGa                                                   >  1:1728269/1‑150 (MQ=255)
                                                                                              tACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATACAACTTAAACAAAACATAGa                                                   <  2:1291014/149‑1 (MQ=255)
                                                                                              tACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATACAACTTAAACAAAACATAGa                                                   >  1:1291014/1‑149 (MQ=255)
                                                                                                  gTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATACAACTTaaa                                                              >  1:4779/1‑134 (MQ=255)
                                                                                                  gTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATACAACTTaaa                                                              <  2:4779/134‑1 (MQ=255)
                                                                                                  gTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGCTCGTGTTTTGTATACAACTTAAACAAAACATAGATacac                                              >  2:1563406/1‑150 (MQ=255)
                                                                                                    cAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATACAACTTAAACAAAACATAGATACACaa                                            <  2:103375/150‑1 (MQ=255)
                                                                                                        ttATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATACAACtt                                                                 <  2:953698/125‑1 (MQ=255)
                                                                                                        ttATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATACAACtt                                                                 >  1:953698/1‑125 (MQ=255)
                                                                                                              tctcTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATACAACTTAAACAAAACATAGATACACAAGAATCCGCGa                                  <  2:1728269/150‑1 (MQ=255)
                                                                                                               ctctATATTTCATTACCTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATACAACTTAAACAAAATATAGATACACAAGAATCCGCGAc                                 >  2:173555/1‑150 (MQ=255)
                                                                                                                tctATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATACAACTTAAACAAACCATAGAt                                                  <  2:1879402/132‑1 (MQ=255)
                                                                                                                tctATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATACAACTTAAACAAAACATAGAt                                                  >  1:1879402/1‑132 (MQ=255)
                                                                                                                      tttCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATACAACTTAAACAAAACATAGATACACAAGAATCCGCGACGCCGAcc                          <  2:120666/150‑1 (MQ=255)
                                                                                                                       ttCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGCTCGTGTTTTGTATACAACTTAAACAAAACATAGATACACAAGAATCCGCGACGCCGAccc                         <  1:1563406/150‑1 (MQ=255)
                                                                                                                         ttcgaaCTTATCTTGTTAAAAGCTAAGTTCGAAATCATCGCCTGTTGACAAGTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATACAACTTAAACAAAACATAGATACACAAGAATCCGCGACGCCGACCCCa                       >  1:49676/5‑150 (MQ=255)
                                                                                                                          gttAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATACAACTTAAACAAAACATAGATACACAAGAATCCGCGACGCCGACCCCaa                      >  2:1327923/2‑150 (MQ=255)
                                                                                                                          aTTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATAc                                                                      <  1:1303305/102‑1 (MQ=255)
                                                                                                                          aTTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATAc                                                                      >  2:1303305/1‑102 (MQ=255)
                                                                                                                           ttAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGtttt                                                                            <  2:1196986/95‑1 (MQ=255)
                                                                                                                           ttAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGtttt                                                                            >  1:1196986/1‑95 (MQ=255)
                                                                                                                           ttAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATACAACTTAAACAAAACATAGATACACAAGAATCCGCGACGCCGACCCCAAt                     <  2:1498056/150‑1 (MQ=255)
                                                                                                                           ttAACTTAGCTTTTAACAACATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATACAACTTAAACAAAACATAGATACACAAGAATCCGCGACGCCGACCCCAAt                     <  1:1181089/150‑1 (MQ=255)
                                                                                                                                    cTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATACAACTTAAACAAAACATAGATACACAAGAATCCGCGACGCCGACCCCAAt                     >  2:497772/1‑141 (MQ=255)
                                                                                                                                    cTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATACAACTTAAACAAAACATAGATACACAAGAATCCGCGACGCCGACCCCAAt                     <  1:497772/141‑1 (MQ=255)
                                                                                                                                    cTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATACAACTTAAACAAAACATAGATACACAAGAATCCGCGACGCCGACCCCAATGTCAgccgc            >  2:690152/1‑150 (MQ=255)
                                                                                                                                     ttttAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTg                                                                                                            >  1:1696305/1‑53 (MQ=255)
                                                                                                                                     ttttAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTg                                                                                                            <  2:1696305/53‑1 (MQ=255)
                                                                                                                                     ttttAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATACAACTTAAACAAAACATAGATACACAAGAATCCGCGACGCCGACCCCAATGTCAgcagcc           >  1:517751/1‑150 (MQ=255)
                                                                                                                                          aaCAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATACAACTTAAACAAAACATAGATACACAAGAATCCGCGACGCCGACCCCAATGTCAGCCGCCgcgc       >  1:1417668/1‑150 (MQ=255)
                                                                                                                                          aaCAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATACAACTTAAACAAAACATAGATACACAAGAATCCGCGACGCCGACCCCAATGTCAGCCGCCgcgc       <  2:1695641/150‑1 (MQ=255)
                                                                                                                                           aCAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATACAAc                                                                   <  2:1046446/89‑1 (MQ=255)
                                                                                                                                           aCAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATACAAc                                                                   >  1:1046446/1‑89 (MQ=255)
                                                                                                                                               gATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATACAACTTAAACAAAACATAGATACACAAGAATCCGCGACGCCGACCCCAATGTCAGCCGCCGCGCCcaca  >  1:1912807/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGATCAACTAATCAATCTTTCAAAAATTTAAACTTGTTCATACGGTCAACTTATATTCTCTATATTTCATTAACTTAGCTTTTAACAAGATAAGTTCGAAATCATCGCCTGTTGAAACCTGCTCTTCCCTTGGCATGGAACTCGCTTTTGTGTGTTCGTGTTTTGTATACAACTTAAACAAAACATAGATACACAAGAATCCGCGACGCCGACCCCAATGTCAGCCGCCGCGCCCACA  >  NC_002947/4900200‑4900491

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