Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,741,995 G→A 7.5% intergenic (‑20/+393) PP_3305 ← / ← PP_3306 TerC family membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,741,9950GA 7.5% 216.2 / 4.4 77intergenic (‑20/+393)PP_3305/PP_3306TerC family membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (35/36);  new base A (3/3);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCCCGCGGTGGAGCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCAGTGAAAA  >  NC_002947/3741851‑3742144
                                                                                                                                                |                                                                                                                                                     
ttCCCGCGGTGGAGCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGtt                                                                                                                                                  <  2:2253654/150‑1 (MQ=255)
   ccGCGGTGGAGCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGaa                                                                                                                                               <  1:2095019/150‑1 (MQ=255)
    cgcgGTGGAGCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGaac                                                                                                                                              <  1:2137063/150‑1 (MQ=255)
    cgcgGTGGAGCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGaac                                                                                                                                              <  1:2137144/150‑1 (MQ=255)
       ggTGGAGCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGaacaac                                                                                                                                           <  2:692891/150‑1 (MQ=255)
        gTGGAGCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCGCCTACTGTACATTCAGCTTTGAACAACt                                                                                                                                          >  2:1879957/1‑150 (MQ=255)
            aGCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCATAGTTGTTCAACTACGa                                                                                                                                      <  2:1763780/150‑1 (MQ=255)
             gCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGAc                                                                                                                                     >  2:340181/1‑150 (MQ=255)
                 acCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTAt                                                                                                                                        <  2:208640/143‑1 (MQ=255)
                 acCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTAt                                                                                                                                        >  1:208640/1‑143 (MQ=255)
                 acCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCg                                                                                                                                 >  2:650894/1‑150 (MQ=255)
                 acCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCg                                                                                                                                 >  1:962960/1‑150 (MQ=255)
                   cTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGa                                                                                                                                      >  2:828572/1‑143 (MQ=255)
                   cTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGa                                                                                                                                      <  1:828572/143‑1 (MQ=255)
                   cTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGAc                                                                                                                               >  2:67087/1‑150 (MQ=255)
                   cTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGAc                                                                                                                               >  1:598397/1‑150 (MQ=255)
                        tCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAg                                                                                                                                                    <  1:1178658/124‑1 (MQ=255)
                        tCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAg                                                                                                                                                    >  2:1178658/1‑124 (MQ=255)
                                 aaCAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCATAGTTGTTCAACTACGACta                                                                                                                                   <  2:99280/132‑2 (MQ=255)
                                 aaCAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCATAGTTGTTCAACTACGACta                                                                                                                                   >  1:99280/1‑131 (MQ=255)
                                    aGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTc                                                                                                              >  1:429434/1‑150 (MQ=255)
                                        ttAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCagag                                                                                                          >  2:827494/1‑150 (MQ=255)
                                         tAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAgg                                                                                                         <  1:1102169/150‑1 (MQ=255)
                                             acaacaGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAg                                                                                                                                                    >  1:1271815/1‑103 (MQ=255)
                                             acaacaGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAg                                                                                                                                                    <  2:1271815/103‑1 (MQ=255)
                                                     gAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCa                                                                                             <  2:1569110/150‑1 (MQ=255)
                                                      aaaaaGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCAc                                                                                            <  2:598397/150‑1 (MQ=255)
                                                        aaaGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTAt                                                                                                                                        <  2:540684/104‑1 (MQ=255)
                                                        aaaGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTAt                                                                                                                                        >  1:540684/1‑104 (MQ=255)
                                                          aGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGt                                                                                                                                                   <  2:1413714/91‑1 (MQ=255)
                                                          aGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGt                                                                                                                                                   >  1:1413714/1‑91 (MQ=255)
                                                           ggAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGCATAACCGTCAGAGGGGCCCGGCGTCACGAATa                                                                                       <  2:1866725/150‑1 (MQ=255)
                                                            gAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGaacaac                                                                                                                                           >  1:1353854/1‑97 (MQ=255)
                                                            gAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGaacaac                                                                                                                                           <  2:1353854/97‑1 (MQ=255)
                                                             aaagccaaagccAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATaaa                                                                                     <  2:429434/150‑1 (MQ=255)
                                                                gccaaagccaGAGGTCAGTCCCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCa                                                                                                             <  1:1116987/123‑1 (MQ=255)
                                                                gccaaagccaGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCa                                                                                                             >  2:1116987/1‑123 (MQ=255)
                                                                gccaaagccaGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGaa                                                                                  >  2:184705/1‑150 (MQ=255)
                                                                 aacacgccaGAGGTCATTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAAt                                                                                 <  2:1357477/145‑1 (MQ=255)
                                                                  caaagccaGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCATAGTTGTTCAACTacga                                                                                                                                      >  1:2082691/1‑93 (MQ=255)
                                                                  caaagccaGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCATAGTTGTTCAACTACGa                                                                                                                                      <  2:2082691/96‑1 (MQ=255)
                                                                   aaagccaGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAAttt                                                                               >  2:460178/1‑150 (MQ=255)
                                                                       ccaGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCg                                                                           <  1:841230/150‑1 (MQ=255)
                                                                              tCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGa                                                                    <  1:460178/150‑1 (MQ=255)
                                                                                 gTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTa                                                                 >  1:410000/1‑150 (MQ=255)
                                                                                 ctACCAAGACGTTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCAc                                                                                            <  1:605911/122‑1 (MQ=255)
                                                                                 ctACCAAGACGTTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCAc                                                                                            >  2:605911/2‑123 (MQ=255)
                                                                                  tACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTcc                                                                                                                                  >  1:1462181/1‑84 (MQ=255)
                                                                                  tACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTcc                                                                                                                                  <  2:1462181/84‑1 (MQ=255)
                                                                                  tACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCgg                                                                                                                        >  1:2014315/1‑94 (MQ=255)
                                                                                  tACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCgg                                                                                                                        <  2:2014315/94‑1 (MQ=255)
                                                                                    ccAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGc                                                              <  2:542795/150‑1 (MQ=255)
                                                                                     cAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCTTTGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCt                                                             >  1:982345/1‑150 (MQ=255)
                                                                                         aaCGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCaa                                                         <  2:1441808/150‑1 (MQ=255)
                                                                                         aaCGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCaa                                                         <  1:67087/150‑1 (MQ=255)
                                                                                            gTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGGCTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTa                                                      <  1:1999893/150‑1 (MQ=255)
                                                                                            gTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTa                                                      >  2:2064911/1‑150 (MQ=255)
                                                                                            gTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTa                                                      >  1:1467118/1‑150 (MQ=255)
                                                                                               ttGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGaa                                                   <  2:410000/150‑1 (MQ=255)
                                                                                                tGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGaaa                                                  <  1:2064911/150‑1 (MQ=255)
                                                                                                 ggTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGaaaa                                                 <  1:340181/150‑1 (MQ=255)
                                                                                                  gTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAAc                                                >  2:677892/1‑150 (MQ=255)
                                                                                                  gTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAAc                                                >  1:1031979/1‑150 (MQ=255)
                                                                                                  gTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAAc                                                >  1:848095/1‑150 (MQ=255)
                                                                                                           cGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAGTGGAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTaacat                                       <  1:1879957/150‑5 (MQ=255)
                                                                                                               cTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAAt                                      >  2:150858/1‑147 (MQ=255)
                                                                                                               cTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAAt                                      <  1:150858/147‑1 (MQ=255)
                                                                                                                       aGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCAt                             <  1:1490853/148‑1 (MQ=255)
                                                                                                                       aGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCAt                             >  2:1490853/1‑148 (MQ=255)
                                                                                                                       aGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGa                           >  2:1489172/1‑150 (MQ=255)
                                                                                                                         aCATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCAt                             <  1:1504527/146‑1 (MQ=255)
                                                                                                                         aCATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCAt                             >  2:1504527/1‑146 (MQ=255)
                                                                                                                           aTGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAg                       >  1:2116386/1‑150 (MQ=255)
                                                                                                                               gCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTaa                   >  1:1132155/1‑150 (MQ=255)
                                                                                                                                     cctGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCt             <  1:827494/150‑1 (MQ=255)
                                                                                                                                      ctGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCtt            >  2:1744813/1‑150 (MQ=255)
                                                                                                                                           taGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGACTAAAGAATTTATCGAGCGCGAGAAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCAGt       >  1:1155361/1‑150 (MQ=255)
                                                                                                                                                gTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCAGTGaaaa  >  2:1630643/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
TTCCCGCGGTGGAGCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACGTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCAGTGAAAA  >  NC_002947/3741851‑3742144

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: