Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,803,350 G→C 7.1% T844T (ACC→ACG) ‡ PP_3361 ← aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,803,3500GC 7.1% 601.1 / 27.2 183T844T (ACC→ACG) ‡PP_3361aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (92/78);  new base C (2/11);  total (94/89)
Fisher's exact test for biased strand distribution p-value = 8.54e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 192 total aligned reads displayed.

CATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTGCCCCGCCACGCTGGGGA  >  NC_002947/3803201‑3803498
                                                                                                                                                      |                                                                                                                                                    
cATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCCTGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTgg                                                                                                                                                      <  1:4088603/150‑1 (MQ=255)
cATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTgg                                                                                                                                                      <  1:2677806/150‑1 (MQ=255)
 aTGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTggg                                                                                                                                                     >  2:4587025/1‑150 (MQ=255)
 aTGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTggg                                                                                                                                                     >  2:4663135/1‑150 (MQ=255)
   gACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGtc                                                                                                                                                   >  1:3530070/1‑150 (MQ=255)
      acCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGTAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCTTCCGGGACCACGCCCTGGGTCTCg                                                                                                                                                >  2:3836236/1‑150 (MQ=255)
      acCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCg                                                                                                                                                >  2:3234878/1‑150 (MQ=255)
      acCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCg                                                                                                                                                >  1:2212589/1‑150 (MQ=255)
      acCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCg                                                                                                                                                >  1:1043398/1‑150 (MQ=255)
        cGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGaa                                                                                                                                              <  1:1882206/150‑1 (MQ=255)
            ttgccaccGAAGGTGTTGAAATATTCCATGCCATTGTTAAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGcc                                                                                                                                          <  1:3021213/143‑1 (MQ=255)
            tagccaccGAAGGTGTTGAAATATTCCATGCCATTGTTAAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGcc                                                                                                                                          <  1:3255659/143‑1 (MQ=255)
             tACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGccc                                                                                                                                         >  1:1457912/1‑150 (MQ=255)
               ccGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCccac                                                                                                                                       <  2:3898604/150‑1 (MQ=255)
               ccGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCccac                                                                                                                                       <  1:2823616/150‑1 (MQ=255)
                 ggCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCccacca                                                                                                                                     <  1:1512929/150‑1 (MQ=255)
                 ggCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCccacca                                                                                                                                     >  1:2575615/1‑150 (MQ=255)
                 accGAAGGTGTTGAAATATTCCATGCCATTGTTAAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGCccacca                                                                                                                                     <  1:2047916/148‑1 (MQ=255)
                     aaGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTgg                                                                                                                                 <  1:82950/150‑1 (MQ=255)
                       ggTGTTGAAATATTCCATGCCATTGTTAAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGCCCACCAGTGGCt                                                                                                                               <  2:1732259/150‑1 (MQ=255)
                       ggTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCt                                                                                                                               <  2:2653797/150‑1 (MQ=255)
                       ggTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCt                                                                                                                               <  1:1700108/150‑1 (MQ=255)
                          aTTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGcc                                                                                                                            >  1:717244/1‑150 (MQ=255)
                              aaaTATTCCATGCCATTGTTAAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGCCCACCAGTGGCTGCCTAcg                                                                                                                        <  2:61277/150‑2 (MQ=255)
                                    tCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCg                                                                                                                  >  1:336368/1‑150 (MQ=255)
                                      cATGCCGTTGTTGAAATCATCGGCAACCTCGCGGGTGTTCACTACCGCCTACATCTGGTGGGCTTTACCAATGTTTTTGCCCATGCTGTCGTCTTCCGGGTCCACTCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGat                                                                                                                <  1:1601852/150‑2 (MQ=255)
                                      cATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGaa                                                                                                                <  1:3477964/150‑1 (MQ=255)
                                      cATGCCATTGTTAAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGCCCACCAGTGGCTGCCTACGCGACCaaa                                                                                                                <  1:1505374/150‑1 (MQ=255)
                                      cATGCCATTGTTAAAGCTGTCTGCCACTTCCCGCGTTGTGACAACCGCAGACATTGGATGACCGTTACCAATGGGTTTGCCCATAGTGACGATGTCAGGGACAACATCCTGCAACTCGAAAGCCCACCAGTGGCTGCCTACGCGACCaaa                                                                                                                <  1:378132/150‑1 (MQ=255)
                                       aTGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAAc                                                                                                               <  1:4167095/150‑1 (MQ=255)
                                       aTGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCAGCCGAAc                                                                                                               <  2:1814305/150‑1 (MQ=255)
                                           cGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGa                                                                                                           >  2:3417207/1‑150 (MQ=255)
                                           cGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGGTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGa                                                                                                           >  2:2935253/1‑150 (MQ=255)
                                            gttgttGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGAc                                                                                                          >  1:4973080/1‑150 (MQ=255)
                                             ttgttgAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGAcc                                                                                                         >  2:1632493/1‑150 (MQ=255)
                                                  gAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAAc                                                                                                    <  1:4684134/150‑1 (MQ=255)
                                                         aTCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATc                                                                                              >  2:345634/1‑150 (MQ=255)
                                                          tCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCg                                                                                             <  1:2902676/150‑1 (MQ=255)
                                                          tCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCg                                                                                             >  2:3734697/1‑150 (MQ=255)
                                                            ggCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGc                                                                                           >  1:3044679/1‑150 (MQ=255)
                                                              cAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCa                                                                                         <  1:5087398/150‑1 (MQ=255)
                                                              cAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCa                                                                                         <  2:2448392/150‑1 (MQ=255)
                                                                   tCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGcac                                                                                    >  2:3470661/1‑150 (MQ=255)
                                                                    cgcgcgTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGcaca                                                                                   <  2:5007019/150‑1 (MQ=255)
                                                                       gcgTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGcc                                                                                <  2:1569631/150‑1 (MQ=255)
                                                                        cgTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGccc                                                                               <  1:1190677/150‑1 (MQ=255)
                                                                        cgTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGccc                                                                               <  1:4774573/150‑1 (MQ=255)
                                                                             tCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCtt                                                                          >  2:842192/1‑150 (MQ=255)
                                                                              cACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTc                                                                         <  1:1274266/150‑1 (MQ=255)
                                                                                 tACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGc                                                                      >  1:2689989/1‑150 (MQ=255)
                                                                                 tACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGc                                                                      >  1:4015837/1‑150 (MQ=255)
                                                                                           aTCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATg                                                            <  2:303239/150‑1 (MQ=255)
                                                                                            tCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATgg                                                           >  1:1753812/1‑150 (MQ=255)
                                                                                              gggTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCg                                                         >  1:3742253/1‑150 (MQ=255)
                                                                                               ggTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGt                                                        <  1:807732/150‑1 (MQ=255)
                                                                                                  ggCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTaca                                                     <  1:1248112/150‑1 (MQ=255)
                                                                                                   gCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTacac                                                    >  2:4576257/1‑150 (MQ=255)
                                                                                                     cGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACAcct                                                  >  1:604736/1‑150 (MQ=255)
                                                                                                     cGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACAcct                                                  >  1:604768/1‑150 (MQ=255)
                                                                                                         aCCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCtt                                              <  2:2207175/150‑1 (MQ=255)
                                                                                                           cAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGa                                            <  2:119331/150‑1 (MQ=255)
                                                                                                           cAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGa                                            >  2:3705474/1‑150 (MQ=255)
                                                                                                           cAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGa                                            <  1:3188332/150‑1 (MQ=255)
                                                                                                            aaTGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAg                                           <  2:790338/150‑1 (MQ=255)
                                                                                                            aaTGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAg                                           <  2:5078093/150‑1 (MQ=255)
                                                                                                            aaTGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAg                                           <  1:1748596/150‑1 (MQ=255)
                                                                                                            aaTGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAg                                           <  2:3398924/150‑1 (MQ=255)
                                                                                                            aaTGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAg                                           <  1:1748504/150‑1 (MQ=255)
                                                                                                            aaTGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAg                                           >  1:973438/1‑150 (MQ=255)
                                                                                                             aTGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCccacca                                                                                                                                     <  2:4260704/59‑1 (MQ=255)
                                                                                                             aTGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAgg                                          >  1:2139912/1‑150 (MQ=255)
                                                                                                              tGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGt                                         <  1:4579391/150‑1 (MQ=255)
                                                                                                               gggTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTa                                        >  2:50998/1‑150 (MQ=255)
                                                                                                               gggTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTa                                        >  1:2641545/1‑150 (MQ=255)
                                                                                                                      cccATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCg                                 >  2:4505312/1‑150 (MQ=255)
                                                                                                                      cccATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCg                                 >  1:1082863/1‑150 (MQ=255)
                                                                                                                             tGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCaaaa                          <  2:1176466/150‑1 (MQ=255)
                                                                                                                               aCGGCAACCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAac                        >  1:3095269/1‑150 (MQ=255)
                                                                                                                                 ggCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAacac                      <  1:1632493/150‑1 (MQ=255)
                                                                                                                                 ggCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAacac                      >  1:3438380/1‑150 (MQ=255)
                                                                                                                                  gCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACAcc                     <  1:4420801/150‑1 (MQ=255)
                                                                                                                                  gCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACAcc                     >  1:4530063/1‑150 (MQ=255)
                                                                                                                                    aTCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTg                   >  2:905156/1‑150 (MQ=255)
                                                                                                                                    aTCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTg                   >  1:4223538/1‑150 (MQ=255)
                                                                                                                                    aTCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGACAAAAACACCTg                   >  2:4025418/1‑150 (MQ=255)
                                                                                                                                     tCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTGc                  >  2:473529/1‑150 (MQ=255)
                                                                                                                                       cGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTGccc                <  1:425489/150‑1 (MQ=255)
                                                                                                                                       cGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTGccc                <  2:268914/150‑1 (MQ=255)
                                                                                                                                       cGGGACAACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGAACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTGccc                >  2:1776081/1‑150 (MQ=255)
                                                                                                                                        gggACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTGcccc               >  2:2261349/1‑150 (MQ=255)
                                                                                                                                              aCGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTGCCCCGCCACg         >  2:3252087/1‑150 (MQ=255)
                                                                                                                                              aCGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTGGAGGTAATGTTCGGGCAAAAACACCTGCCCCGCCACg         >  2:2101884/1‑150 (MQ=255)
                                                                                                                                                gCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTGCCCCGCCACGCt       >  2:1367544/1‑150 (MQ=255)
                                                                                                                                                gCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTGCCCCGCCACGCt       >  2:2877512/1‑150 (MQ=255)
                                                                                                                                                 cccTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACTTTCGGCTCGCACCATGGCGTACACCTCCTCGAGGTAATGTTCGGGCAAAAACACCTGCCCCGCCACGCTg      <  1:4592070/150‑1 (MQ=255)
                                                                                                                                                 cccTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTGCCCCGCCACGCTg      <  2:5077800/150‑1 (MQ=255)
                                                                                                                                                 cccTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTGCCCCGCCACGCTg      <  2:2766207/150‑1 (MQ=255)
                                                                                                                                                  ccagggTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTGCCCCGCCACGCTgg     <  1:3836236/147‑1 (MQ=255)
                                                                                                                                                   cTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTACTGTTCGGGCAAAAACACCTGCCCCGCCACGCTggg    >  2:2735016/1‑150 (MQ=255)
                                                                                                                                                     gggTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTGCCCCGCCACGCTGGGGa  >  1:2337396/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CATGACACCGGGTTACCGGCGAAGGTATTGAAATACTCCATGCCGTTGTTGAAAGCATCGGCAACCTCGCGCGTGGTCACTACCGCCGACATCGGGTGGCCGTTACCAATGGGTTTGCCCATGCTGACGGCATCCGGGACCACGCCCTGGGTCTCGAATGCCCACCAGTGGCTGCCGACCCGGCCGAACCCGACCTGAACCTCATCGGCCAGGCACAGCCCACCTTCGGCTCGCACCATGGCGTACACCTCCTTGAGGTAATGTTCGGGCAAAAACACCTGCCCCGCCACGCTGGGGA  >  NC_002947/3803201‑3803498

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: