Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,039,746 A→C 7.5% P33P (CCT→CCG) ‡ PP_1812 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,039,7460AC 7.5% 222.8 / 10.7 80P33P (CCT→CCG) ‡PP_1812hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (38/36);  new base C (4/2);  total (42/38)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-01

CGTGCCAGTTTGTCGGGGGGAGGCTTACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTACCCGCTAGGGTTGGT  >  NC_002947/2039598‑2039894
                                                                                                                                                    |                                                                                                                                                    
cGTGCCAGTTTGTCGGGGGGAGGCTTACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAg                                                                                                                                                     <  1:1727995/150‑1 (MQ=255)
 tttgcAGTTTGTCGGGGGGAGGTTTACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAgg                                                                                                                                                    <  1:1249808/146‑1 (MQ=255)
 gTGCCAGTTTGTCGGGGGGAGGCTTACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAgg                                                                                                                                                    <  2:1880120/150‑1 (MQ=255)
 gTGCCAGTTTGTCGGGGGGAGGCTTACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAgg                                                                                                                                                    >  1:1705729/1‑150 (MQ=255)
    ccAGTTTGTCGGGGGGAGGCTTACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGcc                                                                                                                                                 >  2:530099/1‑150 (MQ=255)
        tttGTCGGGGGGAGGCTTACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTtgt                                                                                                                                             <  2:1515403/150‑1 (MQ=255)
        tttGTCGGGGGGAGGCTTACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTCGACACTGGGTATCGGACAATGGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTtgt                                                                                                                                             <  1:530099/150‑1 (MQ=255)
          tGTCGGGGGGAGGCTTACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGc                                                                                                                                           <  1:1613469/150‑1 (MQ=255)
                        ttACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAgggccgggc                                                                                                                             >  2:1279199/1‑150 (MQ=255)
                          aCCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGggccgggccg                                                                                                                           >  2:1690468/1‑150 (MQ=255)
                               cTGGGGGAGGGCGGTTTCCCGGGCGGGCCTTGTTCGGCGGGCGTTGGAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCgg                                                                                                                      <  1:603420/150‑1 (MQ=255)
                                tGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCggg                                                                                                                     >  2:841358/1‑150 (MQ=255)
                                  ggggAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGc                                                                                                                   <  1:1279199/150‑1 (MQ=255)
                                         tcGGTTTACCTGGCGGGCCTTTTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACg                                                                                                            <  2:787406/149‑1 (MQ=255)
                                          cGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACgg                                                                                                           <  2:829111/150‑1 (MQ=255)
                                                    ggcgggcCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTtg                                                                                                 <  2:1900505/150‑1 (MQ=255)
                                                       gggcCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTgtgg                                                                                              <  1:1280629/150‑1 (MQ=255)
                                                       gggcCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTgtgg                                                                                              >  1:839425/1‑150 (MQ=255)
                                                          ccTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTgtgggtg                                                                                           >  1:1436803/1‑150 (MQ=255)
                                                           cTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGc                                                                                                                                                  >  1:1871554/1‑94 (MQ=255)
                                                           cTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGc                                                                                                                                                  <  2:1871554/94‑1 (MQ=255)
                                                           cTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGACCGGGGCGGAAACGGCCACGGTTTgtgggtgg                                                                                          <  1:841358/150‑1 (MQ=255)
                                                               ttCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACACGCACAAGGCCCTGGCCTTGTGcctgt                                                                                                                                       <  1:1478507/101‑5 (MQ=38)
                                                               ttCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACACGCACAAGGCCCTGGCCTTGTGcctgt                                                                                                                                       >  2:1478507/1‑97 (MQ=38)
                                                                tCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGaa                                                                                     <  2:311430/150‑1 (MQ=255)
                                                                     gggcgTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTGCAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCTGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGa                                                                                >  2:1397686/1‑150 (MQ=255)
                                                                     gggcgTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGa                                                                                <  2:1705729/150‑1 (MQ=255)
                                                                         gTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGt                                                                            >  2:1757108/1‑150 (MQ=255)
                                                                             tAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAgggccgg                                                                                                                               <  2:815747/95‑1 (MQ=255)
                                                                             tAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAgggccgg                                                                                                                               >  1:815747/1‑95 (MQ=255)
                                                                                 agaATGGCCGGGCAGTGGTGGGCAGTTGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGa                                                                    <  2:451926/150‑1 (MQ=255)
                                                                                    aTGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCtt                                                                                                                                               <  2:400718/72‑1 (MQ=255)
                                                                                    aTGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCtt                                                                                                                                               >  1:400718/1‑72 (MQ=255)
                                                                                     tGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAgggccgg                                                                                                                               <  2:254490/87‑1 (MQ=255)
                                                                                     tGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAgggccgg                                                                                                                               >  1:254490/1‑87 (MQ=255)
                                                                                       gCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGAtt                                                              >  1:1453884/1‑150 (MQ=255)
                                                                                       gCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCGTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGAtt                                                              <  1:28765/150‑1 (MQ=255)
                                                                                           ggCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGt                                                          <  2:839425/150‑1 (MQ=255)
                                                                                              aGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGttt                                                                          >  2:309122/1‑131 (MQ=255)
                                                                                              aGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGttt                                                                          <  1:309122/131‑1 (MQ=255)
                                                                                               gtggtgGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAAAGGAAACGGTTTGTTGGTGGGGGAAGTTGAGCGTTTTGCTGACCGATTCGGTAAcc                                                      >  1:1098745/1‑150 (MQ=255)
                                                                                                tggtggGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCa                                                     >  2:721391/1‑150 (MQ=255)
                                                                                                  gtggGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAg                                                                                    >  1:1776610/1‑117 (MQ=255)
                                                                                                  gtggGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAg                                                                                    <  2:1776610/117‑1 (MQ=255)
                                                                                                              tATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGt                                       >  2:1472006/1‑150 (MQ=255)
                                                                                                              aaTCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGAGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGt                                       >  1:1451389/2‑150 (MQ=255)
                                                                                                                tCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTg                                     >  2:818274/1‑150 (MQ=255)
                                                                                                                tCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTg                                     <  1:818274/150‑1 (MQ=255)
                                                                                                                  ggACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTggg                                   >  1:1136712/1‑150 (MQ=255)
                                                                                                                   gACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGt                                  >  2:1418197/1‑150 (MQ=255)
                                                                                                                         cGGGGCTTTGGCTACCGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCCGGCGGAGCGGGGCGGAAACGGCAACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATtgtg                            <  2:611817/150‑1 (MQ=255)
                                                                                                                           gggATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCg                          >  1:683969/1‑150 (MQ=255)
                                                                                                                            ggATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCgg                         >  2:1336599/1‑150 (MQ=255)
                                                                                                                               tttGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCGTTCGCCGGGGTAGTTGGGTATTGTGCGGCCa                      >  1:1632131/1‑150 (MQ=255)
                                                                                                                               tttGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGccg                      <  2:1451389/150‑2 (MQ=255)
                                                                                                                                ttGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAg                     >  1:1747233/1‑150 (MQ=255)
                                                                                                                                   gCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGtt                                                                                  <  1:1906113/86‑1 (MQ=255)
                                                                                                                                   gCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGtt                                                                                  >  2:1906113/1‑86 (MQ=255)
                                                                                                                                   ccttaacgCACAAGGCCCTGGCCTTGTGCCTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGtt                  >  2:1644814/8‑150 (MQ=255)
                                                                                                                                   cctgaacgCACAAGGCCCTGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGtt                  >  1:878157/8‑150 (MQ=255)
                                                                                                                                    cTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGGATTTTGCGGCCAGGTTa                 >  1:663949/1‑150 (MQ=255)
                                                                                                                                     tACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGTGGAAACGGCCa                                                                                                        <  2:652775/62‑1 (MQ=255)
                                                                                                                                     tACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGTGGAAACGGCCa                                                                                                        >  1:652775/1‑62 (MQ=255)
                                                                                                                                      aCAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTAcc               <  1:1122122/150‑1 (MQ=255)
                                                                                                                                       cAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGCTCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTAccc              >  2:761271/1‑150 (MQ=255)
                                                                                                                                        acgCACAAGGCCCTGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCg                          <  2:927936/135‑1 (MQ=255)
                                                                                                                                        acgCACAAGGCCCTGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCg                          >  1:927936/3‑137 (MQ=255)
                                                                                                                                        aGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTgg                                    <  2:1358058/127‑1 (MQ=255)
                                                                                                                                        aGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTgg                                    >  1:1358058/1‑127 (MQ=255)
                                                                                                                                        aGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTACCCg             >  2:685189/1‑150 (MQ=255)
                                                                                                                                           cacaAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTACCCGCTa          <  1:685189/150‑1 (MQ=255)
                                                                                                                                           aacaAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGCGAAATTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTACCCGCTa          <  1:1397686/149‑1 (MQ=255)
                                                                                                                                            acaAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTACCCGCTAg         <  1:721391/150‑1 (MQ=255)
                                                                                                                                             caAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTa                                 <  2:489960/125‑1 (MQ=255)
                                                                                                                                             caAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTa                                 >  1:489960/1‑125 (MQ=255)
                                                                                                                                              aaGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTACCCGCTAggg       >  2:1734908/1‑150 (MQ=255)
                                                                                                                                               aGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTACCCGCTAGggt      >  2:800829/1‑150 (MQ=255)
                                                                                                                                                   cAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGAtt                                                              <  2:1265624/90‑1 (MQ=255)
                                                                                                                                                   cAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGAtt                                                              >  1:1265624/1‑90 (MQ=255)
                                                                                                                                                   cAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACAATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTACCCGCTAGggttggt  >  1:451315/1‑150 (MQ=255)
                                                                                                                                                    aGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAg                                                                               >  2:1424160/1‑72 (MQ=255)
                                                                                                                                                    aGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAg                                                                               <  1:1424160/72‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGTGCCAGTTTGTCGGGGGGAGGCTTACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTACCCGCTAGGGTTGGT  >  NC_002947/2039598‑2039894

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: