Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,750,735 G→T 7.5% H89N (CAT→AAT)  PP_3315 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,750,7350GT 7.5% 196.6 / 4.2 62H89N (CAT→AAT) PP_3315membrane protein
Reads supporting (aligned to +/- strand):  ref base G (23/34);  new base T (3/2);  total (26/36)
Fisher's exact test for biased strand distribution p-value = 6.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

GCGCCAGCGAGCGCCGCGACTTTGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTACCAGCGACACCACGATAACCACCG  >  NC_002947/3750590‑3750881
                                                                                                                                                  |                                                                                                                                                  
gcgcCAGCGAGCGCCGCGACTTTGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGTATGCtt                                                                                                                                                 <  2:67081/149‑1 (MQ=255)
gcgcCAGCGAGCGCCGCGACTTTGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGcgt                                                                                                                                                 >  3:262177/1‑147 (MQ=255)
gcgcCAGCGAGCGCCGCGACTTTGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTAAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCtt                                                                                                                                                 <  1:231118/149‑1 (MQ=255)
    cAGCGAGCGCCGCGACTTTGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAAt                                                                                                                                             >  5:208080/1‑149 (MQ=255)
     agcgagcgCCGCGACTTTGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAAtt                                                                                                                                            <  4:283492/149‑1 (MQ=255)
          gcgcCGCGACTTTGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATTCAAGCATTCCATc                                                                                                                                       <  2:215751/149‑1 (MQ=255)
          gcgcCGCGACTTTGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATc                                                                                                                                       <  3:94175/149‑1 (MQ=255)
           cgccgcGACTTTGGCCGGTCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAAGACAACGCGCACGATGGAATGCTTGAATTCCATc                                                                                                                                       <  1:252596/149‑1 (MQ=255)
           cgccgcGACTTTGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATcg                                                                                                                                      <  6:208080/149‑1 (MQ=255)
            gccgcGACTTTGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATcgc                                                                                                                                     <  2:352055/149‑1 (MQ=255)
                cgACTTTGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATc                                                                                                                                 <  4:81173/149‑1 (MQ=255)
                    tttGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATcag                                                                                                                               >  2:58903/1‑147 (MQ=255)
                         ccGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTc                                                                                                                        <  5:47283/149‑1 (MQ=255)
                         ccGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTc                                                                                                                        <  5:91387/149‑1 (MQ=255)
                              cTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGAc                                                                                                                   <  8:210304/149‑1 (MQ=255)
                              cTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGAc                                                                                                                   <  8:204753/149‑1 (MQ=255)
                                 aCATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTcc                                                                                                                                          >  6:126593/1‑123 (MQ=255)
                                   aTCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATAc                                                                                                              <  7:291015/149‑1 (MQ=255)
                                   aTCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATAc                                                                                                              >  2:275198/1‑149 (MQ=255)
                                   aTCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATAc                                                                                                              >  6:372520/1‑149 (MQ=255)
                                     cGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATAcc                                                                                                             >  8:85717/1‑148 (MQ=255)
                                              ggATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTc                                                                                                   <  7:85717/149‑1 (MQ=255)
                                              ggATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTc                                                                                                   <  1:4777/149‑1 (MQ=255)
                                                     gATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTgg                                                                                                >  3:71463/1‑145 (MQ=255)
                                                           gcgcGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATc                                                                                          <  8:93351/145‑1 (MQ=255)
                                                           gcgcGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCg                                                                                        <  6:312185/147‑1 (MQ=255)
                                                                   gcgcGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGc                                                                              <  5:372520/149‑1 (MQ=255)
                                                                     gcgATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCCTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCg                                                                             >  6:123085/1‑148 (MQ=255)
                                                                        aTCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGCTGGACCACCGGAtc                                                                         <  5:63514/149‑2 (MQ=255)
                                                                           aaCCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGaa                                                                      >  3:315008/1‑149 (MQ=255)
                                                                           aaCCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGaa                                                                      <  2:81208/149‑1 (MQ=255)
                                                                           aaCCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGCTGGACCACCGGAtcttc                                                                      <  1:58903/149‑5 (MQ=255)
                                                                             cccGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATc                                                                                                                                       <  4:99268/82‑1 (MQ=255)
                                                                             cccGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATc                                                                                                                                       >  3:99268/1‑82 (MQ=255)
                                                                             cccGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAAc                                                                     >  3:103403/1‑148 (MQ=255)
                                                                             cccGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACg                                                                    >  8:84586/1‑149 (MQ=255)
                                                                             cccGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGCTGGACCACCGGAtcttcca                                                                    <  6:110462/149‑7 (MQ=255)
                                                                               cGAGCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCa                                                                  <  4:44132/149‑1 (MQ=255)
                                                                                   cAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAAc                                                                     <  5:160264/142‑1 (MQ=255)
                                                                                   cAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAAc                                                                     >  6:160264/1‑142 (MQ=255)
                                                                                                     cTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGa                                            >  8:338478/1‑149 (MQ=255)
                                                                                                     cTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGa                                            >  7:191298/1‑149 (MQ=255)
                                                                                                               ctgcCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGAt                                  <  7:84586/149‑1 (MQ=255)
                                                                                                                  ccATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAg                               <  5:380843/149‑1 (MQ=255)
                                                                                                                  ccATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAg                               <  1:11460/149‑1 (MQ=255)
                                                                                                                  ccATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAg                               >  2:11460/1‑149 (MQ=255)
                                                                                                                  ccATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGAAGGCTGAAGACTATCTGGATCAg                               >  6:380843/1‑149 (MQ=255)
                                                                                                                   cATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGc                              <  1:275198/149‑1 (MQ=255)
                                                                                                                       gcgcAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGt                           <  2:77495/148‑1 (MQ=255)
                                                                                                                        cgcAAGAGAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATAGCGAAGACATTCTGGAAGATGCGGTGGTCCAGCGGATTGATCGCATCGATGATCAGCTGGCTGAAGACTATCTGGATCAGCTGTAc                         <  5:123085/149‑1 (MQ=255)
                                                                                                                                 aaCGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTACCAGCGaca                 >  5:32044/1‑149 (MQ=255)
                                                                                                                                 aaCGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTACCAGCGaca                 >  1:303398/1‑149 (MQ=255)
                                                                                                                                   cgcgCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAgg                                                   >  5:77151/1‑113 (MQ=255)
                                                                                                                                   cgcgCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAgg                                                   <  6:77151/113‑1 (MQ=255)
                                                                                                                                       cACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTACCAGCGACACCACGa           <  4:71463/149‑1 (MQ=255)
                                                                                                                                         cGATGGAATTCAAGCATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTACCAGCGACACCACGAt          >  1:60897/1‑148 (MQ=255)
                                                                                                                                         cGATGGAATTCAAGCATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTACCAGCGACACCACGAt          >  1:150491/1‑148 (MQ=255)
                                                                                                                                         cGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTACCAGCGACACCACGAt          >  2:310310/1‑148 (MQ=255)
                                                                                                                                          gATGGAATTCAAGCATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTACCAGCGACACCACGATaa        >  6:109899/1‑149 (MQ=255)
                                                                                                                                               aaTTCTTGAAATCCATCGCGATCAGCAGCGTCATGACCAAACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTACCAGCGACACCACGATAaccacc   <  6:32044/149‑1 (MQ=255)
                                                                                                                                               aaTGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTACCAGCGACACCACGATAaccacc   <  2:303398/149‑1 (MQ=255)
                                                                                                                                                aTGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTACCAGCGACACCACGATAACCACcg  >  3:119515/1‑149 (MQ=255)
                                                                                                                                                  gCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTACCAGCGACACCACGATAac       <  1:37935/142‑1 (MQ=255)
                                                                                                                                                  gCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTACCAGCGACACCACGATAac       >  2:37935/1‑142 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCGCCAGCGAGCGCCGCGACTTTGGCCGGGCTGACATCGGGATCAAGGATCACGATCTTGCGCGACAGCGCGATCAACCCGATCAAGAGGACCGTCTTGACCTGAATGATGCTGCCATGGCGCAAGACAACGCGCACGATGGAATGCTTGAATTCCATCGCGATCAGCAGCGTCATGACCATACCGAACACATTCTGGAAGATCCGGTGGTCCAGCGGATTGAACGCATCGATGATCAGCAGGCTGAAGACTATCTGGATCAGCTGTACCAGCGACACCACGATAACCACCG  >  NC_002947/3750590‑3750881

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: