Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,883,651 C→G 9.3% A60G (GCC→GGC)  PP_5157 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,883,6510CG 9.3% 393.3 / 24.9 128A60G (GCC→GGC) PP_5157hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (56/60);  new base G (7/5);  total (63/65)
Fisher's exact test for biased strand distribution p-value = 5.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

Only 100 of 130 total aligned reads displayed.

GTACTGATCAGTCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGATGTAAGCCTG  >  NC_002947/5883503‑5883799
                                                                                                                                                        |                                                                                                                                                    
gTACTGATCAGTCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGcc                                                                                                                                                     >  1:2104268/1‑150 (MQ=255)
 tACTGATCAGTCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACTGGCCGAGAACTCGCCATACGCCTTCCAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGccc                                                                                                                                                    >  1:2340935/1‑150 (MQ=255)
   cTGATCAGTCCATTGGCTCAGGCGGAGCTGATCGACGAACTGGCGAGCCGGGGCGCACTGCGCATGGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTg                                                                                                                                                  >  2:84039/1‑150 (MQ=255)
   cTGATCAGTCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTg                                                                                                                                                  >  2:121229/1‑150 (MQ=255)
         aGTCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGGCCTGGCCaaa                                                                                                                                            <  2:1420800/150‑2 (MQ=255)
              aTTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACCCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACt                                                                                                                                       >  2:379268/1‑150 (MQ=255)
                   cTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGAtt                                                                                                                                  >  2:2376848/1‑150 (MQ=255)
                      aGGCGGAGCTGATCGACGAAATCGGCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTgc                                                                                                                               <  1:2141526/150‑1 (MQ=255)
                           gAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGcc                                                                                                                          <  1:845311/150‑1 (MQ=255)
                            aGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGTGCGCCg                                                                                                                         >  2:2350599/1‑150 (MQ=255)
                               tGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGt                                                                                                                      >  1:1745726/1‑150 (MQ=255)
                                caTCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAACCAAGACCACCACCTGACAGGGTTCGACATCGAATTT‑‑‑TGCCAGGCCCTGGCCAAGGAACTGGATCTGGGCGCCGCGtt                                                                                                                     >  2:2211005/2‑150 (MQ=255)
                                 aTCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTtc                                                                                                                    >  1:1174334/1‑150 (MQ=255)
                                    gacgaAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTtcgtc                                                                                                                 <  2:1745726/150‑1 (MQ=255)
                                    gacgaAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACCGGATTTGCGCGCCGAGTtcgtc                                                                                                                 <  2:1307773/150‑1 (MQ=255)
                                       gaAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGCGTTCGTCGaa                                                                                                              >  1:219319/1‑150 (MQ=255)
                                       gaAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGaa                                                                                                              <  2:242290/150‑1 (MQ=255)
                                        aaaTCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAgg                                                                                                                                           >  1:1611797/1‑120 (MQ=255)
                                        aaaTCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAg                                                                                                             <  2:1813938/150‑1 (MQ=255)
                                          aTCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGcc                                                                                                           >  2:1253330/1‑150 (MQ=255)
                                            cGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GACCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGcccc                                                                                                         <  1:235321/150‑1 (MQ=255)
                                              ccGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACAGCGAATTT‑‑‑GCGCCAGCCCTAGCGAATGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCgg                                                                                                       <  1:1253330/150‑1 (MQ=255)
                                                      ggCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGaa                                                                                               <  1:951862/150‑1 (MQ=255)
                                                       gCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAg                                                                                                                                            <  2:15325/104‑1 (MQ=255)
                                                       gCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAg                                                                                                                                            >  1:15325/1‑104 (MQ=255)
                                                          atcTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTga                                                                                           >  1:1488092/3‑149 (MQ=255)
                                                          aaCTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGc                                                                                           <  2:1590010/150‑1 (MQ=255)
                                                           aCTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCt                                                                                          <  1:796828/150‑1 (MQ=255)
                                                            cTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTg                                                                                         <  1:983173/150‑1 (MQ=255)
                                                            cTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTg                                                                                         >  1:2240339/1‑150 (MQ=255)
                                                             tGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTgg                                                                                        >  1:967762/1‑150 (MQ=255)
                                                                gcATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCgg                                                                                     <  2:2240339/150‑1 (MQ=255)
                                                                gcATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCgg                                                                                     <  2:958953/150‑1 (MQ=255)
                                                                gcATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCgg                                                                                     <  2:967762/150‑1 (MQ=255)
                                                                  aTTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGaa                                                                                               <  1:686313/138‑1 (MQ=255)
                                                                  aTTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCgggg                                                                                   >  1:1012402/1‑149 (MQ=255)
                                                                  aTTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGgt                                                                                   <  2:245458/150‑1 (MQ=255)
                                                                   ttGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCgggcg                                                                                  <  1:217379/150‑3 (MQ=255)
                                                                     gCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTc                                                                                <  1:121229/150‑1 (MQ=255)
                                                                      ccGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGcc                                                                                                                          >  2:796139/1‑107 (MQ=255)
                                                                      ccGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGcc                                                                                                                          <  1:796139/107‑1 (MQ=255)
                                                                           cAGGCCGAGACCTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGTCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGt                                                                          <  1:1704634/150‑1 (MQ=255)
                                                                            aGGCCGAGAACTCGCAAGACGCCTTCAAGCAAGACGACCACCTGACAGG‑‑TCTTCATTGAAACGCACGGCCTGCGGCTGGTTAGGAACTGGATCTGTGAGAAGAGTTCGTCGAAGCACCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTg                                                                         <  1:2350599/150‑1 (MQ=255)
                                                                                cgcgaACTCGCCAAACGCCTGCAAGCAAGAGGACCACCTGACAGGTTTCGCCATCGAATCT‑‑‑GGCCAGGCCCTGACCAAGGAACTCGATTTGCGCACCGACATGGCCGAAGCCCCGGCCACCGAAGTGCTAGCGGTAGTCGAGAGTGCCac                                                                     >  2:472366/4‑149 (MQ=255)
                                                                                cGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCaa                                                                     <  2:1328605/150‑1 (MQ=255)
                                                                                cGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCaa                                                                     >  2:2050322/1‑150 (MQ=255)
                                                                                cGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCa                                                                      >  2:1805643/1‑149 (MQ=255)
                                                                                cGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCa                                                                      <  1:1805643/149‑1 (MQ=255)
                                                                                 gagaACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTgcgc                                                                                                                             >  1:557934/1‑93 (MQ=255)
                                                                                 gagaACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTgcgc                                                                                                                             <  2:557934/93‑1 (MQ=255)
                                                                                   gaACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCTGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTa                                                                  <  2:2044055/150‑1 (MQ=255)
                                                                                      cTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCaa                                                                     >  1:1274705/1‑144 (MQ=255)
                                                                                      cTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGa                                                               >  1:1793631/1‑150 (MQ=255)
                                                                                       tCGCCATACGCCTTCTAGCAATACGACCATCTGACATGTTTGTACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTActtc                                                              <  1:84039/150‑4 (MQ=255)
                                                                                         gCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAATACGACAt                                                            >  1:723285/1‑150 (MQ=255)
                                                                                          ccATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACagc                                                           <  2:1241968/150‑3 (MQ=255)
                                                                                          ccATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATc                                                           <  1:472366/150‑1 (MQ=255)
                                                                                           cATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTCCGACATCg                                                          >  1:360202/1‑150 (MQ=255)
                                                                                            aTACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCga                                                                              <  1:1911797/129‑1 (MQ=255)
                                                                                            aTACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGc                                                         >  2:1539338/1‑150 (MQ=255)
                                                                                             tACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTg                                                                                                                                      <  1:1471311/72‑1 (MQ=255)
                                                                                             tACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTg                                                                                                                                      >  2:1471311/1‑72 (MQ=255)
                                                                                               cGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACa                                                             >  2:535228/1‑143 (MQ=255)
                                                                                               cGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACa                                                             <  1:535228/143‑1 (MQ=255)
                                                                                               cGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGc                                                         >  1:82761/1‑147 (MQ=255)
                                                                                               cGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGc                                                         <  2:82761/147‑1 (MQ=255)
                                                                                               cGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCAtt                                                      <  2:1012402/150‑1 (MQ=255)
                                                                                                        gcaagACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTc                                             <  1:2050322/150‑1 (MQ=255)
                                                                                                           agACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAg                                          <  2:1958708/150‑1 (MQ=255)
                                                                                                               gacCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGaa                                      >  1:1095564/1‑150 (MQ=255)
                                                                                                                 ccaccTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATc                                    <  1:1264110/150‑1 (MQ=255)
                                                                                                                     cTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTAcc                                                                                                  <  1:2262699/84‑1 (MQ=255)
                                                                                                                       gACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTc                                             >  2:792946/1‑135 (MQ=255)
                                                                                                                           ggTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGc                          <  2:2340935/150‑1 (MQ=255)
                                                                                                                            gTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCg                         >  2:607999/1‑150 (MQ=255)
                                                                                                                            gTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCg                         >  1:301685/1‑150 (MQ=255)
                                                                                                                                 gACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAATGCTGGCGGGTGTCGAGAGTGGCAAGTGCGGCAGCGCGTTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGc                    <  1:1855091/150‑1 (MQ=255)
                                                                                                                                 ctcATCGAATTT‑‑‑GGGCAGGCCATGGCCTTGGAATGATATGGGCGCCCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGc                    <  2:1833156/148‑1 (MQ=255)
                                                                                                                                    aTCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCg                 >  1:2106251/1‑150 (MQ=255)
                                                                                                                                    aTCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCCTGTTTACCCCGTCGAGTGAATCGCTCACAGGCGACGGccag                 >  1:1833156/1‑148 (MQ=255)
                                                                                                                                     tCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGa                                           <  1:1000051/123‑1 (MQ=255)
                                                                                                                                     tCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGa                                           >  2:1000051/1‑123 (MQ=255)
                                                                                                                                     tCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGc                >  2:1887380/1‑150 (MQ=255)
                                                                                                                                     tCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGc                >  2:586271/1‑150 (MQ=255)
                                                                                                                                       ttcctt‑‑‑GGCCAGGGCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTg                                                                         <  1:2213725/87‑1 (MQ=255)
                                                                                                                                       ttcctt‑‑‑GGCCAGGGCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTg                                                                         >  2:2213725/5‑91 (MQ=255)
                                                                                                                                       gAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCgaga                                                                            >  2:230670/1‑88 (MQ=255)
                                                                                                                                       gAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTTACAGGCGACGGCCCGCTc              >  1:340710/1‑150 (MQ=255)
                                                                                                                                         aTTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCgt            >  2:2315168/1‑149 (MQ=255)
                                                                                                                                           tt‑‑‑GGCCAGGGCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTg                                                                         >  1:1733945/1‑87 (MQ=255)
                                                                                                                                           tt‑‑‑GGCCAGGGCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTg                                                                         <  2:1733945/87‑1 (MQ=255)
                                                                                                                                           tt‑‑‑GGCCAGGGCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCa                                                                      <  1:534869/90‑1 (MQ=255)
                                                                                                                                           tt‑‑‑GGCCAGGGCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCATAGGCGACGGCCCGCTCGATg          >  1:1168472/1‑150 (MQ=255)
                                                                                                                                           tt‑‑‑GGCCAGGGCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGATg          >  2:419327/1‑150 (MQ=255)
                                                                                                                                           tt‑‑‑GGCCAGGGCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGATg          >  1:75681/1‑150 (MQ=255)
                                                                                                                                           tt‑‑‑GGCCAGGGCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGATg          >  1:1031322/1‑150 (MQ=255)
                                                                                                                                           tt‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGATg          >  1:1367426/1‑150 (MQ=255)
                                                                                                                                                  gCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTtcgtc                                                                                                                 <  1:1905561/44‑1 (MQ=255)
                                                                                                                                                  gCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTtcgtc                                                                                                                 >  2:1905561/1‑44 (MQ=255)
                                                                                                                                                       gCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGATGTAAGCCTg  >  1:226602/1‑150 (MQ=255)
                                                                                                                                                        |                                                                                                                                                    
GTACTGATCAGTCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGCCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTT‑‑‑GGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCAAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGATGTAAGCCTG  >  NC_002947/5883503‑5883799

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: