Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,039,746 A→C 12.8% 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,7460AC12.8% 208.1 / 24.4 86P33P (CCT→CCG) ‡PP_1812hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (34/41);  new base C (6/5);  total (40/46)
Fisher's exact test for biased strand distribution p-value = 7.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.04e-01

GTGCCAGTTTGTCGGGGGGAGGCTTACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTACCCGCTAGGGTTGGT  >  NC_002947/2039599‑2039894
                                                                                                                                                   |                                                                                                                                                    
gTGCCAGTTTGTCGGGGGGAGGCTTACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAgg                                                                                                                                                    >  1:1150281/1‑150 (MQ=255)
gTGCCAGTTTGTCGGGGGGAGGCTTACCGACTGGGGAATGGCTGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTAAAGGCACAAGGCCAgg                                                                                                                                                    <  2:2411631/150‑1 (MQ=255)
 ggCCAGTTTGTCGGGGGGAGGCTTACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAggg                                                                                                                                                   <  2:156686/149‑1 (MQ=255)
    cAGTTTGTCTGGGGGAGGCTTACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCt                                                                                                                                                >  2:2234664/1‑150 (MQ=255)
    cAGTTTGTCGGGGGGAGGCTTACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCt                                                                                                                                                <  1:2176738/150‑1 (MQ=255)
    cAGTTTGTCGGGGGGAGGCTTACCGACTGGGGGAGAGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGAGATTTGGCTACAGGCACAAGGCCAGGGCCt                                                                                                                                                >  2:1276605/1‑150 (MQ=255)
       tttGTCGGGGGGAGGCTTACCGACTGGGGGAGGGCGGTTTACCGGGTGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTtgt                                                                                                                                             <  1:736813/150‑1 (MQ=255)
       tttGTCGGGGGGAGGCTTACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGTCAGCGGGTATCGGATAAACGGGGATTTGGCTACAGGCACAAGGCCCTGGCCTtgt                                                                                                                                             <  1:2127232/150‑1 (MQ=255)
       tgtGTCGGGGGGAGGCTTACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTtgt                                                                                                                                             <  1:1424654/148‑1 (MQ=255)
              gggggAGGCTTACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCa                                                                                                                                      >  2:977119/1‑150 (MQ=255)
               ggggAGGCTTACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAg                                                                                                                                     <  1:977119/150‑1 (MQ=255)
                 ggAGGCTTACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTAAGGGATGCGTGGTAGGAGAATGGGCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAggg                                                                                                                                   >  1:1241215/1‑150 (MQ=255)
                   aGGCTTACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTtgt                                                                                                                                             >  2:462278/1‑138 (MQ=255)
                   aGGCTTACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTtgt                                                                                                                                             <  1:462278/138‑1 (MQ=255)
                     gCTTACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTtgt                                                                                                                                             >  2:2317163/1‑136 (MQ=255)
                     gCTTACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTtgt                                                                                                                                             <  1:2317163/136‑1 (MQ=255)
                     gCTTACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAgggccgg                                                                                                                               >  1:570136/1‑150 (MQ=255)
                         aCCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGggccgggccg                                                                                                                           >  1:1714175/1‑150 (MQ=255)
                            gACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCCTGGCCTTGTGcc                                                                                                                                          >  2:32607/1‑131 (MQ=255)
                            gACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCCTGGCCTTGTGcc                                                                                                                                          <  1:32607/132‑2 (MQ=255)
                               ggggggAGGGCGGTTCACCCGTGGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGGAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACATGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCTGAGCggg                                                                                                                     <  1:2417395/149‑1 (MQ=255)
                                 ggggAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGc                                                                                                                   >  2:751620/1‑150 (MQ=255)
                                 ggggAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGc                                                                                                                   >  1:1792095/1‑150 (MQ=255)
                                 ggggAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGc                                                                                                                   >  2:605899/1‑150 (MQ=255)
                                      gggCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGaaa                                                                                                              <  2:570136/150‑1 (MQ=255)
                                       ggCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAAc                                                                                                             >  1:1399032/1‑150 (MQ=255)
                                            tttACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCa                                                                                                        >  2:1418705/1‑150 (MQ=255)
                                               aCCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACgg                                                                                                     <  1:2340391/150‑1 (MQ=255)
                                                cggggcggggGTTGTTCGGCGGGCGTTGTAGGAGAATGGGCGGGCAGTGGTGGGCAGTGGGTAGCGGACAAACGGGGATTTGGCTACAGTCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGt                                                                                                    <  1:610335/148‑1 (MQ=255)
                                                ccGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGcc                                                                                                                                                 >  2:1487417/1‑105 (MQ=255)
                                                ccGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGcc                                                                                                                                                 <  1:1487417/105‑1 (MQ=255)
                                                ccGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGt                                                                                                    <  2:1150281/150‑1 (MQ=255)
                                                   ggcgggcCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGAACCGTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTtg                                                                                                 <  1:1636693/150‑1 (MQ=255)
                                                    gcgggcCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTtgt                                                                                                <  2:2132328/150‑1 (MQ=255)
                                                      gggcCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTgtgg                                                                                              <  2:1494233/150‑1 (MQ=255)
                                                      gggcCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTgtgg                                                                                              <  1:2234664/150‑1 (MQ=255)
                                                         ccTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTgtgggtg                                                                                           <  2:820474/150‑1 (MQ=255)
                                                          cTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTgtgggtgg                                                                                          <  1:751620/150‑1 (MQ=255)
                                                            tGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTgggtgggt                                                                                        <  1:129769/150‑1 (MQ=255)
                                                            ggTTCGGCGGGCGTTGGAGGAGAATGGCCGGGCGGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTgggtgggt                                                                                        <  2:701323/149‑1 (MQ=255)
                                                              ttCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGggtgggtgg                                                                                      >  1:2356376/1‑149 (MQ=255)
                                                                 ggcgggcgTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGt                                                                                   <  1:605899/150‑1 (MQ=255)
                                                                     ggcgTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAg                                                                               <  2:1714175/150‑1 (MQ=255)
                                                                      gcgTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGc                                                                              <  1:425691/150‑1 (MQ=255)
                                                                        gTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGt                                                                            <  2:225568/150‑1 (MQ=255)
                                                                             aGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTtgt                                                                                                                                             <  1:317311/80‑1 (MQ=255)
                                                                             aGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTtgt                                                                                                                                             >  2:317311/1‑80 (MQ=255)
                                                                              ggAGAATGGCCGGGCCGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCt                                                                      >  1:2264745/1‑150 (MQ=255)
                                                                               gagaATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTg                                                                     >  2:1361032/1‑150 (MQ=255)
                                                                                        cGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCg                                                            <  1:183832/150‑1 (MQ=255)
                                                                                        cGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCg                                                            <  1:1689323/150‑1 (MQ=255)
                                                                                           gCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGCCCGATTCGGTa                                                         <  2:1676561/150‑1 (MQ=255)
                                                                                             aGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAAc                                                       <  1:1891593/150‑1 (MQ=255)
                                                                                             aGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAAc                                                       >  2:2143731/1‑150 (MQ=255)
                                                                                              gtggtgGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAAcc                                                      <  1:1276605/150‑1 (MQ=255)
                                                                                                  tggGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTc                                                  <  1:463111/150‑1 (MQ=255)
                                                                                                       aGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAAAGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCgg                                             <  1:1478252/150‑1 (MQ=255)
                                                                                                               tCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGACGGGCGGGGCGGAAACGGACACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGGTGACCGATTCGGTAACCATTCGCCGGGGTAGTTg                                     <  2:802553/150‑1 (MQ=255)
                                                                                                                 ggACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTggg                                   >  1:797555/1‑150 (MQ=255)
                                                                                                                     aaaCGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTAtt                               >  2:393226/1‑150 (MQ=255)
                                                                                                                       aCGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTgtgggtgt                                                                                          <  1:462345/89‑2 (MQ=255)
                                                                                                                       aCGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTgtgggtgg                                                                                          >  2:462345/1‑89 (MQ=255)
                                                                                                                       aCGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATGCGCCGGGGTAGTTGGGTATtgt                             >  2:302280/1‑150 (MQ=255)
                                                                                                                               ttGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAg                     <  2:1399032/150‑1 (MQ=255)
                                                                                                                               ttGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAg                     <  1:302280/150‑1 (MQ=255)
                                                                                                                                 ggCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAg                                                                                    >  1:847070/1‑85 (MQ=255)
                                                                                                                                 ggCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAg                                                                                    <  2:847070/85‑1 (MQ=255)
                                                                                                                                    tACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTAc                >  2:1990894/1‑150 (MQ=255)
                                                                                                                                    tACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTAc                >  2:1768066/1‑150 (MQ=255)
                                                                                                                                    tACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGATGAGCGTTTTACTGAACGGGTCGGTAACCATTCGCCGGGGTTGTTGGGTATTGTGCGGCAAGGTTaa                >  2:561915/1‑149 (MQ=255)
                                                                                                                                       aGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTACCCg             >  2:1439011/1‑150 (MQ=255)
                                                                                                                                        cgCACAAGGCCCTGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTACCCGc            >  2:1363933/2‑150 (MQ=255)
                                                                                                                                         gCACAAGGCCCTGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGtt                                                                                  >  2:1980978/1‑79 (MQ=39)
                                                                                                                                         gCACAAGGCCCTGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGtt                                                                                  <  1:1980978/79‑1 (MQ=39)
                                                                                                                                         gCACAAGGCCCTGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAg                     >  1:1512735/1‑140 (MQ=255)
                                                                                                                                         gCACAAGGCCCTGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAg                     <  2:1512735/140‑1 (MQ=255)
                                                                                                                                         gCACAAGGCCCTGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTACCCGCt           >  1:2075035/1‑150 (MQ=255)
                                                                                                                                          cacaAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTACCCGCTa          >  1:443370/1‑150 (MQ=255)
                                                                                                                                           acaAGGCCCTGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGc                                                                                                                   >  1:877291/1‑44 (MQ=25)
                                                                                                                                           acaAGGCCCTGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGc                                                                                                                   <  2:877291/44‑1 (MQ=25)
                                                                                                                                           acaAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTACCCGCTAg         >  2:2235814/1‑150 (MQ=255)
                                                                                                                                               ggCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTACCCGCTAGggtt     <  2:2379927/150‑1 (MQ=255)
                                                                                                                                                gCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTACCCGCTAGggttg    >  2:772732/1‑150 (MQ=255)
                                                                                                                                                 ccAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGGAGTTGGGTATTGTGCGGCCAGGTTACCCGCTAGggttgg   >  1:2228922/1‑150 (MQ=255)
                                                                                                                                                  cAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTACCCGCTAGggttggt  >  2:690210/1‑150 (MQ=255)
                                                                                                                                                  cAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTACCCGCTAGggttggt  <  1:1418705/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GTGCCAGTTTGTCGGGGGGAGGCTTACCGACTGGGGGAGGGCGGTTTACCGGGCGGGCCTTGTTCGGCGGGCGTTGTAGGAGAATGGCCGGGCAGTGGTGGGCAGTGGGTATCGGACAAACGGGGATTTGGCTACAGGCACAAGGCCAGGGCCTTGTGCGTTCAGGGCCGGGCCGAGCGGGGCGGAAACGGCCACGGTTTGTGGGTGGGTGAAGTTGAGCGTTTTGCTGACCGATTCGGTAACCATTCGCCGGGGTAGTTGGGTATTGTGCGGCCAGGTTACCCGCTAGGGTTGGT  >  NC_002947/2039599‑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: