Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,621,330 G→C 7.6% D49H (GAT→CAT) ‡ PP_2294 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,621,3300GC 7.6% 204.9 / 6.0 66D49H (GAT→CAT) ‡PP_2294hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/29);  new base C (2/3);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

TGTCGAGACGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCTAGACGGGAAAG  >  NC_002947/2621190‑2621477
                                                                                                                                            |                                                                                                                                                    
tgTCGAGACGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAg                                                                                                                                              >  4:209480/1‑149 (MQ=255)
 gTCGAGACGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAgg                                                                                                                                             >  5:56126/1‑149 (MQ=255)
  tCGAGACGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGa                                                                                                                                                     <  3:301972/140‑1 (MQ=255)
  tCGAGACGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGa                                                                                                                                                     >  4:301972/1‑140 (MQ=255)
    gagaCGGGCGATTGGGATGGTGTTGGCCCTTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCATCAAGCTGAgcata                                                                                                                                          <  4:14706/149‑5 (MQ=255)
     agaCGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTgg                                                                                                                                        >  8:249706/1‑149 (MQ=255)
     agaCGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCATCAAGCTGAgcatag                                                                                                                                         <  1:275703/149‑6 (MQ=255)
       aCGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCt                                                                                                                                      <  4:147156/149‑1 (MQ=255)
         gggCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTACCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTg                                                                                                                                     >  1:183094/1‑148 (MQ=255)
           gCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCaa                                                                                                                                  <  4:223327/149‑1 (MQ=255)
               ttGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGc                                                                                                                              <  3:209480/149‑1 (MQ=255)
                    aTCGTGTTGGCCATTGGAGTTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATa                                                                                                                         <  4:326888/149‑1 (MQ=255)
                    aTCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATa                                                                                                                         <  2:312276/149‑1 (MQ=255)
                          ttGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTg                                                                                                                    <  6:56126/148‑1 (MQ=255)
                          ttGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGa                                                                                                                   >  7:118912/1‑149 (MQ=255)
                               cATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGAt                                                                                                                  <  6:158120/145‑1 (MQ=255)
                               cATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGAt                                                                                                                  >  5:158120/1‑145 (MQ=255)
                                aTTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCagaa                                                                                                             <  8:118912/149‑1 (MQ=255)
                                          ctTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCAGAAAGCGCAGAGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGa                                                                                                    >  7:32992/1‑148 (MQ=255)
                                           ttATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATAAGAAGAAGGATGACa                                                                                                  >  5:390136/1‑149 (MQ=255)
                                              tgttgtCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACAttt                                                                                               >  6:88022/1‑149 (MQ=255)
                                                  gtCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAag                                                                                                                                 >  4:375739/1‑111 (MQ=255)
                                                  gtCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAag                                                                                                                                 <  3:375739/111‑1 (MQ=255)
                                                    ccTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGtt                                                                                                                                           <  7:255703/100‑1 (MQ=255)
                                                    ccTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGtt                                                                                                                                           >  8:255703/1‑100 (MQ=255)
                                                     cTCAAAGGAGTAGATAACGGGCCTAGCTCTTCCTCTTCAAGAGATCCAGAGGTTTATCCCGAAAGCGCAGTGACCAATGCTAAGCTTGAGGCTGAGGTTGGATGAAAGATAATATACCGTGATCAGAAAAAGGATCACATTTTGGATTc                                                                                        >  2:166813/1‑149 (MQ=255)
                                                      tCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTcc                                                                                       >  2:94911/1‑149 (MQ=255)
                                                      tCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTcc                                                                                       >  3:187635/1‑149 (MQ=255)
                                                        aaaGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATACCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCctc                                                                                     >  1:369404/1‑149 (MQ=255)
                                                          aGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGAt                                                                                                      <  2:271397/130‑1 (MQ=255)
                                                          aGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGAt                                                                                                      >  1:271397/1‑130 (MQ=255)
                                                                  cTAACGCGCCTAGCTCTTCCTCTTCAAATGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACaaaaa                                                                            >  1:208691/1‑148 (MQ=255)
                                                                    aaCGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATc                                                                          >  6:345234/1‑148 (MQ=255)
                                                                     aCGCGCCTAGCTCTTCCTCTTCAATTTCTCCAGAGGTTTATGCCGAAAGAGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCAc                                                                        <  6:329854/149‑1 (MQ=255)
                                                                     aCGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGACCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAAt                                                                           <  1:298532/147‑1 (MQ=255)
                                                                              gCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGAcc                                                               <  5:345234/149‑1 (MQ=255)
                                                                                tcttcctcttcAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTcc                                                                                       >  6:383608/1‑123 (MQ=255)
                                                                                tcttcctcttcAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTcc                                                                                       <  5:383608/123‑1 (MQ=255)
                                                                                  ttcctcttcAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCt                                                                                                                                      <  1:148626/74‑1 (MQ=255)
                                                                                  ttcctcttcAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCt                                                                                                                                      >  2:148626/1‑74 (MQ=255)
                                                                                       cttcAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCAt                                                             <  5:99747/142‑1 (MQ=255)
                                                                                       cttcAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCAt                                                             >  6:99747/1‑142 (MQ=255)
                                                                                         tcAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCt                                                     >  4:301230/1‑148 (MQ=255)
                                                                                               gCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTg                                               <  1:94911/148‑1 (MQ=255)
                                                                                                cTCCAGAGTTTTATCCCTAAAGCGCAGTGACCTATTCTCAGATTTATGATTAGGTTGGCTGCAATAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGAc                                             <  1:271842/149‑1 (MQ=255)
                                                                                                   cAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCtt                                          >  2:63515/1‑149 (MQ=255)
                                                                                                     gagGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTtgt                                        <  5:88022/149‑1 (MQ=255)
                                                                                                     gagGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTtgt                                        <  3:301230/149‑1 (MQ=255)
                                                                                                         tttATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGAt                                    <  4:176295/149‑1 (MQ=255)
                                                                                                           tATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGCAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGAtt                                   >  8:267369/1‑148 (MQ=255)
                                                                                                               ccGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCt                                                     >  5:358500/1‑126 (MQ=255)
                                                                                                               ccGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCt                                                     <  6:358500/126‑1 (MQ=255)
                                                                                                                cGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGAcc                                                               <  7:265334/115‑1 (MQ=255)
                                                                                                                cGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGAcc                                                               >  8:265334/1‑115 (MQ=255)
                                                                                                                 gAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGGAAGATCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACaa                            <  2:183094/149‑1 (MQ=255)
                                                                                                                       gcAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATc                      >  6:314470/1‑149 (MQ=255)
                                                                                                                                    cTCAGCATCAAGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGa                                     >  6:129064/1‑121 (MQ=255)
                                                                                                                                    cTCAGCATCAAGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGa                                     <  5:129064/121‑1 (MQ=255)
                                                                                                                                    cTCAGCATCAAGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCTAGAc         >  1:111458/1‑149 (MQ=255)
                                                                                                                                     tCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGATGAAGGATGACATTTTCGATTCCTCTTTCAAAAATCACTGGATGACCATGAAAGGCTAAATTGACCTTGTGGATTCATCCAACGTATAGCTCAACCTAGACg        >  7:110871/1‑149 (MQ=255)
                                                                                                                                     tCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTtgt                                        >  6:376119/1‑117 (MQ=255)
                                                                                                                                     tCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTtgt                                        <  5:376119/117‑1 (MQ=255)
                                                                                                                                       aGCTTGATGCTTAGGTTGGCTGCAAGAGCAAATACAGTGATCAGATTAATGATTACATTTTCGATTCCTCTTACAAAAATCAATGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCTAGACgg       <  7:158411/148‑1 (MQ=255)
                                                                                                                                        gCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTtgt                                        >  2:23845/1‑114 (MQ=255)
                                                                                                                                        gCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTtgt                                        <  1:23845/114‑1 (MQ=255)
                                                                                                                                         cTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCTAGACGGGaa    <  7:249706/149‑1 (MQ=255)
                                                                                                                                            gATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCTAGACGGGAAAg  <  2:369404/148‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
TGTCGAGACGGGCGATTGGGATCGTGTTGGCCATTGGAGCTGCTTATGTTGTCCTCAAAGGAGTAGCTAACGCGCCTAGCTCTTCCTCTTCAAGTGCTCCAGAGGTTTATCCCGAAAGCGCAGTGACCTATGCTCAGCTTGATGCTGAGGTTGGCTGCAAGAGCAAATACAGTGATCAGAAGAAGGATGACATTTTCGATTCCTCTTACAAAAATCACTGGATGACCATGAAAGGCTCAATTGACCTTGTGGATTCAGACAACGTATCGCTCAACCTAGACGGGAAAG  >  NC_002947/2621190‑2621477

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