Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,228,341 C→T 5.3% H379H (CAC→CAT) ‡ PP_2828 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,228,3410CT 5.3% 282.0 / 3.1 93H379H (CAC→CAT) ‡PP_2828hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (46/42);  new base T (2/3);  total (48/45)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

GTGGCAGCAAGCAGCGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCATCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCAGTGGCTTGCCGTTGGGGCAACTGCGGCAAATTTATGATGGGCAGTTGCCGTTGGCGG  >  NC_002947/3228202‑3228489
                                                                                                                                            |                                                                                                                                                    
gTGGCAGCAAGCAGCGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCATCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCGTGCCGACactg                                                                                                                                            <  1:546059/150‑4 (MQ=255)
 tGGCAGCAAGCAGCGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCATCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACAtggt                                                                                                                                           <  2:1538938/150‑1 (MQ=255)
  ggCAGCAAGCAGCGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCATCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACAtggtg                                                                                                                                          <  2:189212/150‑1 (MQ=255)
     agcaagcaGCGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCATCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACAtggtggtg                                                                                                                                       >  1:375292/1‑150 (MQ=255)
     agcaagcaGCGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCATCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACAtggtggtg                                                                                                                                       <  2:281013/150‑1 (MQ=255)
     agcaagcaGCGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCATCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACAtggtggtg                                                                                                                                       >  2:683725/1‑150 (MQ=255)
     agcaagcaGCGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCATCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACAtggtggtg                                                                                                                                       >  2:443652/1‑150 (MQ=255)
      gcaagcaGCGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCATCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCGTGCCGACactg                                                                                                                                            <  2:597234/144‑4 (MQ=255)
      gcaagcaGCGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCATCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCGTGCCGACactg                                                                                                                                            >  1:597234/1‑141 (MQ=255)
       caagcaGCGTGACGCCCTGCCGTTGCTATTTGGCAAAGCCTACGGCATCCACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGc                                                                                                                                     <  1:217519/150‑1 (MQ=255)
          gcagcGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCATCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGtt                                                                                                                                  >  1:1461842/1‑150 (MQ=255)
          gaagcGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCATCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGtt                                                                                                                                  >  1:74302/3‑150 (MQ=255)
           acgagTGAAGCCGTGTGGTCGGTATTTGGGAAAGCCGACGGCATCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTg                                                                                                                                 <  2:417945/146‑1 (MQ=255)
            agcGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCATCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGc                                                                                                                                <  2:1461842/150‑1 (MQ=255)
            agcGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCATCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGc                                                                                                                                >  2:305696/1‑150 (MQ=255)
            agcGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCATCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGc                                                                                                                                <  1:166703/150‑1 (MQ=255)
            agcGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCATCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGc                                                                                                                                <  2:1881312/150‑1 (MQ=255)
             gcGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCATCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGcc                                                                                                                               <  1:150472/150‑1 (MQ=255)
                tGACGCAGTGCCGTTGCTATTTGGCAAAGCCTACGGCATCAACCAGGCCAAACCCGCGTGAGTAAACGATCA‑CGCTTGTGCCAGGCAGCTTAAACACCTAGCGAGGTTCTTGCCGTGGCGTCTCGCCGACATGGTGGTGGCGTTGCCGTa                                                                                                                            <  1:443652/150‑1 (MQ=255)
                   cGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCATCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTTCCGTACCt                                                                                                                         <  2:573610/150‑1 (MQ=255)
                        tgccgtTGCTATTTGGCAAAGCCTACGGCATCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGcaca                                                                                                                    <  2:816166/150‑1 (MQ=255)
                                 tATTTGGCAAAGCCTACGGCATCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGGTGCCGTACCTGCACAGGGCCTTGc                                                                                                           >  2:1122328/1‑150 (MQ=255)
                                   tttGGCAAAGCCTACGGCATCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGa                                                                                                                                                 >  2:49275/1‑110 (MQ=255)
                                   tttGGCAAAGCCTACGGCATCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGa                                                                                                                                                 <  1:49275/110‑1 (MQ=255)
                                       gCAAAGCCTACGGCATCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCa                                                                                                     >  2:1705196/1‑150 (MQ=255)
                                        cAAAGCCTACGGCATCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGcaca                                                                                                                    <  1:74374/134‑1 (MQ=255)
                                        cAAAGCCTACGGCATCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGcaca                                                                                                                    >  2:74374/1‑134 (MQ=255)
                                           aGCCTACGGCATCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAAcc                                                                                                 <  2:629671/150‑1 (MQ=255)
                                             ccTACGGCATCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACGCt                                                                                               >  1:545643/1‑150 (MQ=255)
                                              cTACGGCATCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTg                                                                                              >  1:1033851/1‑150 (MQ=255)
                                                  ggCATCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGcccc                                                                                          >  1:1296819/1‑150 (MQ=255)
                                                  ggCATCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCCACCTGGTGGTGCCGTTGCCGTACCTGCACCGGGCCTCGCAAGGCAACCCCTGcccc                                                                                          >  2:831396/1‑150 (MQ=255)
                                                     aTCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGAt                                                                                       >  1:605718/1‑150 (MQ=255)
                                                          ccAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGTGCCTTGCAAGGCAAACCCTGCCCCGATGCAGt                                                                                  >  2:941833/1‑150 (MQ=255)
                                                           cAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGCTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGGTGCAGtg                                                                                 >  2:1334007/1‑149 (MQ=255)
                                                              gCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCaaa                                                                              >  1:1673868/1‑150 (MQ=255)
                                                                ccAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTAc                                                                                                                           >  1:593058/1‑103 (MQ=255)
                                                                ccAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTAc                                                                                                                           <  2:593058/103‑1 (MQ=255)
                                                                       gtgtGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGAt                                                                                       <  1:1760364/132‑1 (MQ=255)
                                                                       gtgtGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGAt                                                                                       >  2:1760364/1‑132 (MQ=255)
                                                                       gtgtGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCctggctg                                                                     >  2:1900935/1‑150 (MQ=255)
                                                                        tgtgCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCctggctgg                                                                    >  1:121004/1‑150 (MQ=255)
                                                                        tgtgCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCctggctgg                                                                    <  1:1705196/150‑1 (MQ=255)
                                                                         gtgCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGt                                                                   >  1:492597/1‑150 (MQ=255)
                                                                          tgCCT‑GCTGGTCAGCCCGTTTGACAGGCAGGGGAAAACCCTCGAGCGGTTCTTGCAGTGTCGGCCCGCCGACATGGTGGTGGCGTTGCCGTACCTGCACATGGCCTCGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCAGGCTGGTg                                                                  <  2:229982/150‑1 (MQ=255)
                                                                           gCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGa                                                                 <  2:121004/150‑1 (MQ=255)
                                                                           gCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCCGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGa                                                                 >  2:1233206/1‑150 (MQ=255)
                                                                            ccT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGaa                                                                >  1:789262/1‑150 (MQ=255)
                                                                                gctgTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTTCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATAc                                                             >  1:1445949/1‑150 (MQ=255)
                                                                                gctgTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATAc                                                             >  2:1817355/1‑150 (MQ=255)
                                                                                gctgTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATAc                                                             <  1:249422/150‑1 (MQ=255)
                                                                                gctgTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTACCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATAc                                                             >  2:233921/1‑150 (MQ=255)
                                                                                    ttCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCAGt                                                         >  2:978503/1‑150 (MQ=255)
                                                                                     tCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCAGTg                                                        <  1:1900935/150‑1 (MQ=255)
                                                                                       aGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCaa                                                                                                         >  1:1033667/1‑99 (MQ=255)
                                                                                       aGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCaa                                                                                                         <  2:1033667/99‑1 (MQ=255)
                                                                                        gCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCAGTGGCt                                                     >  2:1862593/1‑150 (MQ=255)
                                                                                         cccGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGATGGTGAATACCAGTGGCtt                                                    <  2:375292/150‑1 (MQ=255)
                                                                                          ccGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCt                                                                                                              >  1:819789/1‑91 (MQ=255)
                                                                                          ccGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCt                                                                                                              <  2:819789/91‑1 (MQ=255)
                                                                                          ccGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCAGTGGCTTg                                                   <  2:605718/150‑1 (MQ=255)
                                                                                            gTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCAGTGGCTTGcc                                                 <  2:1296819/150‑1 (MQ=255)
                                                                                                  caggcaggGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCAGTGGCTTGCCGTTggg                                           <  2:1749177/150‑1 (MQ=255)
                                                                                                  caggcaggGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCAGTGGCTTGCCGTTggg                                           <  2:492597/150‑1 (MQ=255)
                                                                                                    ggcaggGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCAGTGGCTTGCCGTTGGGGc                                         >  1:1941117/1‑150 (MQ=255)
                                                                                                       aggGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCAGTGGCTTGCCGTTGGGGCAAc                                      >  1:1323920/1‑150 (MQ=255)
                                                                                                        ggggAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCAGTGGCTTGCCGTTGGGGCAACt                                     >  1:373000/1‑150 (MQ=255)
                                                                                                         gggAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCAGTGGCTTGCCGTTGGGGCAACTg                                    >  1:1243797/1‑150 (MQ=255)
                                                                                                          ggAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCAGTGGCTTGCCGTTGGGGCAACTGc                                   <  2:789262/150‑1 (MQ=255)
                                                                                                           gAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCAGTGGCTTGCCGTTGGGGCAACTGCg                                  >  1:1056786/1‑150 (MQ=255)
                                                                                                            aaaaaCCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCAGTGGCTTGCCGTTGGGGCAACTGCgg                                 <  2:373000/150‑1 (MQ=255)
                                                                                                             aaaaCCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCAGTGGCTTGCCGTTGGGGCAACTGCGGc                                <  2:1323920/150‑1 (MQ=255)
                                                                                                              aaaCCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGcc                                                                                                                               <  1:1164225/54‑1 (MQ=255)
                                                                                                              aaaCCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGcc                                                                                                                               >  2:1164225/1‑54 (MQ=255)
                                                                                                               aaCCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGa                                                                                        >  2:1459659/1‑92 (MQ=255)
                                                                                                               aaCCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGa                                                                                        <  1:1459659/92‑1 (MQ=255)
                                                                                                                aCCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCAGTGGCTTGCCGTTGGGGCAACTGCGGCaaa                             <  1:1233206/150‑1 (MQ=255)
                                                                                                                   tCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCAGTGGCTTGCCGTTGGGGCAACTGCGGCAAAttt                          >  1:617813/1‑150 (MQ=255)
                                                                                                                   tCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCAGTGGCTTGCCGTTGGGGCAACTGCGGCAAAttt                          >  1:1229190/1‑150 (MQ=255)
                                                                                                                    cGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCAGTGGCTTGCCGTTGGGGCAACTGCGGCAAATTTa                         <  2:1673868/150‑1 (MQ=255)
                                                                                                                       gCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCCGTGGCGTGCCGTTGGGGCAACTGCGGCAAATTTatga                      <  1:941833/150‑1 (MQ=255)
                                                                                                                          gTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATATCAGTGGCTTGCCGTTGGGGAAACCGCGGCAATTCTATGAggg                   <  2:74302/150‑1 (MQ=255)
                                                                                                                          gTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCAGTGGCTTGCCGTTGGGGCAACTGCGGCAAAt                            >  2:1922577/1‑141 (MQ=255)
                                                                                                                          gTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCAGTGGCTTGCCGTTGGGGCAACTGCGGCAAAt                            <  1:1922577/141‑1 (MQ=255)
                                                                                                                             cTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCAGTGGCTTGCCGTTGGGGCAACTGCGGCAAATTTATGATGGGCa                <  1:1862593/150‑1 (MQ=255)
                                                                                                                               caccatgtCGGCGTGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCAGTGGCTTGCCGTTGGGGCAACTGCGGCAAATTTatga                      >  2:1786397/6‑142 (MQ=255)
                                                                                                                               caccatgtCGGCGTGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCAGTGGCTTGCCGTTGGGGCAACTGCGGCAAATTTatga                      <  1:1786397/137‑1 (MQ=255)
                                                                                                                                gCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCAGTGGCTTGCCGTTGGGGCAACTGCGGCAAATTTATGATGGGCAGtt             >  2:1662689/1‑150 (MQ=255)
                                                                                                                                 cAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCAGTGGCTTGCCGTTGGGGCAACTg                                    >  2:1067782/1‑126 (MQ=255)
                                                                                                                                 cAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCAGTGGCTTGCCGTTGGGGCAACTg                                    <  1:1067782/126‑1 (MQ=255)
                                                                                                                                        ggCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCAGTGGCTTGCCGTTGGGGCAACTGCGGCAAATTTATGATGGGCAGTTGCCGTTgg     >  2:1601472/1‑150 (MQ=255)
                                                                                                                                           aCGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCAGTGGCTTGCCGTTGGGGCAACTGCGGCAAATTTATGATGGGCAGTTGCCGTTGGCgg  >  2:1132003/1‑150 (MQ=255)
                                                                                                                                           aCGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCAGTGGCTTGCCGTTGGGGCAACTGCGGCAAATTTATGATGGGCAGTTGCCGTTGGCgg  <  1:978503/150‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
GTGGCAGCAAGCAGCGTGACGCCGTGCCGTTGCTATTTGGCAAAGCCTACGGCATCAACCAGGCCCAACTGGTGTGCCT‑GCTGTTCAGCCCGTTTGCCAGGCAGGGGAAAAACCTCGAGCGGTTCTTGCAGTGTCGGCACGCCGACATGGTGGTGGCGTTGCCGTACCTGCACAGGGCCTTGCAAGGCAAACCCTGCCCCGATGCAGTCAAAGCCTGGCTGGTGAATACCAGTGGCTTGCCGTTGGGGCAACTGCGGCAAATTTATGATGGGCAGTTGCCGTTGGCGG  >  NC_002947/3228202‑3228489

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: ‑