Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 872,089 C→A 5.9% W359L (TGG→TTG)  PP_0756 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947872,0890CA 5.9% 217.6 / 2.4 67W359L (TGG→TTG) PP_0756hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (29/34);  new base A (2/2);  total (31/36)
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.11e-01

GCATACTGCAACACCAGGTACGGCCGCTGCTCTTGCGGGCAGGCATCCAGGCGCCGTTGCCATTCTTCCTTGGCCTTGGCCAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATCTCTGCCGCCAGCTGTTTTGTGTCTTCGTGGTACTCGGTTACCGGCGGCGCGAAGTTGATGTGGACTCGGCCCTTGTAGCCGG  >  NC_002947/871945‑872236
                                                                                                                                                |                                                                                                                                                   
gCATACTGCAACACCAGGTACGGCCGCTGCTCTTGCGGGCAGGCATCCAGGCGCCGTTGCCATTCTTCCTTGGCCTTGGCCAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATc                                                                                                                                                >  2:116601/1‑150 (MQ=255)
  aTACTGCAACACCAGGTACGGCCGCTGCTCTTGCGGGCAGGCATCCAGGCGCCGTTGCCATTCTTCCTTGGCCTTGGCCAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGc                                                                                                                                              >  1:108843/1‑150 (MQ=255)
    aCTGCAACACCAGGTACGGCCGCTGCTCTTGCGGGCAGGCATCCAGGCGCCGTTGCCATTCTTCCTTGGCCTTGGCCAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCat                                                                                                                                            <  2:283811/150‑1 (MQ=255)
      tGCAACACCAGGTACGGCCGCTGCTCTTGCGGGCAGGCATCCAGGCGCCGTTGCCATTCTTCCTTGGCCTTGGCCAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCatat                                                                                                                                          <  1:99346/150‑1 (MQ=255)
       gCAACACCAGGTACGGCCGCTGCTCTTGCGGGCAGGCAGCCAGGCGCCGTTGCCATTCTTCCTTGGCCTTGGCCAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATg                                                                                                                                         <  1:1374257/150‑1 (MQ=255)
            acCAGGTACGGCCGCTGCTCTTGCGGGCAGGCATCCAGGCGCCGTTGCCATTCTTCCTTGGCCTTGGCCAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAgg                                                                                                                                    >  1:254583/1‑150 (MQ=255)
            acCAGGTACGGCCGCGGCTCTTGCGGCCAGGCATCCAGGCGCCGTTGCCATTCTTCCTTGGCCTTGGCCAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAgg                                                                                                                                    >  2:668606/1‑150 (MQ=255)
             ccAGGTACGGCCGCTGCTCTTGCGGGCAGGCATCCAGGCGCCGTTGCCATTCTTCCTTGGCCTTGGCCAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGt                                                                                                                                   >  2:572637/1‑150 (MQ=255)
             ccAGGTACGGCCGCTGCTCTTGCGGGCAGGCATCCAGGCGCCGTTGCCATTCTTCCTTGGCCTTGGCCAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGt                                                                                                                                   >  1:243464/1‑150 (MQ=255)
               aGGTACGGCCGCTGCTCTTGCGGGCAGGCATCCAGGCGCCGTTGCCATTCTTCCTTGGCCTTGACCAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTaa                                                                                                                                 <  2:182529/150‑1 (MQ=255)
                ggTACGGCCGCTGCTCTTGCGGGCAGGCATCCAGGCGCCGTTGCCATTCTTCCTTGGCCTTGGCCAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAAt                                                                                                                                <  2:861829/150‑1 (MQ=255)
                      gCCGCTGCTCTTGCGGGCAGGCATCCAGGCGCCGTTGCCATTCTTCCTTGGCCTTGGCCAGCTCATCCGCCGGGAACCCCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCg                                                                                                                          <  2:400757/150‑1 (MQ=255)
                       ccGCTGCTCTTGCGGGCAGGCATCCAGGCGCCGTTGCCATTCTTCCTTGGCCTTGGCCAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCgg                                                                                                                         <  1:869224/150‑1 (MQ=255)
                               ctTGCGGGCAGGCATCCAGGCGCCGTTGCCATTCTTCCTTGGCCTTGGCCAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTc                                                                                                                 <  1:819727/150‑1 (MQ=255)
                               ctTGCGGGCAGGCATCCAGGCGCCGTTGCCATTCTTCCTTGGCCTTGGCCAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTc                                                                                                                 >  1:445618/1‑150 (MQ=255)
                                 tGCGGGCAGGCATCCAGGCGCCGTTGCCATTCTTCCTTGGCCTTGGCCAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCgg                                                                                                               <  1:116601/150‑1 (MQ=255)
                                          gcaTCCAGGCGCCGTTGCCATTCTTCCTTGGCCTTGGCCAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGAAACAGTCGGTAGcccccc                                                                                                      >  2:220894/1‑150 (MQ=255)
                                           caTCCAGGCGCCGTTGCCATTCTTCCTTGGCCTTGGCCAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCa                                                                                                     <  2:656962/150‑1 (MQ=255)
                                           caTCCAGGCGCCGTTGCCATTCTTCCTTGGCCTTGGCCAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCa                                                                                                     <  2:445618/150‑1 (MQ=255)
                                            aTCCAGGCGCCGTTGCCATTCTTCCTTGGCCTTGGCCAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAg                                                                                                    >  1:1208549/1‑150 (MQ=255)
                                                     ccGTTGCCATTCTTCCTTGGCCTTGGCCAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCg                                                                                           <  1:668606/150‑1 (MQ=255)
                                                            cATTCTTCCTTGGCCTTGGCCAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTc                                                                                        >  1:1312715/1‑146 (MQ=255)
                                                            cATTCTTCCTTGGCCTTGGCCAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTc                                                                                        <  2:1312715/146‑1 (MQ=255)
                                                             aTTCTTCCTTGGCCTTGGCCAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAAtct                                                                                   >  2:652674/1‑150 (MQ=255)
                                                              ttcttcCTTGGCCCTGGCCAGCTCATCCGCCGGTAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAAtctc                                                                                  <  2:12706/150‑1 (MQ=255)
                                                                    cttggccttggcCAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATCTCTgccgc                                                                            >  1:188153/1‑150 (MQ=255)
                                                                      tggccttggccAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATCTCTGCCGCCa                                                                          <  1:1076642/150‑1 (MQ=255)
                                                                         ccttggccAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTATCCCTGCAGAATCTGCCGGTCAATCTCTGCCGCCAGCt                                                                       >  2:114488/1‑150 (MQ=255)
                                                                            tggccAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATCTCTGCCGCCAGCTGtt                                                                    <  2:243464/150‑1 (MQ=255)
                                                                             ggccAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCAACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTa                                                                                                             >  2:1313821/1‑108 (MQ=255)
                                                                             ggccAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCAACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTa                                                                                                             <  1:1314164/108‑1 (MQ=255)
                                                                             ggccAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCAACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTa                                                                                                             <  1:1313821/108‑1 (MQ=255)
                                                                             ggccAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCAACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTa                                                                                                             >  2:1314164/1‑108 (MQ=255)
                                                                               ccAGCTCATACGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATCTCTGCCGCCAGCTGTTTtg                                                                 >  1:56746/1‑150 (MQ=255)
                                                                                cAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATCTCTGCCGCCAGCTGTTTtgt                                                                <  2:1208549/150‑1 (MQ=255)
                                                                                  gCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCa                                                                                                                             <  2:1361132/87‑1 (MQ=255)
                                                                                  gCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCa                                                                                                                             >  1:1361132/1‑87 (MQ=255)
                                                                                   cTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATctct                                                                                 >  2:536388/1‑130 (MQ=255)
                                                                                   cTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATctct                                                                                 <  1:536388/130‑1 (MQ=255)
                                                                                    tCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATCTCTGCCGCCAGCTGTTTTGTGTCt                                                            >  2:1098938/1‑150 (MQ=255)
                                                                                        ccgccgGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATCTCTGCCGCCAGCTGTTTTGTGTCTTCGt                                                        <  1:572637/150‑1 (MQ=255)
                                                                                         cgccgGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATCTCTGCCGCCAGCTGTTTTGTGTCTTCGTg                                                       <  1:482069/150‑1 (MQ=255)
                                                                                          gccgGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGCCAATCTCTGCCGCCAGCTGTTTTGTGTCTTCGTgg                                                      <  1:114488/150‑1 (MQ=255)
                                                                                             gggAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATCTCTGCCGCCAGCTGTTTTGTGTCTTCGTGGTAc                                                   >  1:495460/1‑150 (MQ=255)
                                                                                             gggAACACCCTCTCGGCGCGTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATCTCTGCCGCCAGCTGTTTTGTGTCTTCGTGGTAc                                                   >  2:379211/1‑150 (MQ=255)
                                                                                                aaCACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATCTCTGCCGCCAGCTGTTTTGTGTCTTCGTGGTACTCg                                                <  2:56746/150‑1 (MQ=255)
                                                                                                aaCACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATCTCTGCCGCCAGCTGTTTTGTGTCTTCGTGGTACTCg                                                >  1:962843/1‑150 (MQ=255)
                                                                                                  cacCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATCTCTGCCGCCAGCTGTTTTGTGTCTTCGTGGTACTCGGt                                              >  1:1332458/1‑150 (MQ=255)
                                                                                                    cccTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATCTCTGCCGCCAGCTGTTTTGTGTCTTCGTGGTACTCGGTTa                                            >  1:1406110/1‑150 (MQ=255)
                                                                                                      ctctCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGt                                                                                                                  >  2:168038/1‑78 (MQ=255)
                                                                                                      ctctCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGt                                                                                                                  <  1:168038/78‑1 (MQ=255)
                                                                                                         tcGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCg                                                                                           <  2:1354888/98‑1 (MQ=255)
                                                                                                         tcGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCg                                                                                           >  1:1354888/1‑98 (MQ=255)
                                                                                                           ggCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATCTCTGCCGCCAGc                                                                        <  2:700249/115‑1 (MQ=255)
                                                                                                           ggCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATCTCTGCCGCCAGc                                                                        >  1:700249/1‑115 (MQ=255)
                                                                                                              gcTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATCTCTGCCGCCAGCTGTTTTGTGTCTTCGTGGTACTCGGTTACCGGCGgcgc                                  <  2:1332458/150‑1 (MQ=255)
                                                                                                                   gCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATCTCTGCCGCCAGCTGTTTTGTGTCTTCGTGGTACTCGGTTACCGGCGGCGCGAAGt                             >  2:436247/1‑150 (MQ=255)
                                                                                                                       cTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATCTCTGCCGCCAGCTGTTTTGTGTCTTCGTGGTACTCGGTTACCGGCGGCGCGAAGTTGAt                         <  1:831653/150‑1 (MQ=255)
                                                                                                                            gggCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATCTCTGCCGCCAGCTGTTTTGTGTCTTCGTGGTACTCGGTTACCGGCGGCGCGAAGTTGATGTGGa                    <  2:1406110/150‑1 (MQ=255)
                                                                                                                                  cGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATCTCTGCCGCCAGCTGTTTTGTGTCTTCGTGGTACTCGGTTACCGGCGGCGCGAAGTTGATGTGGACTCGGc              <  1:379211/150‑1 (MQ=255)
                                                                                                                                      cTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATCTCTGCCGCCAGCTGTTTTGTGTCTTCGTGGTACTCg                                                >  1:480642/1‑112 (MQ=255)
                                                                                                                                      cTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATCTCTGCCGCCAGCTGTTTTGTGTCTTCGTGGTACTCg                                                <  2:480642/112‑1 (MQ=255)
                                                                                                                                      cTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATCTCTGCCGCCAGCTGTTTTGTGTCTTCGTGGTACTCGGTTACCGGCGGCGCGAAGTTGATGTGGACTCGGCCCtt          <  2:495460/150‑1 (MQ=255)
                                                                                                                                       ttGCCCTCCCACATCGCATATGCTAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATCTCTGCCGCCAGCTGTTTTGTGTCTTCGTGGTACTCGGTTACCGGCGGCGCGAAGTTGATGTGGACTCGGCCCTTg         <  2:108843/150‑1 (MQ=255)
                                                                                                                                        tGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATCTCTGCCGCCAGCTGTTTTGTGTCTTCGTGGTACTCGGTTACCGGCGGCGCGAAGTTGATGTGGACTCGGCCCTTGt        <  1:652674/150‑1 (MQ=255)
                                                                                                                                           ccTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATCTCTGCCGCCAGCTGTTTTGTGTCTTCGTGGTACTCGGTTACCGGCGGCGCGAAGTTGATGTGGACTCGGCCCTTGTAGc     >  1:794967/1‑150 (MQ=255)
                                                                                                                                             tCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATCTCTGCCGCCAGCTGTTTTGTGTCTTCGTGGTACTCGGTTACCGGCGGCGCGAAGTTGATGTGGACTCGGCCCTTGTAGCCg   >  2:268325/1‑150 (MQ=255)
                                                                                                                                              cccACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATCTCTGCCGCCAGCTGTTTTGTGTCTTCGTGGTACTCGGTTACCGGCGGCGCGAAGTTGATGTGGACTCGGCCCTTGTAGCCgg  <  1:1098938/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GCATACTGCAACACCAGGTACGGCCGCTGCTCTTGCGGGCAGGCATCCAGGCGCCGTTGCCATTCTTCCTTGGCCTTGGCCAGCTCATCCGCCGGGAACACCCTCTCGGCGCTTGGCACCTGCAGGGCTTCGTCCTTGCCCTCCCACATCGCATATGCCAGGTAATGCACCGGGAACAGTCGGTAGCCCCCCAGAATCTGCCGGTCAATCTCTGCCGCCAGCTGTTTTGTGTCTTCGTGGTACTCGGTTACCGGCGGCGCGAAGTTGATGTGGACTCGGCCCTTGTAGCCGG  >  NC_002947/871945‑872236

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: