Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,514,850 T→A 6.2% L1159M (TTG→ATG)  PP_3108 → rhs family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,514,8500TA 6.2% 252.8 / 6.7 80L1159M (TTG→ATG) PP_3108rhs family protein
Reads supporting (aligned to +/- strand):  ref base T (35/40);  new base A (3/2);  total (38/42)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.83e-01

TGGGGAGGCGAATAGCCAAGCAGGCTGAAATCAATGGCGAGGTCGAGCAGAAGCGCTTCCTCTGGCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGAGCTGACGGATAGCGAGGGTGAGATTGTTTGGCAGGCGACGTATCGCTCTTGGGGCACGATAGAAACTTT  >  NC_002947/3514707‑3514998
                                                                                                                                               |                                                                                                                                                    
tGGGGAGGCGAATAGCCAAGCAGGCTGAAATCAATGGCGAGGTCGAGCAGAAGCGCTTCCTCTGGCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGgcgc                                                                                                                                                <  2:1554764/150‑1 (MQ=255)
  gggAGGCGAATAGCCAAGCAGGCTGAAATCAATGGCGAGGTCGAGCAGAAGCGCTTCCTCTGGCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCgg                                                                                                                                              <  1:1198118/150‑1 (MQ=255)
       gCGAATAGCCAAGCAGGCTGAAATCAATGGCGAGGTCGAGCAGAAGCCCTTCCTCTGGCAAGGGCTAAGAATGCTGCGGGACGAACCGCCCGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGa                                                                                                                                         >  2:684912/1‑150 (MQ=39)
          aaTAGCCAAGCAGGCTGAAATCAATGGCGAGGTCGAGCAGAAGCGCTTCCTCTGGCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCa                                                                                                                                      <  1:516875/150‑1 (MQ=255)
           aTAGCCAAGCAGGCTGAAATCAATGGCGAGGTCGAGCAGAAGCGCTTCCTCTGGCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCggg                                                                                                                                             >  1:1213755/1‑142 (MQ=255)
           aTAGCCAAGCAGGCTGAAATCAATGGCGAGGTCGAGCAGAAGCGCTTCCTCTGGCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCggg                                                                                                                                             <  2:1213755/142‑1 (MQ=255)
           aTAGCCAAGCAGGCTGAAATCAATGGCGAGGTCGAGCAGAAGCGCTTCCTCTGGCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAg                                                                                                                                     <  1:454336/150‑1 (MQ=255)
                 aaGCAGGCTGAAATCAATGGCGAGGTCGAGCAGAAGCGCTTCCTCTGGCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGaa                                                                                                                               >  2:815257/1‑150 (MQ=255)
                   gCAGGCTGAAATCAATGGCGAGGTCGAGCAGAAGCGCTTCCTCTGGCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAgg                                                                                                                             >  1:1078183/1‑150 (MQ=255)
                    cAGGCTGAAATCAATGGCGAGGTCGAGCAGAAGCGCTTCCTCTGGCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGga                                                                                                                            <  1:1778458/150‑1 (MQ=255)
                       gCTGAAATCAATGGCGAGGTCGAGCAGAAGCGCTTCCTCTGGCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGaa                                                                                                                         >  1:2040420/1‑150 (MQ=255)
                         tGAAATCAATGGCGAGGTCGAGCAGAAGCGCTTCCTCTGGCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAgg                                                                                                                       >  2:1609323/1‑150 (MQ=255)
                            aaTCAATGGCGAGGTCGAGCAGAAGCGCTTCCTCTGGCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCa                                                                                                                    <  2:1914602/150‑1 (MQ=255)
                              tCAATGGCGAGGTCGAGCAGAAGCGCTTCCTCTGGCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAg                                                                                                                        >  2:480362/1‑144 (MQ=255)
                              tCAATGGCGAGGTCGAGCAGAAGCGCTTCCTCTGGCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAg                                                                                                                        <  1:480362/144‑1 (MQ=255)
                                aaTGGCGAGGTCGAGCAGAAGCGCTTCCTCTGGCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAggg                                                                                                                      <  1:698527/144‑1 (MQ=255)
                                aaTGGCGAGGTCGAGCAGAAGCGCTTCCTCTGGCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAggg                                                                                                                      >  2:698527/1‑144 (MQ=255)
                                aaTGGCGAGGTCGAGCAGAAGCGCTTCCTCTGGCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAg                                                                                                                >  2:947756/1‑150 (MQ=255)
                                  tGGCGAGGTCGAGCAGAAGCGCTTCCTCTGGCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGt                                                                                                              >  2:718536/1‑150 (MQ=255)
                                   ggCGAGGTCGAGCAGAAGCGCTTCCTCTGTCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTc                                                                                                             >  1:63971/1‑150 (MQ=255)
                                     cGAGGTCGAGCAGAAGCGCTTCCTCTGGCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTa                                                                                                           <  1:815257/150‑1 (MQ=255)
                                      gAGGTCGAGCAGAAGCGCTTCCTCTGGCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTAt                                                                                                          <  2:1498024/150‑1 (MQ=255)
                                        ggTCGAGCAGAAGCGCTTCCTCTGGCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGAAGAAGGTCTATTa                                                                                                        >  2:12729/1‑150 (MQ=255)
                                            gAGCAGAAGCGCTTCCTCTGGCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTc                                                                                                    <  1:1417885/150‑1 (MQ=255)
                                                 gAAGCGCTTCCTCTGGCCAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCACAGCAGGAGCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATAc                                                                                               >  1:263211/1‑150 (MQ=37)
                                                 gAAGCGCTTCCTCTGGCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATAc                                                                                               <  2:1219028/150‑1 (MQ=255)
                                                  aaGCGCTTCCTCTGGCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATAcc                                                                                              <  1:718536/150‑1 (MQ=255)
                                                      gcTTCCTCTGGCAAGGGCTAAGGATGCTGCGGGAGAAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATc                                                                                          >  1:2267793/1‑150 (MQ=255)
                                                          ccTCTGGCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGAt                                                                                      <  1:1369635/150‑1 (MQ=255)
                                                            tctGGCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAggg                                                                                                                      <  2:312553/116‑1 (MQ=255)
                                                            tctGGCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAggg                                                                                                                      >  1:312553/1‑116 (MQ=255)
                                                             ctGGCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTgg                                                                                   <  2:1617560/150‑1 (MQ=255)
                                                             ctGGCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTgg                                                                                   <  1:947756/150‑1 (MQ=255)
                                                              tGGCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTggg                                                                                  >  2:1366240/1‑150 (MQ=255)
                                                                gCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCgg                                                                                                                                              <  2:1460371/88‑1 (MQ=255)
                                                                gCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCgg                                                                                                                                              >  1:1460371/1‑88 (MQ=255)
                                                                       cTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATGGGCGCGATTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCCATTACTTCCATACCGATCAGATTGGGACGCCGCtt                                                                         >  1:1250436/1‑150 (MQ=255)
                                                                         aaGGATGCTGCGGGAGGAAACGGCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGa                                                                       <  2:44238/150‑1 (MQ=255)
                                                                            gATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGAGCt                                                                    >  2:2050338/1‑150 (MQ=255)
                                                                                  gcGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGAGCTGACGGa                                                              <  2:1078183/150‑1 (MQ=255)
                                                                                   cGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTcc                                                                                                   >  1:526622/1‑112 (MQ=255)
                                                                                   cGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTcc                                                                                                   <  2:526622/112‑1 (MQ=255)
                                                                                              gCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCtt                                                                         <  1:1486557/127‑1 (MQ=255)
                                                                                              gCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCtt                                                                         >  2:1486557/1‑127 (MQ=255)
                                                                                              gCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGAGCTGACGGATAGCGAGGGTga                                                  >  1:366790/1‑150 (MQ=255)
                                                                                              gCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGAGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGAGCTGACGGATAGCGAGGGTga                                                  >  2:1532842/1‑150 (MQ=255)
                                                                                                  gggCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGAGCTGACGGATAGCGAGGGTGAGAtt                                              >  1:709190/1‑150 (MQ=255)
                                                                                                    gCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGAGCTGACGGATAGCGAGGGTGAGATTGt                                            >  2:1519353/1‑150 (MQ=255)
                                                                                                     cAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGAGCTGACGGATAGCGAGGGTGAGATTGtt                                           <  1:1609323/150‑1 (MQ=255)
                                                                                                      agagCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGAGCTGACGGATAGCGAGGGTGAGATTGttt                                          <  1:1519353/150‑1 (MQ=255)
                                                                                                      agagCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGAGCTGACGGATAGCGAGGGTGAGATTGttt                                          <  2:2294752/150‑1 (MQ=255)
                                                                                                       gagCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGAGCTGACGGATAGCGAGGGTGAGATTGTTTg                                         <  2:2267793/150‑1 (MQ=255)
                                                                                                        agCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGAGCTGACGGATAGCGAGGGTGAGATTGTTTgg                                        >  1:1962954/1‑150 (MQ=255)
                                                                                                             cTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGGGTGGGTTGTTGACAGTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGAGCTGACGGATAGCGAGGGTGAGATTGTTTGGCAGGc                                   <  2:63971/150‑1 (MQ=37)
                                                                                                             cTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGAGCTGACGGATAGCGAGGGTGAGATTGTTTGGCAGGc                                   <  1:684912/150‑1 (MQ=255)
                                                                                                                    cTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGt                                                                                                                                   >  2:2108510/1‑47 (MQ=255)
                                                                                                                    cTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGt                                                                                                                                   <  1:2108510/47‑1 (MQ=255)
                                                                                                                    cTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGAGCTGACGGATAGCGAGGGTGAGATTGTTTGGCAGGCGACatat                            <  1:1366240/150‑1 (MQ=255)
                                                                                                                     tGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGa                                                                       <  2:2254097/106‑1 (MQ=255)
                                                                                                                     tGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGa                                                                       >  1:2254097/1‑106 (MQ=255)
                                                                                                                     tGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGAGCTGACGGATAGCGAGGGTGAGATTGTTTGGCAGGCGACGTATc                           >  1:2144427/1‑150 (MQ=255)
                                                                                                                          gAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAg                                                                                        >  1:1033036/1‑84 (MQ=255)
                                                                                                                          gAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAg                                                                                        <  2:1033036/84‑1 (MQ=255)
                                                                                                                           aGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAg                                                                                                                >  2:2155246/1‑59 (MQ=255)
                                                                                                                           aGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAg                                                                                                                <  1:2155246/59‑1 (MQ=255)
                                                                                                                               gggTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGAGCTGACGGATAGCGAGGGTGAGATTGTTTGGCAGGCGACGTATCGCTCTTgggg                 <  1:2050338/150‑1 (MQ=255)
                                                                                                                               gggTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGAGCTGACGGATAGCGAGGGTGAGATTGTTTGGCAGGCGACGTATCGCTCTTgggg                 <  1:1532842/150‑1 (MQ=255)
                                                                                                                                ggTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGAGCTGACGGATAGCGAGGGTGAGATTGTTTGGCAGGCGACGTATCGCTCTTGGGGc                >  1:2296267/1‑150 (MQ=255)
                                                                                                                                 gTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGAGCTGACGGATAGCGAGGGTGAGATTGTTTGGCAGGCGACGTATCGCTCTTGGGGCa               >  1:2175910/1‑150 (MQ=255)
                                                                                                                                    gCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGAGCTGACGGATAGCGAGGGTGAGATTGTTTGGCAGGCGACGTATCGCTCTTGGGGCACGa            >  1:1337021/1‑150 (MQ=255)
                                                                                                                                     acccgcgcCAATGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGAGCTGACg                                                                <  1:183491/93‑1 (MQ=255)
                                                                                                                                     acccgcgcCAATGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGAGCTGACg                                                                >  2:183491/5‑97 (MQ=255)
                                                                                                                                     acccgcgcCAATGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGAGCTGACGGATAGCGAGGGTGAGATTGttt                                          >  1:2303524/5‑119 (MQ=255)
                                                                                                                                     acccgcgcCAATGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGAGCTGACGGATAGCGAGGGTGAGATTGttt                                          <  2:2303524/115‑1 (MQ=255)
                                                                                                                                      cccgcgcCAATGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGAGCTGACGGATAGCGAGGGTGAGATTGTTTGGCAGGCGACGTATCGCTCTTGGGGCACGATa          >  2:1310263/4‑150 (MQ=255)
                                                                                                                                        tGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGAGCTGACGGATAGCGAGGGTGAGATTGTTTGGCAGGCGACGTATCGCTCTTGGGGCACGATAGa        >  1:595666/1‑150 (MQ=255)
                                                                                                                                          cgcCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGAGCTGACGGATAGCGAGGGTGAGATTGTTTGGCAGGCGACGTATCGCTCTTGGGGCACGATAGaaa      <  2:2040420/150‑1 (MQ=255)
                                                                                                                                           gcCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGAGCTGACGGATAGCGAGGGTGAGATTGTTTGGCAGGCGACGTATCGCTCTTGGGGCACGATAGAAAc     <  2:1337021/150‑1 (MQ=255)
                                                                                                                                              aTTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGAGCTGACGGATAGCGAGGGTGAGATTGTTTGGCAGGCGACGTATCGCTCTTGGGGCACGATAGAAACttt  <  2:2144427/150‑1 (MQ=255)
                                                                                                                                              aTTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGAGCTGACGGATAGCGAGGGTGAGATTGTTTGGCAGGCGACGTATCGCTCTTGGGGCACGATAGAAACttt  <  2:1962954/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
TGGGGAGGCGAATAGCCAAGCAGGCTGAAATCAATGGCGAGGTCGAGCAGAAGCGCTTCCTCTGGCAAGGGCTAAGGATGCTGCGGGAGGAAACGCCAGGGCAGAGCATCTTGTACCTGTACGAGCCGGGTAGCTATGCGCCATTGGCGCGGGTTGATCAGGTTGAAGGAGAAGGGCAGAAGGTCTATTACTTCCATACCGATCAGATTGGGACGCCGCTTGAGCTGACGGATAGCGAGGGTGAGATTGTTTGGCAGGCGACGTATCGCTCTTGGGGCACGATAGAAACTTT  >  NC_002947/3514707‑3514998

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