Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,030,792 A→C 6.9% M364L (ATG→CTG)  PP_4432 → Xaa‑Pro aminopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,030,7920AC 6.9% 955.9 / 45.9 290M364L (ATG→CTG) PP_4432Xaa‑Pro aminopeptidase
Reads supporting (aligned to +/- strand):  ref base A (133/137);  new base C (12/8);  total (145/145)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

Only 100 of 320 total aligned reads displayed.

AGTTCGATATCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAG  >  NC_002947/5030643‑5030941
                                                                                                                                                     |                                                                                                                                                     
  ttCGATATCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTatggatg                                                                                                                                                     >  2:4624887/1‑150 (MQ=255)
   tCGATATCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTatggatgg                                                                                                                                                    <  2:6759618/150‑1 (MQ=255)
   tCGATATCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTatggatgg                                                                                                                                                    <  1:4402409/150‑1 (MQ=255)
     gATATCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATggatggat                                                                                                                                                  <  2:7529714/150‑1 (MQ=255)
        atCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGgatgat                                                                                                                                               >  1:6368792/1‑150 (MQ=255)
         tCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGAtt                                                                                                                                              <  2:7190787/150‑1 (MQ=255)
           cacaAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTgg                                                                                                                                            >  2:3415740/1‑150 (MQ=255)
            acaAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTggg                                                                                                                                           >  2:6096574/1‑150 (MQ=255)
             caAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTgggg                                                                                                                                          >  2:1081991/1‑150 (MQ=255)
                 gACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTc                                                                                                                                      <  1:4696656/150‑1 (MQ=255)
                 gACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTc                                                                                                                                      >  2:7923516/1‑150 (MQ=255)
                   cAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCga                                                                                                                                    <  2:4000735/150‑1 (MQ=255)
                   cAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCga                                                                                                                                    >  2:223865/1‑150 (MQ=255)
                          tGCGGTTACCCAATCGGTATCAGCTATCCACTAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTAcc                                                                                                                             <  2:6184044/150‑1 (MQ=255)
                            cGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCga                                                                                                                           <  1:3418729/150‑1 (MQ=255)
                            cGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCga                                                                                                                           <  1:3883184/150‑1 (MQ=255)
                            cGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCga                                                                                                                           <  1:3070635/150‑1 (MQ=255)
                              gTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCgaga                                                                                                                         <  1:1158012/150‑1 (MQ=255)
                              gTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGCATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCgaga                                                                                                                         <  1:8622103/150‑1 (MQ=255)
                                   ccAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATc                                                                                                                    >  2:6000289/1‑150 (MQ=255)
                                   ccAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATc                                                                                                                    >  1:5094827/1‑150 (MQ=255)
                                    cAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATcc                                                                                                                   >  2:4400639/1‑150 (MQ=255)
                                    cAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATcc                                                                                                                   <  1:3287875/150‑1 (MQ=255)
                                     aaTCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCt                                                                                                                  >  2:6960470/1‑150 (MQ=255)
                                     aaTCGGTATCAGCTATCCACAAGATTGGGGCGAGCGCACAATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGAACTTCCACTTCATGCCTGGCCAATGGATGGATGATTGGGGGCTCGAGATTAAAGAGAGCATCCt                                                                                                                  >  2:5059265/1‑150 (MQ=255)
                                       tCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCGTgt                                                                                                                >  1:5147180/1‑149 (MQ=255)
                                       tCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGa                                                                                                                >  1:574595/1‑150 (MQ=255)
                                        cGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGAt                                                                                                               <  1:2344498/150‑1 (MQ=255)
                                         ggTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATc                                                                                                              >  1:7671189/1‑150 (MQ=255)
                                           cctgtccTATCCGCCGGACTGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCAc                                                                                                            >  1:7004957/7‑150 (MQ=255)
                                             tCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTg                                                                                                          <  2:1879845/150‑1 (MQ=255)
                                               tccTATCCGCCGGACTGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCACCGaa                                                                                                        >  1:7433908/3‑149 (MQ=255)
                                                gCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATa                                                                                                       >  1:8508856/1‑150 (MQ=255)
                                                 cTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATAc                                                                                                      >  2:3040321/1‑150 (MQ=255)
                                                  tATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACt                                                                                                     >  2:4930929/1‑150 (MQ=255)
                                                   aTCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTg                                                                                                    >  2:544094/1‑150 (MQ=255)
                                                    tCCGCCGGACTGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCACCGAAACAgg                                                                                                   >  1:6695905/1‑150 (MQ=255)
                                                    tCCGCCGGACTGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCACCGAAACAgg                                                                                                   >  1:2537756/1‑150 (MQ=255)
                                                     ccaccaGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTggg                                                                                                  <  1:9126539/150‑1 (MQ=255)
                                                     ccaccaGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTggg                                                                                                  <  2:1284564/150‑1 (MQ=255)
                                                     ccaccaGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTggg                                                                                                  <  1:6329305/150‑1 (MQ=255)
                                                     ccaccaGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTggg                                                                                                  <  2:5604057/150‑1 (MQ=255)
                                                     ccaccaGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTggg                                                                                                  <  2:8109853/150‑1 (MQ=255)
                                                       accaGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGt                                                                                                >  1:8360337/1‑150 (MQ=255)
                                                        ccggactGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCACCGAAACAGGCGTc                                                                                               >  1:7199583/7‑150 (MQ=255)
                                                         cggactGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCACCGAAACAGGCGTCg                                                                                              >  1:1566088/6‑150 (MQ=255)
                                                         caGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCg                                                                                              <  1:8117381/150‑1 (MQ=255)
                                                            actGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCACCGAAACAGGCGTCgaat                                                                                           <  2:4149420/148‑3 (MQ=255)
                                                              tGGGGCGAACGCACCATGAGCCTGCGTCCCGGCGACCGCACGGAATTGCAGCCCGGCATGACCTTTCATTTCATGACAGGTCTCTGGCTGGAAACCATGGGGCTCGAGATTACCGAGAGCATCGTGATCACCGAAACAGGCGTCgaatgc                                                                                         <  1:1656029/150‑5 (MQ=255)
                                                                  gCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCt                                                                                     >  1:548101/1‑150 (MQ=255)
                                                                         cACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGAc                                                                              <  2:7332081/150‑1 (MQ=255)
                                                                         cACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGAc                                                                              <  2:8854860/150‑1 (MQ=255)
                                                                         cACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGAc                                                                              <  2:1142458/150‑1 (MQ=255)
                                                                          aCCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGAct                                                                             >  2:5276929/1‑150 (MQ=255)
                                                                          aCCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGAct                                                                             >  2:9098080/1‑150 (MQ=255)
                                                                          aCCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGAct                                                                             >  1:8214616/1‑150 (MQ=255)
                                                                             aTGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGActcct                                                                          >  1:8707244/1‑150 (MQ=255)
                                                                                  ccTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCaa                                                                     <  1:7383089/150‑1 (MQ=255)
                                                                                    tGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGc                                                                   >  1:2704373/1‑150 (MQ=255)
                                                                                    tGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGc                                                                   <  2:2636366/150‑1 (MQ=255)
                                                                                         ccGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTc                                                              <  2:1710194/150‑1 (MQ=255)
                                                                                         ccGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTc                                                              <  1:5407942/150‑1 (MQ=255)
                                                                                         ccGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTc                                                              <  2:6876233/150‑1 (MQ=255)
                                                                                         ccGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTc                                                              >  2:3782811/1‑150 (MQ=255)
                                                                                         ccGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTc                                                              <  2:2352279/150‑1 (MQ=255)
                                                                                         ccGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTc                                                              <  2:2139831/150‑1 (MQ=255)
                                                                                           gAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGt                                                            <  1:6456449/150‑1 (MQ=255)
                                                                                              cGACAACAGCGTGCTAAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCaa                                                         <  1:2822205/150‑1 (MQ=255)
                                                                                                    caGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTg                                                   >  2:5747733/1‑150 (MQ=255)
                                                                                                     aGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGa                                                  <  1:891958/150‑1 (MQ=255)
                                                                                                     aGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGa                                                  >  2:5860046/1‑150 (MQ=255)
                                                                                                        gTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGa                                               <  1:4211957/150‑1 (MQ=255)
                                                                                                        gTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGa                                               <  1:4475404/150‑1 (MQ=255)
                                                                                                            tCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTc                                                                                                                                      <  2:8371466/59‑1 (MQ=255)
                                                                                                            tCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAg                                           >  2:3321255/1‑150 (MQ=255)
                                                                                                              aaaCCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGtt                                         >  2:8634492/1‑150 (MQ=255)
                                                                                                              aaaCCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGtt                                         >  1:6899836/1‑150 (MQ=255)
                                                                                                               aaCCTGGTATGACCTTCCACTTCCTGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGNCTCCTCGCAAACTTTTCGTCAGTCCTTGACGAAAAGctt                                        <  2:6366832/150‑1 (MQ=255)
                                                                                                               aaCCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGttt                                        <  2:2534896/150‑1 (MQ=255)
                                                                                                                 ccTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGc                                      >  2:177869/1‑150 (MQ=255)
                                                                                                                     gTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGt                                  >  2:2735569/1‑150 (MQ=255)
                                                                                                                     gTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCAATGATACTGGGGTCTAGACGATCTGCTCGACTCAACGCAAGAATTTCGTCAAGGACTGAAGAAAAGTTAGAAAGt                                  >  1:1966295/1‑150 (MQ=255)
                                                                                                                      aaTGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTc                                 <  1:8933819/149‑1 (MQ=255)
                                                                                                                       aTGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGTTACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCg                                >  2:6556738/1‑150 (MQ=255)
                                                                                                                        tGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGc                               >  1:8028161/1‑150 (MQ=255)
                                                                                                                          aCCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTg                             <  1:7093033/150‑1 (MQ=255)
                                                                                                                            cTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGaaa                                             <  2:6671096/132‑1 (MQ=255)
                                                                                                                            cTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGa                           >  1:5328171/1‑150 (MQ=255)
                                                                                                                            cTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGa                           <  1:3664484/150‑1 (MQ=255)
                                                                                                                            cTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGa                           <  1:7332921/150‑1 (MQ=255)
                                                                                                                            cTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGa                           >  1:6754432/1‑150 (MQ=255)
                                                                                                                               ccACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGtt                        >  2:232283/1‑150 (MQ=255)
                                                                                                                                  cTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTtata                     <  1:5885918/150‑1 (MQ=255)
                                                                                                                                    tCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGt                   >  1:8175295/1‑150 (MQ=255)
                                                                                                                                     cATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTC‑‑CAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGcc                <  1:2967253/150‑1 (MQ=255)
                                                                                                                                           tGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACgg            >  1:1761254/1‑150 (MQ=255)
                                                                                                                                               cTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTAt        >  2:438390/1‑150 (MQ=255)
                                                                                                                                                tATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATc       >  1:7916474/1‑150 (MQ=255)
                                                                                                                                                 atggatggATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCt      >  1:2752882/1‑150 (MQ=255)
                                                                                                                                                    gatggatgATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCa   >  2:4503142/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGTTCGATATCCACAAGGACAGCCGTTGCGGTTACCCAATCGGTATCAGCTATCCACCAGATTGGGGCGAGCGCACCATGAGCCTGCGCCCGAGCGACAACAGCGTGCTCAAACCTGGTATGACCTTCCACTTCATGCCTGGCCTATGGATGGATGATTGGGGGCTCGAGATTACCGAGAGCATCCTGATCACTGATACTGGGGTCGAGACGTTCTGCTCGACTCCTCGCAAGCTTTTCGTCAAGGACTGACGAAAAGTTTGCAAGTCGCTGGAGTTATAGTGCCACGGCTATCTTCAG  >  NC_002947/5030643‑5030941

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: