Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,507,789 T→G 11.6% F32C (TTC→TGC) ‡ PP_3105 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,507,7890TG11.6% 168.7 / 17.4 69F32C (TTC→TGC) ‡PP_3105membrane protein
Reads supporting (aligned to +/- strand):  ref base T (27/34);  new base G (4/4);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.92e-01

ATTTTTCGGGGCTGGGCACTACCCTGGGTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTCGGTTGCACCAATGTTCGCATCG  >  NC_002947/3507642‑3507928
                                                                                                                                                    |                                                                                                                                           
aTTTTTCGGGGCTGGGCACTACCCTGGGTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCa                                                                                                                                           <  2:547708/150‑1 (MQ=255)
  ttttCGGGGCTGGGCACTACCCTGGGTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATg                                                                                                                                         <  2:560021/150‑1 (MQ=255)
  ttttCGGGCCTGGGACTTACCTTGGTAATTTTTTGGAGGGGAAAGAATAAAAGGTTTCCACAAAACCGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATg                                                                                                                                         <  2:1660639/150‑1 (MQ=255)
   tttCGGGGCTGGGCACTACCCTGGGTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATgg                                                                                                                                        <  2:557782/150‑1 (MQ=255)
   ttcaggggCTGGAGACTAATCTGGGTAATTTATTGACGCTGAAGGATTAAATAATTACCAAAGAACGCTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATgg                                                                                                                                        <  2:558310/146‑1 (MQ=255)
       ggggCTGGGCACTACCCTGGGTAATTTCTTGACGGGTAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATGAATGGCggc                                                                                                                                    <  1:197383/150‑2 (MQ=255)
       ggggCTGGGCACTACCCTGGGTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGt                                                                                                                                    >  2:385241/1‑150 (MQ=255)
       ggggCTGGGCACTACCCTGGGTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGt                                                                                                                                    >  2:1742051/1‑150 (MQ=255)
           cTGGGCACTACCCTGGGTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCg                                                                                                                                <  2:1560378/150‑1 (MQ=255)
               gCACTACCCTGGGTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCAt                                                                                                                                          <  2:685224/136‑1 (MQ=255)
               gCACTACCCTGGGTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCAt                                                                                                                                          >  1:685224/1‑136 (MQ=255)
                     cccTGGGTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGGGCCGCCATTCCTGGCGGTGCCGGCGATCCTGa                                                                                                                      >  1:889687/1‑150 (MQ=255)
                       cTGGGTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATc                                                                                                                    <  2:1847082/150‑1 (MQ=255)
                         gggTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCAt                                                                                                                  >  1:39009/1‑150 (MQ=255)
                          ggTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCAtt                                                                                                                 >  1:353506/1‑150 (MQ=255)
                                ttCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGAGGTGCCGGCGATCCTGATCATTATCGtt                                                                                                           <  2:475540/150‑1 (MQ=255)
                                 tCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTa                                                                                                          <  2:353506/150‑1 (MQ=255)
                                 tCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTa                                                                                                          <  1:1012766/150‑1 (MQ=255)
                                  cTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAACTCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCTGGCGATCCTGATCATTATCGTTAt                                                                                                         <  1:1786099/150‑1 (MQ=255)
                                  cTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTAt                                                                                                         <  1:15583/150‑1 (MQ=255)
                                      aCGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACg                                                                                                     >  1:76123/1‑150 (MQ=255)
                                                       gTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTAt                                                                                                               <  1:1599375/123‑1 (MQ=255)
                                                       gTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTAt                                                                                                               >  2:1599375/1‑123 (MQ=255)
                                                           aCCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATtgtg                                                                                  >  2:1511964/1‑148 (MQ=255)
                                                           aCCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATtgtg                                                                                  <  1:1511964/148‑1 (MQ=255)
                                                              acaGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCa                                                                                           <  1:1363072/136‑1 (MQ=255)
                                                              acaGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCa                                                                                           >  2:1363072/1‑136 (MQ=255)
                                                                aGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCg                                                                           <  1:1742051/150‑1 (MQ=255)
                                                                     cGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTGGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGt                                                                       <  1:1476969/150‑1 (MQ=255)
                                                                      gTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAACTGTGCTTACCGCAGTg                                                                      >  1:1931881/1‑150 (MQ=255)
                                                                            ttAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACt                                                                <  2:76123/150‑1 (MQ=255)
                                                                              aaaaaaaTCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATggag                                                                                                                                      <  1:559684/78‑3 (MQ=255)
                                                                              aaaaaaaTCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCg                                                                                                                                      >  2:559684/1‑78 (MQ=255)
                                                                                     tCTGGTGCGGGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGACCACTTTGTAAGtt                                                       <  2:889687/150‑1 (MQ=255)
                                                                                     tCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGtt                                                       >  2:615280/1‑150 (MQ=255)
                                                                                        ggTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAAtt                                                                                     <  2:1047845/117‑1 (MQ=255)
                                                                                        ggTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAAtt                                                                                     >  1:1047845/1‑117 (MQ=255)
                                                                                             gtgtGGCTGCTGGTCGTGATATTAGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGaga                                               >  1:1868083/1‑148 (MQ=255)
                                                                                               gtgGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAAttt                                             <  2:1909958/150‑1 (MQ=255)
                                                                                                 ggCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGt                                                            <  2:331111/133‑1 (MQ=255)
                                                                                                 ggCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGt                                                            >  1:331111/1‑133 (MQ=255)
                                                                                                  gctgctGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCACCGCCATGAATGGCGGCGcca                                                                                                                                <  2:421458/64‑2 (MQ=21)
                                                                                                  gctgctGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCACCGCCATGAATGGCGGCGcca                                                                                                                                >  1:421458/1‑63 (MQ=25)
                                                                                                    tgctgGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCa                                        >  1:1612525/1‑150 (MQ=255)
                                                                                                           cGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGttctt                                 <  1:615280/150‑1 (MQ=255)
                                                                                                               atatTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAat                             <  1:385241/150‑1 (MQ=255)
                                                                                                                 atTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAatat                           >  2:986573/1‑150 (MQ=255)
                                                                                                                 atTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAACTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAatat                           <  2:1931881/150‑1 (MQ=255)
                                                                                                                    ggCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGAGCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTc                        >  2:1061937/1‑150 (MQ=255)
                                                                                                                      cccTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTAt                                                                                                         >  1:1044722/1‑67 (MQ=255)
                                                                                                                      cccTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTAt                                                                                                         <  2:1044722/67‑1 (MQ=255)
                                                                                                                         tgatgatgaCCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTCGGTTg                   >  1:1028903/1‑150 (MQ=255)
                                                                                                                         tgatgatgaCCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTCGGTTg                   >  2:752418/1‑150 (MQ=255)
                                                                                                                          gatgatgaCCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTCGGTTGc                  <  2:39009/150‑1 (MQ=255)
                                                                                                                           atgatgaCCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATc                                                                                                              >  1:322687/1‑57 (MQ=255)
                                                                                                                           atgatgaCCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATc                                                                                                              <  2:322687/57‑1 (MQ=255)
                                                                                                                                   cGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGt                         <  1:920664/134‑1 (MQ=255)
                                                                                                                                   cGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGt                         >  2:920664/1‑134 (MQ=255)
                                                                                                                                       gTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAg                                                                        <  1:727375/83‑1 (MQ=255)
                                                                                                                                       gTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAg                                                                        >  2:727375/1‑83 (MQ=255)
                                                                                                                                        cggcaccgccATGAATGGCGGCGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTCGGTTGCACCAATGtt         <  1:636075/140‑1 (MQ=255)
                                                                                                                                        cggcaccgccATGAATGGCGGCGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTCGGTTGCACCAATGtt         >  2:636075/6‑145 (MQ=255)
                                                                                                                                          gcaccgccATGAATGGCGGCGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTaa                                                          >  1:1435198/4‑94 (MQ=255)
                                                                                                                                          gcaccgccATGAATGGCGGCGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTaa                                                          <  2:1435198/91‑1 (MQ=255)
                                                                                                                                          gcaccgccATGAATGGCGGCGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTCGGTTGCACCAATGTTCGCATCg  >  1:1354713/4‑150 (MQ=255)
                                                                                                                                            gccgccATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGt                                                            >  2:1059142/1‑90 (MQ=255)
                                                                                                                                            gccgccATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGt                                                            <  1:1059142/90‑1 (MQ=255)
                                                                                                                                               gccATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTCgg                      <  1:911363/125‑1 (MQ=255)
                                                                                                                                               gccATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTCgg                      >  2:911363/1‑125 (MQ=255)
                                                                                                                                                    |                                                                                                                                           
ATTTTTCGGGGCTGGGCACTACCCTGGGTAATTTCTTGACGGGGAAGGATTAAATGTTTACCACAGAACGTTTTTTTAAAAAAATCTGGTCGGTGTGGCTGCTGGTCGTGATATTGGCCCTGATGATGACCGGAGTGGCGCCGCCATTCATGGCGGTGCCGGCGATCCTGATCATTATCGTTATGACGCTGTGGTGCATCAATTGTGCTTACCGCAGTGAGCACTTTGTAAGTTTCGCGAATTTGCGCATGTTCTTCAATATGTCGGTTGCACCAATGTTCGCATCG  >  NC_002947/3507642‑3507928

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