Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,310,130 A→T 5.5% L470Q (CTG→CAG)  PP_1144 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,310,1300AT 5.5% 283.8 / 5.1 90L470Q (CTG→CAG) PP_1144membrane protein
Reads supporting (aligned to +/- strand):  ref base A (46/39);  new base T (3/2);  total (49/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

GCAATGCCGACAAAACGGTCGGGCCCGAGCAGGCCGTGTTGAGGATGCTGCCAGCGTACCAGCGCTTCCAGTTTGGCCACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGATTGGTAGAGGGCCAGGTCGGCATGCTTGAGCAGCTCTTCGGCTG  >  NC_002947/1309983‑1310268
                                                                                                                                                   |                                                                                                                                          
gCAATGCCGACAAAACGGTCGGGCCCGAGCAGGCCGTGTTGAGGATGCTGCCAGCGTACCAGCGCTTCCCGTTTGGCCACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCCCGCTCTTCGCGCAACGCCAgg                                                                                                                                          >  1:447442/1‑150 (MQ=255)
         aCAAAACGGTCGGGCCCGAGCAGGCCGTGTTGAGGATGCTGCCAGCGTACCAGCGCTTCCAGTTTGGCCACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGctct                                                                                                                                 >  1:711330/1‑150 (MQ=255)
         aCAAAACGGTCGGGCCCGAGCAGGCCGTGTTGAGGATGCTGCCAGCGTACCAGCGCTTCCAGTTTGGCCACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTCCACGTACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGctct                                                                                                                                 >  1:1233125/1‑150 (MQ=255)
         aCAAAACGGTCGGGCCCGAGCAGGCCGTGTTGAGGATGCTGCCAGCGTACCAGCGCTTCCAGTTTGGCCACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTACTGCCCACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGctct                                                                                                                                 >  2:438407/1‑150 (MQ=255)
          cAAAACGGTCGGGCCCGAGCAGGCCGTGTTGAGGATGCTGCCAGCGTACCAGCGCTTCCAGTTTGGCCACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCtctt                                                                                                                                <  1:1723563/150‑1 (MQ=255)
           aaaaCGGTCGGGCCCGAGCAGGCCGTGTTGAGGATGCTGCCAGCGTACCAGCGCTTCCAGTTTGGCCACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCTGCGCCTGGCGCTGctcgaa                                                                                                                               <  1:1475395/150‑4 (MQ=255)
           aaaaCGGTCGGGCCCGAGCAGGCCGTGTTGAGGATGCTGCCAGCGTACCAGCGCTTCCAGTTTGGCCACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCtcttc                                                                                                                               >  1:1445672/1‑150 (MQ=255)
                 gTCGGGCCCGAGCAGGCCGTGTTGAGGATGCTGCCAGCGTACCAGCGCTTCCAGTTTGGCCACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAg                                                                                                                         <  2:1445672/150‑1 (MQ=255)
                      gCCCGAGCAGGCCGTGTTGAGGATGCTGCCAGCGTACCAGCGCTTCCAGTTTGGCCACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCa                                                                                                                    >  1:1414598/1‑150 (MQ=255)
                                 ccGTGTTGAGGATGCTGCCAGCGTACCAGCGCTTCCAGTTTGGCCACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGAc                                                                                                         <  2:1924196/150‑1 (MQ=255)
                                     gtTGAGGATGCTGCCAGCGTACCAGCGCTTCCAGTTTGGCCACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGAcgcgc                                                                                                     <  1:1059913/150‑1 (MQ=255)
                                        gAGGATGCTGCCAGCGTACCAGCGCTTCCAGTTTGGCCACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAg                                                                                                  >  1:1657130/1‑150 (MQ=255)
                                              gctgcCAGCGTACCAGCGCTTCCAGTTTGGCCACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCa                                                                                            <  1:730832/150‑1 (MQ=255)
                                              gctgcCAGCGTACCAGCGCTTCCAGTTTGGCCACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCa                                                                                            <  1:186519/150‑1 (MQ=255)
                                               ctgcCAGCGTACCAGCGCTTCCAGTTTGGCCACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAg                                                                                           <  1:399000/150‑1 (MQ=255)
                                                ttccAGCGTACCAGCGCTTCCAGTTTGGCCACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAg                                                                                           >  1:1290428/3‑149 (MQ=255)
                                                tgcCAGCGTACCAGCGCTTCCAGTTTGGCCACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAg                                                                                           <  2:1290428/149‑1 (MQ=255)
                                                  ccAGCGTACCAGCGCTTCCAGTTTGGCCACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCt                                                                                        <  2:1414598/150‑1 (MQ=255)
                                                    agggTACCAGCGCTTCCAGTTTGGCCACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGt                                                                                      >  1:1908323/4‑150 (MQ=255)
                                                                 ttCCAGTTTGGCCACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTTCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTATAGTTTCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAAtt                                                                         >  1:2330692/1‑150 (MQ=255)
                                                                   ccAGTTTGGCCACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCa                                                                       <  1:438407/150‑1 (MQ=255)
                                                                    cagggtGGCCACTTGCCCGCTGCGAAGGTCGAAAATCGAATGGGTATGCAGACACGGCCCGCGCTCTGCGAGCAGCGGCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAg                                                                      <  2:1233125/145‑1 (MQ=255)
                                                                    cAGTTTGGCCACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAg                                                                      <  1:668181/150‑1 (MQ=255)
                                                                     aGTTTGGCCACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAgg                                                                     >  2:38261/1‑150 (MQ=255)
                                                                     aGTTTGGCCACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGAAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAgg                                                                     >  1:164014/1‑150 (MQ=255)
                                                                      gTTTGGCCACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGc                                                                    >  1:455766/1‑150 (MQ=255)
                                                                      gTTTGGCCACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGCACCATCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGTCTTGAGGCGTTCATTGAGGAAATTCAGGc                                                                    >  1:43335/1‑150 (MQ=255)
                                                                          ggCCACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGtt                                                                <  1:1420931/150‑1 (MQ=255)
                                                                            ccACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCACGCCCTTCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAAATCAGGCTGTTgc                                                              >  1:1728267/1‑150 (MQ=255)
                                                                            ccACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTc                                                                                             <  1:1250316/119‑1 (MQ=255)
                                                                            ccACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTc                                                                                             >  2:1250316/1‑119 (MQ=255)
                                                                            ccACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCa                                                                       <  1:1865280/141‑1 (MQ=255)
                                                                            ccACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCa                                                                       >  2:1865280/1‑141 (MQ=255)
                                                                            ccACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTgc                                                              >  1:1860615/1‑150 (MQ=255)
                                                                             cACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGAcgcgc                                                                                                     >  2:2319685/1‑110 (MQ=255)
                                                                             cACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGAcgcgc                                                                                                     <  1:2319685/110‑1 (MQ=255)
                                                                                ttGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGccc                                                          >  2:203699/1‑150 (MQ=255)
                                                                                ttGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGccc                                                          >  2:182573/1‑150 (MQ=255)
                                                                                     cGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGAcgc                                                                                                       >  2:128446/1‑100 (MQ=255)
                                                                                     cGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGAcgc                                                                                                       <  1:128446/100‑1 (MQ=255)
                                                                                         gcgcAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCTCCCGTTGTCCtt                                                 <  2:1552423/150‑1 (MQ=255)
                                                                                         gcgcAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCtt                                                 <  2:1519471/150‑1 (MQ=255)
                                                                                           gcAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTgg                                               >  1:622446/1‑150 (MQ=255)
                                                                                            cAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGAc                                                                                                         <  2:1292507/91‑1 (MQ=255)
                                                                                            cAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGAc                                                                                                         >  1:1292507/1‑91 (MQ=255)
                                                                                               gTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGATTg                                           >  2:1112712/1‑150 (MQ=255)
                                                                                                 cGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCACGAGGCTTATATCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGATTGGt                                         >  1:1143907/1‑150 (MQ=255)
                                                                                                      aTCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGATTGGTAGAgg                                    <  2:1860615/150‑1 (MQ=255)
                                                                                                       tCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGATTGGTagag                                     <  2:1389066/148‑1 (MQ=255)
                                                                                                       tCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGATTGGTagag                                     >  1:1389066/1‑148 (MQ=255)
                                                                                                          gCTGGTAGTGCACAAACAGCCCGCGCTCTTAGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTTTTGCGCCCGTTGTCCTTGGATTGGTAGAGGGCCa                                >  2:762578/1‑150 (MQ=255)
                                                                                                                aGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGATTGGTAGAGGGCCAGGTCgg                          <  2:1425669/150‑1 (MQ=255)
                                                                                                                  tGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGATTGGTAGAGGGCCAGGTCGGCa                        >  1:2343385/1‑150 (MQ=255)
                                                                                                                   gaaaacacaGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGATTGGTAGAGGGCCAGGTCGGCAt                       <  2:447442/146‑1 (MQ=255)
                                                                                                                    cacacacaGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTcc                                                   <  1:92710/121‑1 (MQ=255)
                                                                                                                    cacacacaGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTcc                                                   >  2:92710/1‑121 (MQ=255)
                                                                                                                    cacacaaaGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGAc                                                                                                         <  1:231455/67‑1 (MQ=255)
                                                                                                                    cacacaaaGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGAc                                                                                                         >  2:231455/1‑67 (MQ=255)
                                                                                                                     acacacaGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAg                                                                                           >  2:1451660/1‑80 (MQ=255)
                                                                                                                     acacacaGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAg                                                                                           <  1:1451660/80‑1 (MQ=255)
                                                                                                                     acacacaGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGATTGGTAGAGGGCCAGGTCGGCATGc                     >  2:465286/1‑150 (MQ=255)
                                                                                                                       acacaGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCtt                                                                                              <  1:113295/75‑1 (MQ=255)
                                                                                                                       acacaGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCtt                                                                                              >  2:113295/1‑75 (MQ=255)
                                                                                                                       acacaGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGATTGGTAGAGGGCCAGGTCGGCATGCtt                   <  2:1908323/150‑1 (MQ=255)
                                                                                                                        cacaGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAg                                                                                           >  1:1099223/1‑77 (MQ=255)
                                                                                                                        cacaGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAg                                                                                           <  2:1099223/77‑1 (MQ=255)
                                                                                                                        cacaGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGATTg                                           >  1:349967/1‑125 (MQ=255)
                                                                                                                        cacaGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGATTg                                           <  2:349967/125‑1 (MQ=255)
                                                                                                                         acaGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGATTGGTAGAGGGCCAGGTCGGCATGCTTGa                 >  1:1826332/1‑150 (MQ=255)
                                                                                                                         acaGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGATTGGTAGAGGGCCAGGTCGGCATGCTTGa                 >  1:2379694/1‑150 (MQ=255)
                                                                                                                          caGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCtcttct                                                                                                                              <  2:2198088/40‑1 (MQ=255)
                                                                                                                          caGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCtcttct                                                                                                                              >  1:2198088/1‑40 (MQ=255)
                                                                                                                            gCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGATTGGTAGAGGGCCAGGTCGGCATGCTTGagca              <  1:38261/150‑1 (MQ=255)
                                                                                                                                 cgcTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGAtt                                            >  2:463384/1‑115 (MQ=255)
                                                                                                                                 cgcTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGAtt                                            <  1:463384/115‑1 (MQ=255)
                                                                                                                                  gcTCTTCGAGCAGCGCCAGGTGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGATTGGTa                                        >  1:2241459/1‑118 (MQ=255)
                                                                                                                                  gcTCTTCGAGCAGCGCCAGGTGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGATTGGTa                                        <  2:2241459/118‑1 (MQ=255)
                                                                                                                                  gcTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGATTGGTAGAGGGCCAGGTCGGCATGCTTGAGCAGCTCtt        >  1:2221615/1‑150 (MQ=255)
                                                                                                                                   ctctTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGATTGGTAGAGGGCCAGGTCGGCATGCTTGAGCAGCTCTTc       <  1:203699/150‑1 (MQ=255)
                                                                                                                                    tctTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGATTGGTa                                        >  2:1309699/1‑116 (MQ=255)
                                                                                                                                    tctTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGATTGGTa                                        <  1:1309699/116‑1 (MQ=255)
                                                                                                                                    tctTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGATTGGTAGAGGGCCAGGTCGGCATGCTTGAGCAGCTCTTCg      >  2:2150517/1‑150 (MQ=255)
                                                                                                                                     agaagAGCTGCGCCTGGCGCTGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGATTGGTAGAGGGCCAGGTCGGCATGCTTGAGCAGCTCTTCgg     >  2:920695/5‑150 (MQ=255)
                                                                                                                                      ttCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGATTgg                                          <  2:1725884/112‑1 (MQ=255)
                                                                                                                                      ttCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGATTgg                                          >  1:1725884/1‑112 (MQ=255)
                                                                                                                                       tCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGATTGGTAGAGGGCCAGGTCGGCATGCTTGAGCAGCTCTTCGGCt   <  1:465286/150‑1 (MQ=255)
                                                                                                                                        cGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGATTGGTAGAGGGCCAGGTCGGCATGCTTGAGCAGCTCTTCGGCtg  >  2:2259085/1‑150 (MQ=255)
                                                                                                                                        agAGCTGCGCCTGGCGCTGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGATTGGTAGAGGGCCAGGTCGGCATGCTTGAGCAGCTCTTCGGCtg  >  1:1818879/2‑150 (MQ=255)
                                                                                                                                                gcCTGGCGCTGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGATTGGTAGAGGGCCAGGTCGGCATGc                     <  1:1133952/123‑1 (MQ=255)
                                                                                                                                                gcCTGGCGCTGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGATTGGTAGAGGGCCAGGTCGGCATGc                     >  2:1133952/1‑123 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
GCAATGCCGACAAAACGGTCGGGCCCGAGCAGGCCGTGTTGAGGATGCTGCCAGCGTACCAGCGCTTCCAGTTTGGCCACTTGCCCGCTGCGCAGGTCGAAGATCGGCTGGTAGTGCACACACAGCCCGCGCTCTTCGAGCAGCGCCAGGCGCAGCTCTTCTTCCAGTTGCAAGGCGATCGACGCGCGAGCCTTCAGGCTGTCATTGAAGAAATTCAGGCTGTTGCGCCCGTTGTCCTTGGATTGGTAGAGGGCCAGGTCGGCATGCTTGAGCAGCTCTTCGGCTG  >  NC_002947/1309983‑1310268

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: