Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,317,344 G→C 9.1% Q252H (CAG→CAC)  PP_1150 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,317,3440GC 9.1% 199.0 / 10.9 66Q252H (CAG→CAC) PP_1150membrane protein
Reads supporting (aligned to +/- strand):  ref base G (30/30);  new base C (3/3);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.46e-01

GAAAGCCTGGGGCGCATCGTTGCCATCGACAGCCAACGTGAACATGCCTGGTTCGAGGGCAATCGCGGGCGCCAACGTGCCCGTGCCATTCGTGGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGCCCAACCCGATCAGCCCAGCCTGCTGTTGTTGCGTCAGCGAATCTGGGATGCCGCCCACGATGAGCAAATCAGT  >  NC_002947/1317195‑1317488
                                                                                                                                                     |                                                                                                                                                
gAAAGCCTGGGGCGCATCGTTGCCATCGACAGCCAACGTGAACATGCCTGGTTCGAGGGCAATCGCGGGCGCCAACGTGCCCGTGCCATTCGTGGCCTTAGCCAGAAACTGATGGTGATGTTCAGTATTTCCCGTTCGGTGCACCGCCAg                                                                                                                                                  >  1:786020/1‑150 (MQ=255)
 aaaGCCTGGGGCGCATCGTTGCCATCGACAGCCAACGTGAACATGCCTGGTTCGAGGGCAATCGCGGGCGCCAACGTGCCCGTGCCATTCGTGGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGt                                                                                                                                                 >  1:166273/1‑150 (MQ=255)
      cTGGGGCGCATCGTTGCCATCGACAGCCAACGTGAACATGCCTGGTTCGAGGGCAATCGCGGGCGCCAACGTGCCCGTGCCATTCGTGGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTg                                                                                                                                                >  1:589764/1‑146 (MQ=255)
      cTGGGGCGCATCGTTGCCATCGACAGCCAACGTGAACATGCCTGGTTCGAGGGCAATCGCGGGCGCCAACGTGCCCGTGCCATTCGTGGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTg                                                                                                                                                <  2:589764/146‑1 (MQ=255)
      cTGGGGCGCATCGTTGCCATCGACAGCCAACGTGAACATGCCTGGTTCGAGGGCAATCGCGGGCGCCAACGTGCCCGTGCCATTCGTGGCCTTAGCCAGAAACTGATGGTCTTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTggcgg                                                                                                                                            >  1:1082063/1‑150 (MQ=255)
      cTGGGGCGCATCGTTGCCATCGACAGCCAACGTGAACATGCCTGGCTAGAGGGCAATCGCGGGCGCCAACGTGCCCGTGCCATTCGTGGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTggcgg                                                                                                                                            >  2:709849/1‑150 (MQ=255)
           gcgcATCGTTGCCATCGACAGCCAACGTGAACATGCCTGGTTCGAGGGCAATCGCGGGCGCCAACGTGCCCGTGCCATTCGTGGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCACTggcggcgcac                                                                                                                                       <  1:386744/150‑5 (MQ=255)
            cgcATCGTTGCCATCGACAGCCAACGTGAACATGCCTGGTTCGAGGGCAATCGCGGGCGCCAACGTGCCCGTGCCATTCGTGGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTg                                                                                                                                      >  2:427514/1‑150 (MQ=255)
                 cGTTGCCATCGACAGCCAACGTGAACATGCCTGGTTCGAGGGCAATCGCGGGCGCCAACGTGCCCGTGCCATTCGTGGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGa                                                                                                                                 <  2:754666/150‑1 (MQ=255)
                  gTTGCCATCGACAGCCAACGTGAACATGCCTGGTTCGAGGGCAATCGCGGGCGCCAACGTGCCCGTGCCATTCGTGGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAg                                                                                                                                >  1:1026684/1‑150 (MQ=255)
                   ttGCCATCGACAGCCAACGCGAACATGCCTGGTTAGAGGGCAATCGACCTCGTCGACGTTGCCGTGCCATTCGTGGCCTTAGCCAGAAACTCATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAgg                                                                                                                               <  2:578388/150‑1 (MQ=255)
                      ccATCGACAGCCAACGTGAACATGCCTGGTTCGAGGGCAATCGCGGGCGCCAACGTGCCCGTGCCATTCGTGGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTg                                                                                                                            >  2:139509/1‑150 (MQ=255)
                       cATCGACAGCCAACGTGAACATGCCTGGTTCGAGGGCAATCGCGGGCGCCAACGTGCCCGTGCCATTCGTGGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGa                                                                                                                           >  2:228873/1‑150 (MQ=255)
                        aTCGACAGCCAACGTGAACATGCCTGGTTCGAGGGCAATCGCGGGCGCCAACGTGCCCGTGCCATTCGTGGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAg                                                                                                                          >  1:768451/1‑150 (MQ=255)
                           gACAGCCAACGTGAACATGCCTGGTTCGAGGGCAATCGCGGGCGCCAACGTGCCCGTGCCATTCGTGGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTa                                                                                                                       <  1:4897/150‑1 (MQ=255)
                             cAGCCAACGTGAACATGCCTGGTTCGAGGGCAATCGCGGGCGCCAACGTGCCCGTGCCATTCGTGGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGACAGTGGCGGCAACTGGATGAGGGTGAGGTaga                                                                                                                     <  1:422669/150‑1 (MQ=255)
                                   aCGTGAACATGCCTGGTTCGAGGGCAATCGCGGGCGCCAACGTGCCCGTGCCATTCGTGGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCt                                                                                                               <  2:166273/150‑1 (MQ=255)
                                        aaCATGCCTGGTTCGAGGGCAATCGCGGGCGCCAACGTGCCCGTGCCATTCGTGGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGc                                                                                                          <  1:313568/150‑1 (MQ=255)
                                           aTGCCTGGTTAGAGGGACATCGAGGGAGCAAACATGCCCGTGCCATTCGTGATCTTAGACAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCt                                                                                                       <  2:57414/150‑1 (MQ=255)
                                            tGCCTGGTTCGAGGGCAATCGCGGGCGCCAACGTGCCCGTGCCATTCGTGGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGa                                                                                                             <  1:1127274/143‑1 (MQ=255)
                                            tGCCTGGTTCGAGGGCAATCGCGGGCGCCAACGTGCCCGTGCCATTCGTGGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGa                                                                                                             >  2:1127274/1‑143 (MQ=255)
                                              ccTGGTTCGAGGGCAATCGCGGGCGCCAACGTGCCCGTGCCATTCGTGGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGt                                                                                                    >  2:1072112/1‑150 (MQ=255)
                                                  gTTCGAGGGCAATCGCGGGCGCCAACGTGCCCGTGCCATTCGTGGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCa                                                                                                <  1:427514/150‑1 (MQ=255)
                                                  gTTCGAGGGCAATCGCGGGCGCCAACGTGCCCGTGCCATTCGTGGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCa                                                                                                <  1:923905/150‑1 (MQ=255)
                                                  gTTCGAGGGCAATCGCGGGCGCCAACGTGCCCGTGCCATTCGTGGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCa                                                                                                >  2:923905/1‑150 (MQ=255)
                                                  gTTCGAGGGCAATCGCGGGCGCCAACGTGCCCGTGCCATGCGGGGCGGTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCa                                                                                                >  2:770802/1‑150 (MQ=255)
                                                              tCGCGGGCGCCAACGTGCCCGTGCCATTCGTGGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTg                                                                                                                            <  2:650223/110‑1 (MQ=255)
                                                              tCGCGGGCGCCAACGTGCCCGTGCCATTCGTGGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTg                                                                                                                            >  1:650223/1‑110 (MQ=255)
                                                                    gcgcCAACGTGCCCGTGCCATTCGTCGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACgctgct                                                                              >  1:567866/1‑150 (MQ=255)
                                                                      gcCAACGTGCCCGTGCCATTCGTGGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTgg                                                                            <  1:1072112/150‑1 (MQ=255)
                                                                       ccAACGTGCCCGTGCCATTCGTGGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGc                                                                           <  1:353594/150‑1 (MQ=255)
                                                                       ccAACGTGCCCGTGCCATTCGTGGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGc                                                                           <  1:228873/150‑1 (MQ=255)
                                                                        cAACGTGCCCGTGCCATTCGTGGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGc                                                                                                 >  1:43105/1‑127 (MQ=255)
                                                                        cAACGTGCCCGTGCCATTCGTGGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGc                                                                                                 <  2:43105/127‑1 (MQ=255)
                                                                             tgcccgtgccATTCGTCGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGCCCAAcc                                                                     <  2:567866/150‑1 (MQ=255)
                                                                                   tgccATTCGTGGCCTTAGCCCGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGCCCAACCCGATca                                                               >  1:565114/1‑150 (MQ=255)
                                                                                     ccATTCGTTGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTCCTGCGACGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGCCCAACCCGATcagc                                                             <  1:68407/150‑1 (MQ=255)
                                                                                      cATTCGTGGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGCCCAACCCGATcaacc                                                            >  2:1100343/1‑150 (MQ=255)
                                                                                      cATTCGGGGCCTGAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGCCCAACCCGATcagcc                                                            >  2:856439/1‑150 (MQ=255)
                                                                                       aTTCGTGGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGCCCAACCCGATcagccc                                                           >  1:236223/1‑150 (MQ=255)
                                                                                                  tAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGCCCAACCCGATCAGCCCAGCCTGCtgtt                                                <  1:974045/150‑1 (MQ=255)
                                                                                                     ccAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGGTGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGCCCAACCCGATCAGCCCAGCCTGCTgttgtt                                             >  1:423963/1‑150 (MQ=255)
                                                                                                     ccAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGACGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGCCCAACCCGATCAGCCCAGCCTGCTgttgtt                                             <  2:768451/150‑1 (MQ=255)
                                                                                                             tGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGCCCAACCCGATCAGCCCAGCCTGCTGTTGTTGCGTCAGc                                     >  2:859095/1‑150 (MQ=255)
                                                                                                              gATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGCCCAACCCGATCAGCCCAGCCTGCTGTTGTTGCGTCAGCg                                    >  2:494837/1‑150 (MQ=255)
                                                                                                                    gCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGCCCAACCCGATCAGCCCAGCCTGCTGTTGTTGCGTCAGCGAATCTg                              <  1:709849/150‑1 (MQ=255)
                                                                                                                     cTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGCCCAACCCGATCAGCCCAGCCTGCTGTTGTTGCGTCAGCGAATCTgg                             <  2:1026684/150‑1 (MQ=255)
                                                                                                                     cTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGCCCAACCCGATCAGCCCAGCCTGCTGTTGTTGCGTCAGCGAATCTgg                             >  1:710820/1‑150 (MQ=255)
                                                                                                                     cTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGCCCAACCCGATCAGCCCAGCCTGCTGTTGTTGCGTCAGCGAATCTgg                             <  1:859095/150‑1 (MQ=255)
                                                                                                                     cTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGCCCAACCCGATCAACCCAGCCTGCTGTTGTTGCGTCAGCGAATCTgg                             <  1:1100343/150‑1 (MQ=255)
                                                                                                                     cTGTGGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGCCCAACCCGATCAGCCCAGCCTGCTGTTGTTGCGTCAGCGAATCTgg                             <  1:856439/150‑1 (MQ=255)
                                                                                                                      tGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGCCCAACCCGATCAGCCCAGCCTGCTGTTGTTGCGTCAGCGAATCTggg                            <  2:190635/150‑1 (MQ=255)
                                                                                                                      tGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGGGCGCACGCTGCTGGCCCAACCCGATCAGCCCAGCCTGCTGTTGTTGCGTCAGCGAATCTggg                            <  2:191219/150‑1 (MQ=255)
                                                                                                                          gCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGCCCAACCCGATCAGCCCAGCCTGCTGTTGTTGCGTCAGCGAATCTGGGATgc                        >  1:1109722/1‑150 (MQ=255)
                                                                                                                              aTTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGATACTGGCCCAACCCGATCAGCCCAGCCTGCTGTTGTTGCGTCAGCGAATCTGGGATgccgcc                    >  2:145853/1‑150 (MQ=255)
                                                                                                                               tttCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGCCCAACCCGATCAGCCCAGCCTGCTGTTGTTGCGTCAGCGAATCTGGGATGCCGccc                   <  1:494837/150‑1 (MQ=255)
                                                                                                                                  cccGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGCCCAACCCGATCAGCCCAGCCTGCTGTTGTTGCGTCAGCGAATCTGGGATGCCGCCCACg                >  1:396525/1‑150 (MQ=255)
                                                                                                                                      ttCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGCCCAACCCGATCAGCCCAGCCTGCTGTTGTTGCGTCAGCGAATCTGGGATGCCGCCCACTATGa            <  1:632878/150‑1 (MQ=255)
                                                                                                                                          gTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGCCCAACCCGATCAGCCCAGCCTGCTGTTGTTGCGTCAGCGAATCTGGGATGCCGCCCACGATGAGCaa        >  1:797973/1‑150 (MQ=255)
                                                                                                                                           gttgccgccACTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGCCCAACCCGATCAGCCCAGCCTGCTGTTGTTGCGTCAGCGAATCTGGGATGCCGCCCACGATGAGCaaa       >  2:615565/4‑150 (MQ=255)
                                                                                                                                            tttccgccACTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGCCCAACCCGATCAGCCCAGCCTGCtgttgt                                              <  1:378621/107‑1 (MQ=255)
                                                                                                                                            tttccgccACTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGCCCAACCCGATCAGCCCAGCCTGCtgttgt                                              >  2:378621/4‑110 (MQ=255)
                                                                                                                                              gccgccAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGCCCAACCCGATCAGCCCAGCCTGCTGTTGTTGCGTCAGCGAATCTGGGATGCCGCCCACGATGAGCAAATCa    <  2:236223/150‑1 (MQ=255)
                                                                                                                                               ccgccACTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCt                                                                                                               <  2:732685/42‑1 (MQ=255)
                                                                                                                                               ccgccACTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCt                                                                                                               >  1:732685/1‑42 (MQ=255)
                                                                                                                                                cgccAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGCCCAACCCGATCAGCCCAGCCTGCTGTTGTTGCGTCAGCGAATCTGGGATGCCGCCCACGATGCGCAAATCAGt  >  1:81660/1‑150 (MQ=255)
                                                                                                                                                cgccAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGCCCAACCCGATCAGCCCAGCCTGCTGTTGTTGCGTCAGCGAATCTGGGATGCCGCCCACGATGAGCAAATCAGt  >  2:798524/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GAAAGCCTGGGGCGCATCGTTGCCATCGACAGCCAACGTGAACATGCCTGGTTCGAGGGCAATCGCGGGCGCCAACGTGCCCGTGCCATTCGTGGCCTTAGCCAGAAACTGATGGTGCTGTTGCGTATTTCCCGTTCGGTGCGCCGCCAGTGGCGGCAACTGGATGAGGGTGAGGTAGAGCATCTGACGCCCTGGTTGCAGGAAGTGCGCACGCTGCTGGCCCAACCCGATCAGCCCAGCCTGCTGTTGTTGCGTCAGCGAATCTGGGATGCCGCCCACGATGAGCAAATCAGT  >  NC_002947/1317195‑1317488

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