Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,697,446 T→C 7.2% S15G (AGC→GGC) ‡ PP_5595 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,697,4460TC 7.2% 236.0 / 6.6 83S15G (AGC→GGC) ‡PP_5595hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (39/38);  new base C (4/2);  total (43/40)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

TGTCCTTGCTTTCGTTGCTGCATGACCAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAGGATCAGGTAGGTCACCACTGCC  >  NC_002947/3697299‑3697592
                                                                                                                                                   |                                                                                                                                                  
tGTCCTTGCTTTCGTTGCTGCATGACCAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTg                                                                                                                                                   >  7:222367/1‑149 (MQ=255)
  tCCTTGCTTTCGTTGCTGCATGACCAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGc                                                                                                                                                 >  5:205895/1‑149 (MQ=255)
   ccTTGCTTTCGTTGCTGCATGACCAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGc                                                                                                                                                 >  4:168104/1‑148 (MQ=255)
    cTTGCTTTCGTTGCTGCATGACCAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAG‑CCTGGCCAGCGGCCAg                                                                                                                                              >  2:35796/1‑149 (MQ=255)
           tCGTTGCTGCATGACCAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAG‑CCTGGCCAGCGGCCAggtctg                                                                                                                                         <  1:35796/147‑5 (MQ=255)
             gTTGCTGCATGACCAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCaa                                                                                                                                      <  2:265566/149‑1 (MQ=255)
                gctgcATGACCAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCAGCGGCCAggt                                                                                                                                            >  8:158424/1‑139 (MQ=255)
                gctgcATGACCAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCAGCGGCCAggt                                                                                                                                            <  7:158424/140‑2 (MQ=255)
                 ctgcATGACCAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACca                                                                                                                                   <  7:405616/148‑1 (MQ=255)
                          cAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCatcatcatc                                                                                                                           >  5:354456/1‑147 (MQ=255)
                          cAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCatcatcatc                                                                                                                           >  2:228331/1‑147 (MQ=255)
                          cAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCatcatcatc                                                                                                                           >  2:163116/1‑147 (MQ=255)
                                aTCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAg                                                                                                                    <  5:272429/148‑1 (MQ=255)
                                aTCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAg                                                                                                                    <  2:280444/148‑1 (MQ=255)
                                 tCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAgg                                                                                                                   >  4:116623/1‑148 (MQ=255)
                                  cGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGaa                                                                                                                 >  3:200519/1‑149 (MQ=255)
                                   gCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGaa                                                                                                                 >  1:237393/1‑148 (MQ=255)
                                    cTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAaca                                                                                                               >  1:30114/1‑149 (MQ=255)
                                        ggcggGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGc                                                                                                            <  8:222367/148‑1 (MQ=255)
                                           gggTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCAttt                                                                                                        <  2:237393/149‑1 (MQ=255)
                                           gggTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCAttt                                                                                                        >  1:236710/1‑149 (MQ=255)
                                           gggTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCAttt                                                                                                        >  1:2030/1‑149 (MQ=255)
                                              tATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACcatcatc                                                                                                                              <  5:137547/124‑1 (MQ=255)
                                              tATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACcatcatc                                                                                                                              >  6:423813/1‑124 (MQ=255)
                                              tATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACcatcatc                                                                                                                              <  5:423813/124‑1 (MQ=255)
                                              tATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACcatcatc                                                                                                                              >  6:137547/1‑124 (MQ=255)
                                                 gTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGAAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCAcc                                                                                                  >  1:295398/1‑149 (MQ=255)
                                                    gctgtagctgtaACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTc                                                                                                <  4:200519/148‑1 (MQ=255)
                                                    gctgtagctgtaACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTc                                                                                                >  4:233287/1‑148 (MQ=255)
                                                      tgtagctgtaACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAg                                                                                             >  4:27815/1‑149 (MQ=255)
                                                       gtagctgtaACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAg                                                                                                                    <  7:14669/125‑1 (MQ=255)
                                                       gtagctgtaACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAg                                                                                                                    >  8:14669/1‑125 (MQ=255)
                                                        tagctgtaACAGGTGTGCCGCTCCACCGGCGTCACGTCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAgg                                                                                            <  2:295398/148‑1 (MQ=255)
                                                          gctgtaACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGca                                                                                         >  3:350445/1‑149 (MQ=255)
                                                            tgtaACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGcac                                                                                        >  1:405118/1‑148 (MQ=255)
                                                              taACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACaa                                                                                      <  1:163116/148‑1 (MQ=255)
                                                              taACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAg                                                                                     >  6:271293/1‑149 (MQ=255)
                                                              taACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGAGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAg                                                                                     >  4:21714/1‑149 (MQ=255)
                                                               aaCAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAg                                                                                     <  3:233287/148‑1 (MQ=255)
                                                                aCAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAgg                                                                                    <  7:250968/148‑1 (MQ=255)
                                                                aCAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGc                                                                                   >  8:282188/1‑149 (MQ=255)
                                                                aCAGGTGTGCCGCCCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGc                                                                                   >  7:157422/1‑149 (MQ=255)
                                                                  aGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGccc                                                                                 <  1:1331/149‑1 (MQ=255)
                                                                  aGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGccc                                                                                 <  3:168104/149‑1 (MQ=255)
                                                                      gtgCCGCTCAACCGGCGTCACGGCGTCGCTGGAGGAGCGCGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGGCAGGCTGGCACCCATGATCAGCAAGGCAGGCACCCGCAGATTCAACCTCCCGTGCACACGGCCCTCGt                                                                             >  8:317587/1‑149 (MQ=255)
                                                                       tgCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGt                                                                             <  2:30114/148‑1 (MQ=255)
                                                                        gCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGt                                                                             >  1:60959/1‑147 (MQ=255)
                                                                         ccGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAg                                                                          >  4:102695/1‑149 (MQ=255)
                                                                           gCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAgga                                                                        >  1:155245/1‑148 (MQ=255)
                                                                               cACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCt                                                                      <  4:209666/147‑1 (MQ=255)
                                                                               cACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCt                                                                      <  2:60959/147‑1 (MQ=255)
                                                                               cACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGa                                                                    <  8:157422/149‑1 (MQ=255)
                                                                                 ccGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGc                                                                  <  3:27815/149‑1 (MQ=255)
                                                                                    gCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAgg                                                               <  2:405118/149‑1 (MQ=255)
                                                                                     cGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGc                                                              <  3:21714/149‑1 (MQ=255)
                                                                                             cGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTcgcg                                                      <  1:228331/149‑1 (MQ=255)
                                                                                              gTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCgcgc                                                     <  4:350445/149‑1 (MQ=255)
                                                                                              gTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCgcgc                                                     >  4:413306/1‑149 (MQ=255)
                                                                                              gTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCgcgc                                                     <  5:271293/149‑1 (MQ=255)
                                                                                               tcccTGGAGGATCGAGTCATCACATTGCAAGCTTCGGCCCACACCTGCCCGCTGGCCAGGCTGGCACCAATCATCATCAAGGCAGGAACACGCAATTTCCAACTCCAGGGCACAAGGCCCTCGGAAGGACTGAGCAGGCTTGTCGCGcc                                                    >  2:385546/4‑149 (MQ=255)
                                                                                                     gaggagCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGgcgc                                               >  8:379219/1‑148 (MQ=255)
                                                                                                        gagCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATgg                                           <  7:282188/149‑1 (MQ=255)
                                                                                                                 aGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGgcg                                   <  3:116623/148‑1 (MQ=255)
                                                                                                                      aTTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCCAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGCGCAGGCTTGTCGCGCCGGCGCATGGCGTTGGCGCGtttt                             <  1:385546/149‑1 (MQ=255)
                                                                                                                          cAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTg                                      >  3:300327/1‑136 (MQ=255)
                                                                                                                          cAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTg                                      <  4:300327/136‑1 (MQ=255)
                                                                                                                          cAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCa                         <  1:404661/149‑1 (MQ=255)
                                                                                                                           aaGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTgg                                     <  2:135345/136‑1 (MQ=255)
                                                                                                                           aaGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTgg                                     >  1:135345/1‑136 (MQ=255)
                                                                                                                           aaGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAg                        >  4:429588/1‑149 (MQ=255)
                                                                                                                                tCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAGGAt                     >  1:104277/1‑147 (MQ=255)
                                                                                                                                  ggCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAGGATCag                  <  3:102695/148‑1 (MQ=255)
                                                                                                                                   gCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAGGATCagg                 >  8:210356/1‑148 (MQ=255)
                                                                                                                                     ccAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAGGATCaggtag              >  5:243346/1‑149 (MQ=255)
                                                                                                                                       cagccTGGCCAGCGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAGGATCaggtaggt            >  2:57836/4‑149 (MQ=255)
                                                                                                                                        gACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAGGATCAGGTAGGTc           <  6:354456/149‑1 (MQ=255)
                                                                                                                                         gccTGGCCAGCGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAGGATCAGGTAGGTca          >  6:387048/2‑149 (MQ=255)
                                                                                                                                         aCCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAGGATCAGGTAGGTca          >  6:187426/1‑149 (MQ=255)
                                                                                                                                              gCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATgg                                           >  6:371711/1‑111 (MQ=255)
                                                                                                                                              gCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATgg                                           <  5:371711/111‑1 (MQ=255)
                                                                                                                                              gCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAGGATCAGGTAGGTcacca       <  2:2030/147‑1 (MQ=255)
                                                                                                                                              gCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAGGATCAGGTAGGTCACCACt     <  7:379219/149‑1 (MQ=255)
                                                                                                                                                 gCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAGGATCAGGTAGGTCACCACTGcc  <  7:133875/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGTCCTTGCTTTCGTTGCTGCATGACCAGTTGATCGCTTCGGCGGGTATGTTGCTGTAGCTGTAACAGGTGTGCCGCTCCACCGGCGTCACGGCGTCGCTGGAGGAGCGAGTCAGCACATTGCAAGCTTCGGCCCAGACCTGGCCGCTGGCCAGGCTGGCAACCATCATCATCAAGGCAGGAACACGCATTTTCCACCTCCAGGGCACAAGGCCCTCGTCAGGTCTGAGCAGGCTTGTCGCGCCGGCGCATGGCTTTGGCGCGTTTTTCCAGGATCAGGTAGGTCACCACTGCC  >  NC_002947/3697299‑3697592

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑