Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,407,287 G→A 5.2% H80H (CAC→CAT PP_3886 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,407,2870GA 5.2% 612.4 / 14.6 193H80H (CAC→CATPP_3886hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (97/86);  new base A (8/2);  total (105/88)
Fisher's exact test for biased strand distribution p-value = 1.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 195 total aligned reads displayed.

GATAGCCAGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGG  >  NC_002947/4407140‑4407436
                                                                                                                                                   |                                                                                                                                                     
     ccAGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAg                                                                                                                                                >  2:1427593/1‑150 (MQ=255)
       aGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGc                                                                                                                                              <  2:1017954/150‑1 (MQ=255)
       aGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAAGTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGc                                                                                                                                              >  2:244935/1‑150 (MQ=255)
          gAACTGCCTGGGAACAGGGCGACCCGGGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTg                                                                                                                                           <  2:3426353/150‑1 (MQ=255)
            aCTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCTCGGCACATGCCGGGGCttttt                                                                                                                                         <  1:4286260/150‑4 (MQ=255)
             cTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCTCGGCACATGCCGGGGCtttttc                                                                                                                                        <  1:3480149/150‑5 (MQ=255)
              tGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAAt                                                                                                                                       >  1:357530/1‑150 (MQ=255)
              tGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAAt                                                                                                                                       >  1:1827668/1‑150 (MQ=255)
                ccTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGACTTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATgg                                                                                                                                     >  2:4842869/1‑150 (MQ=255)
                ccTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATgg                                                                                                                                     <  1:4713371/150‑1 (MQ=255)
                ccTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATgg                                                                                                                                     <  1:1774536/150‑1 (MQ=255)
                ccTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATgg                                                                                                                                     <  2:1442602/150‑1 (MQ=255)
                ccTGGGAACAGGGCGACCCGAGGTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGTGCATCTAACGGCGCGACGCTATCCGTGGCTGGCTCCGGGGAACGAAAATGCCGCGGCCTGTGCCGTGGCTGGGAAggg                                                                                                                                     >  2:417837/1‑150 (MQ=255)
                 cTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATggg                                                                                                                                    <  1:2069094/150‑1 (MQ=255)
                 cTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATgg                                                                                                                                     >  1:244433/1‑149 (MQ=255)
                  tGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGt                                                                                                                                   <  2:5104438/150‑1 (MQ=255)
                   gggAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTg                                                                                                                                  <  2:2177343/150‑1 (MQ=255)
                   gggAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTg                                                                                                                                  <  2:2302380/150‑1 (MQ=255)
                       aCAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGgaga                                                                                                                              <  1:3678403/150‑1 (MQ=255)
                       aCAGGGCGACCCGAGGTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGgaga                                                                                                                              <  1:5130308/150‑1 (MQ=255)
                        cAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGAt                                                                                                                             >  1:1243635/1‑150 (MQ=255)
                         aGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATg                                                                                                                            >  1:2709387/1‑150 (MQ=255)
                          gggCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATgg                                                                                                                           >  2:3992104/1‑150 (MQ=255)
                              gACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAAcc                                                                                                                      <  2:2731556/150‑1 (MQ=255)
                               aCCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGATGCTTGGAATGGGTGGAGATGGCAAcc                                                                                                                      >  2:1012086/1‑150 (MQ=255)
                               aCCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAAcc                                                                                                                      <  2:4460252/150‑1 (MQ=255)
                                cccGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAAccc                                                                                                                     >  1:1722252/1‑150 (MQ=255)
                                 ccGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCt                                                                                                                    >  1:936110/1‑150 (MQ=255)
                                  cGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTg                                                                                                                   >  1:3987243/1‑150 (MQ=255)
                                     gTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCa                                                                                                                <  2:924506/150‑1 (MQ=255)
                                         cgccgATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTc                                                                                                            >  1:1634667/1‑150 (MQ=255)
                                         cgccgATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTc                                                                                                            <  1:499031/150‑1 (MQ=255)
                                          gccgATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCg                                                                                                           >  1:4169170/1‑150 (MQ=255)
                                          gccgATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCg                                                                                                           >  1:1873415/1‑150 (MQ=255)
                                             gATGACCCGTGACGGACTTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTg                                                                                                        >  2:2426118/1‑150 (MQ=255)
                                             gATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACTAAAAAGCCACGTCATGTGCCGAGGCTTGGAATGGGTGGAGATTGCAACCCTGTCAGGTCGCTg                                                                                                        >  1:1500750/1‑150 (MQ=255)
                                             gATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTg                                                                                                        <  1:2729456/150‑1 (MQ=255)
                                                  cccGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTg                                                                                                   >  1:5432584/1‑150 (MQ=255)
                                                     gTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGgcg                                                                                                >  2:3727733/1‑150 (MQ=255)
                                                         cGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTc                                                                                            >  2:4632900/1‑150 (MQ=255)
                                                         cGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTc                                                                                            >  2:4633091/1‑150 (MQ=255)
                                                           gAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTg                                                                                                                                  >  2:1825819/1‑110 (MQ=255)
                                                           gAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTccc                                                                                          >  1:373848/1‑150 (MQ=255)
                                                            aaTTGGGAGGGACGGGGCAGAAATCTGAGCATCTACCGTCTCGACGCTATCCATGACGGGCTCCATGGAACGAAAAAGCCCCGGCCTGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCt                                                                                         <  2:731625/150‑1 (MQ=255)
                                                             aTTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTc                                                                                        <  2:977191/150‑1 (MQ=255)
                                                             aTTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCAGGACTGGCTCCATTGAACGAACAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGACGCTGTCGTGGCGCTTCCCTc                                                                                        <  2:1089649/150‑1 (MQ=255)
                                                              tttGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTcc                                                                                       <  2:3119361/150‑1 (MQ=255)
                                                              tttGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTcc                                                                                       >  1:4123462/1‑150 (MQ=255)
                                                               ttGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCa                                                                                      <  1:1355162/150‑1 (MQ=255)
                                                               ttGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCa                                                                                      <  1:3382160/150‑1 (MQ=255)
                                                                tGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAg                                                                                     <  2:3568261/150‑1 (MQ=255)
                                                                    gggACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTccc                                                                                 <  2:4040044/150‑1 (MQ=255)
                                                                       aCTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCAc                                                                              >  1:3584181/1‑150 (MQ=255)
                                                                         tGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGc                                                                            >  2:2975619/1‑150 (MQ=255)
                                                                           ggCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTg                                                                          >  2:3303147/1‑150 (MQ=255)
                                                                              aGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGc                                                                                                                                                    >  2:5198248/1‑73 (MQ=255)
                                                                              aGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACt                                                                       >  1:662939/1‑150 (MQ=255)
                                                                                     tGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCAGTAAACGAAAATGACCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGAAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTAccc                                                                >  2:5396434/1‑150 (MQ=255)
                                                                                             tAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGt                                                        >  2:1114266/1‑150 (MQ=255)
                                                                                               aCGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTg                                                      >  2:4016160/1‑150 (MQ=255)
                                                                                               aCGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTg                                                      >  2:2355768/1‑150 (MQ=255)
                                                                                                      gACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCgg                                               >  1:693303/1‑150 (MQ=255)
                                                                                                      gACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCGTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCgg                                               <  1:2570506/150‑1 (MQ=255)
                                                                                                       tcGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCggg                                              <  1:4687917/149‑1 (MQ=255)
                                                                                                       aCGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCggg                                              <  2:2634867/150‑1 (MQ=255)
                                                                                                           tATCCATGACTTGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCttt                                          <  2:4048733/150‑1 (MQ=255)
                                                                                                            aTCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTg                                         <  2:1527414/150‑1 (MQ=255)
                                                                                                              ccATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGat                                       <  1:3818280/150‑1 (MQ=255)
                                                                                                               cATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATg                                                                                                                            <  1:530865/64‑1 (MQ=255)
                                                                                                               cATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGatc                                      <  1:3318552/150‑1 (MQ=255)
                                                                                                                aTGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGccgccg                                                <  1:1105389/139‑1 (MQ=255)
                                                                                                                aTGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGatca                                     >  1:2256126/1‑150 (MQ=255)
                                                                                                                     tGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTcc                                >  2:2817247/1‑150 (MQ=255)
                                                                                                                      ggCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTcccc                                                                                <  1:434814/101‑1 (MQ=255)
                                                                                                                      ggCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTcccc                                                                                >  2:434814/1‑101 (MQ=255)
                                                                                                                       gCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCaa                              >  1:3552159/1‑150 (MQ=255)
                                                                                                                        cTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCaaa                             >  2:4105577/1‑150 (MQ=255)
                                                                                                                           cATTTAACGAAAAAGCCCCGGCATGTGCCGTGGCTTGGAATGGGTTGAGATGGCAACCGTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTTTTTCACCTGCACGTTGCCGCCGGGCTTTGATCTTCTCAAAAAAc                          <  2:2565086/150‑1 (MQ=255)
                                                                                                                            aTTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTc                                                                                                      >  2:188692/1‑73 (MQ=255)
                                                                                                                            aTTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACg                         >  1:563652/1‑150 (MQ=255)
                                                                                                                                 aCGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACaaa                    >  1:1145823/1‑150 (MQ=255)
                                                                                                                                  cGAAAAAGCCCCGGCATGTGCCGAGGCTTGTAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTCCAAGTTGCCGCCTGGCTTTGATCATCTCCACACACGACaaac                   <  2:1311518/150‑2 (MQ=255)
                                                                                                                                  cGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAg                   <  1:1199402/150‑1 (MQ=255)
                                                                                                                                   gAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTcc                                                                                           <  2:815662/77‑1 (MQ=255)
                                                                                                                                   gAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGc                  >  2:3554952/1‑150 (MQ=255)
                                                                                                                                   gAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGc                  >  1:4303642/1‑150 (MQ=255)
                                                                                                                                   gAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGc                  >  1:1939220/1‑150 (MQ=255)
                                                                                                                                     aaaaGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGccc                >  1:2849820/1‑150 (MQ=255)
                                                                                                                                      caaGCCTCGGCACATGCCGGGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCa               >  2:122014/2‑150 (MQ=255)
                                                                                                                                      caaGCCTCGGCACATGCCGGGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCa               >  2:2157526/2‑150 (MQ=255)
                                                                                                                                      caaGCCTCGGCACATGCCGGGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCa               >  2:4892591/2‑150 (MQ=255)
                                                                                                                                      aaaGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCa               <  2:1722252/150‑1 (MQ=255)
                                                                                                                                       aaGCCTCGGCACATGCCGGGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAg              >  2:2982469/1‑150 (MQ=255)
                                                                                                                                       aaGCCTCGGCACATGCCGGGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAg              >  1:4842522/1‑150 (MQ=255)
                                                                                                                                         gCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCt            >  2:2154870/1‑150 (MQ=255)
                                                                                                                                           cccGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTgg          >  2:5295987/1‑150 (MQ=255)
                                                                                                                                            ccGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGatcatc                                   <  2:456082/124‑1 (MQ=255)
                                                                                                                                            ccGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGatcatc                                   >  1:456082/1‑124 (MQ=255)
                                                                                                                                             cGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGt        <  1:4105577/150‑1 (MQ=255)
                                                                                                                                                cATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGc     <  1:4417757/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GATAGCCAGCGAACTGCCTGGGAACAGGGCGACCCGAGTTCCGCCGATGACCCGTGACGGAATTTGGAGGGACTGGGCAGAAATCTGAGCATCTAACGTCTCGACGCTATCCATGACTGGCTCCATTTAACGAAAAAGCCCCGGCATGTGCCGAGGCTTGGAATGGGTGGAGATGGCAACCCTGTCAGGTCGCTGTCGTGGCGCTTCCCTCCAGTCCCCACGCTGACTGTTACCCCTGCACGTTGCCGCCGGGCTTTGATCATCTCCAAAAACGACAAAGCCCAGCTGGGTGGCCGG  >  NC_002947/4407140‑4407436

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