Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 898,152 C→T 9.9% Y11Y (TAC→TAT PP_0781 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947898,1520CT 9.9% 175.6 / 8.7 69Y11Y (TAC→TATPP_0781hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (34/28);  new base T (4/3);  total (38/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAACAACTGAACAGCTGAACTTCGCCGAGCATGCCGGAAACGAAGTCGCATTTTCCTACAACTGACTCAGAAACCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCT  >  NC_002947/898005‑898300
                                                                                                                                                   |                                                                                                                                                    
gCAACAACTGAACAGCTGAACTTCGCCGAGCATGCCGGAAACGAAGTCGCATTTTCCTACAACTGACTCAGAAACCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCt                                                                                                                                                    <  2:1304791/150‑1 (MQ=255)
     aacTGAACAGCTGAACTTCGCCGAGCATGCCGGAAACGAAGTCGCATTTTCCTAACAGTGACTCAGAAACCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGc                                                                                                                                               >  1:610925/1‑150 (MQ=255)
      acTGAACAGCTGAACTTCGCCGAGCATGCCGGAAACGAAGTCGCATTTTCCTACAACTGACTAAGAAATCACTGATCGGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGcc                                                                                                                                              <  1:196866/150‑1 (MQ=17)
       cTGAACAGCTGAACTTCGCCGAGCATGCCGGAAACGAAGTCGCATTTTCCTACAACTGACTCAGAAACCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGccc                                                                                                                                             >  1:475823/1‑150 (MQ=255)
           aCAGCTGAACTTCGCCGAGCATGCCGGAAACGAAGTCGCATTTTCCTACAACTGACTCAGAAACCTCCGATCCGCCTCAGTGATGGCATTTAACCTCTCGCAGCTATGCATGTAGTCACCTCATGTACCTGGATTACCTCGTGCcctcct                                                                                                                                         <  1:1265278/150‑1 (MQ=16)
              gcctggtagaaggctgagcgcacagcgtcgatccggcgcaggttgaggcccagcgaggcgacttcttcattgaggctggtgcggaactggccaccgcttgggtagcccaccatggcgatgacacgttactcgATCTCGTGCCCTCCTGGc                                                                                                                                      <  2:231039‑M1/18‑1 (MQ=255)
              gCTGAACTTCGCCGAGCATGCCGGAAACGAAGTCGCATTTTCCTACAACTGACTCAGAAACCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGc                                                                                                                                      >  2:281210/1‑150 (MQ=255)
                tGAACTTCGCCGAGCATGCCGGAAACGAAGTCGCATTTTCCTACAACTGACTCAGAAACCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCAc                                                                                                                                    >  2:1392179/1‑150 (MQ=255)
                    tagaaggctgagcgcacagcgtcgatccggcgcaggttgaggcccagcgaggcgacttcttcattgaggctggtgcggaactggccaccgcttgggtagctcaccatggcgatgaaacgtgactcaATCTCGTGCCCTCCTGGCACgaga                                                                                                                                >  2:60652‑M1/127‑150 (MQ=255)
                     agaaggctgagcgcacagcgtcgatccggcgcaggttgaggcccagcgaggcgacttcttcattgaggctggtgcggaactggccaccgcttgggtagcccaccatggcgatgacacgtgactcgATCTCGTGCCCTCCTGGCACGAGAt                                                                                                                               >  1:231039‑M1/126‑150 (MQ=255)
                     agaaggctgagcgcacagcgtcgatccggcgcaggttgaggcccagcgaggcgacttcttcattgaggctggtgcggaactggccaccgcttgggtagcccaccatggcgatgacacgtgactcgACCTCGTGCCCTCCTGGCACGAGAt                                                                                                                               >  1:65338‑M1/126‑150 (MQ=255)
                         ccGAGCATGCCGGAAACGAAGTCGCATTTTCCTACAACTGACTCAGAAACCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAg                                                                                                                           <  1:761779/150‑1 (MQ=255)
                          cGAGCATGCCGGAAACGAAGTCGCATTTTCCTACAACTGACTCAGAAACCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGt                                                                                                                          <  1:1392179/150‑1 (MQ=255)
                                            aGTCGCATTTTCCTACAACTGACTCAGAAACCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGAACTACCTCGTTCCGTCCTGGCACGAGATCGAGTCACGTGTCATCGCCAtgg                                                                                                        <  1:32112/150‑1 (MQ=255)
                                                gCATTTTCCTACAACTGACTCAGAAACCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTggg                                                                                                    >  1:1153844/1‑150 (MQ=39)
                                                gCATTTTCCTACAACTGACTCAGAAACCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTAACCTGTCCACGCCATGTTggg                                                                                                    >  2:41522/1‑150 (MQ=11)
                                                 gcgcaggttgaggcccagcgaggcgacttcttcattgaggctggtgcggaactggccaccgcttgggtagcccaccatggcgatgacacgtgactcgATCTCGTGCCCTCCTGGCACgag                                                                                                                                 <  2:592450‑M1/23‑1 (MQ=255)
                                                 gcgcaggttgaggcccagcgaggcgacttcttcattgaggctggtgcggaactggccaccgcttgggtagcccaccatggcgatgacacgtgactcgATCTCGTGCCCTCCTGGCACgag                                                                                                                                 >  1:592450‑M1/98‑120 (MQ=255)
                                                 cATTTTCCTACAACTGACTCAGAAACCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGc                                                                                                   >  2:96655/1‑150 (MQ=39)
                                                     ttCCTACAACTGACTCAGAAACCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTAcc                                                                                               <  1:1246578/150‑1 (MQ=38)
                                                        cTACAACTGACTCAGAAACCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGcc                                                                                                            <  1:823045/134‑1 (MQ=40)
                                                        cTACAACTGACTCAGAAACCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGcc                                                                                                            >  2:823045/1‑134 (MQ=40)
                                                        cTACAACTGACTCAGAAACCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCaa                                                                                            >  1:1317214/1‑150 (MQ=255)
                                                          aCAACTGACTCAGAAACCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGc                                                                                          <  2:475823/150‑1 (MQ=40)
                                                            aaCTGACTCAGAAACCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCgg                                                                                        <  1:500496/150‑1 (MQ=255)
                                                                   tCAGAAACCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAg                                                                                 >  1:853419/1‑150 (MQ=35)
                                                                   tCAGAAACCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAg                                                                                 <  1:1183628/150‑1 (MQ=35)
                                                                     aGAAACCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGtt                                                                               >  2:404050/1‑150 (MQ=35)
                                                                       aaaCCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTcc                                                                             <  1:404050/150‑1 (MQ=35)
                                                                       aaaCCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTcc                                                                             >  1:158704/1‑150 (MQ=35)
                                                                           cTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCa                                                                                                                                     <  1:428311/90‑1 (MQ=35)
                                                                           cTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCa                                                                                                                                     >  2:428311/1‑90 (MQ=33)
                                                                            tccgatccgaATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCg                                                                            >  2:893380/1‑146 (MQ=35)
                                                                            tccgatccgaATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCg                                                                            <  1:893380/146‑1 (MQ=34)
                                                                            tccgatccgaATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCAcc                                                                        >  2:1245362/1‑150 (MQ=35)
                                                                            tccgatccgaATCCATGATGGGATTTCCCCCTGAACACGTTGTGAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCAcc                                                                        >  2:747235/1‑150 (MQ=25)
                                                                             ccgatccgaATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCa                                                                       <  2:1153844/150‑1 (MQ=35)
                                                                                 tccgaATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCTCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCt                                                                   <  1:378381/150‑1 (MQ=25)
                                                                                       tCCATGATGGGTTGTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATga                                                             >  2:1096331/1‑150 (MQ=25)
                                                                                           tgatgGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATgaagaa                                                         <  1:96655/150‑1 (MQ=37)
                                                                                            gatgGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGaagaag                                                        <  2:610925/150‑1 (MQ=37)
                                                                                             atgGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGt                                                       >  1:1413873/1‑150 (MQ=38)
                                                                                              tgGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGCGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTCCCGCAGCAGCCTCAATGAAGAAGtg                                                      >  1:346328/1‑149 (MQ=255)
                                                                                                ggATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGc                                                    >  1:547038/1‑150 (MQ=39)
                                                                                                 gATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCCACACCTCCTGGCACGAGATGGATTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGcc                                                   >  1:610182/1‑150 (MQ=11)
                                                                                                  aTTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCt                                                  <  1:281210/150‑1 (MQ=39)
                                                                                                   tttCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTc                                                 >  2:1132770/1‑150 (MQ=39)
                                                                                                       gccaccgctttggtagcccaccatggcgatgacacgtgactcgACCTCGTGCCCTCCTGGCACgaga                                                                                                                                >  2:654438‑M1/44‑67 (MQ=255)
                                                                                                       gccaccgctttggtagcccaccatggcgatgacacgtgactcgACCTCGTGCCCTCCTGGCACgaga                                                                                                                                <  1:654438‑M1/24‑1 (MQ=255)
                                                                                                        cccTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTgg                                            <  2:1317214/150‑1 (MQ=255)
                                                                                                          cTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCGATCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGc                                          <  1:41522/150‑1 (MQ=255)
                                                                                                             gctttggtagcccaccatggcgatgacacgtgactcgACCTCGTGCCCTCCTGGCACGAGAt                                                                                                                               <  1:999761‑M1/25‑1 (MQ=255)
                                                                                                             gctttggtagcccaccatggcgatgacacgtgactcgACCTCGTGCCCTCCTGGCACGAGAt                                                                                                                               >  2:999761‑M1/38‑62 (MQ=255)
                                                                                                                acGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAAc                                    >  1:1199819/1‑150 (MQ=255)
                                                                                                                  gTTGACAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCGAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCt                                  >  1:1439895/1‑150 (MQ=255)
                                                                                                                   gtagcccaccatggcgatgacacgtgactcgATCTCGTGCcctcc                                                                                                                                          >  2:305373‑M1/32‑45 (MQ=255)
                                                                                                                   gtagcccaccatggcgatgacacgtgactcgATCTCGTGCcctcc                                                                                                                                          <  1:305373‑M1/14‑1 (MQ=255)
                                                                                                                      tCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTgcgc                              <  2:547038/150‑1 (MQ=255)
                                                                                                                         gggAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAAt                                                               >  1:1474351/1‑114 (MQ=38)
                                                                                                                         gggAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAAt                                                               <  2:1474351/114‑1 (MQ=38)
                                                                                                                            agaCACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATc                        <  1:1245362/150‑1 (MQ=255)
                                                                                                                                acCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATgaa                                                            >  2:953148/1‑110 (MQ=38)
                                                                                                                                acCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATgaa                                                            <  1:953148/110‑1 (MQ=38)
                                                                                                                                   tCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCtg                 >  1:566324/1‑150 (MQ=255)
                                                                                                                                       gTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCACCCTGCGCCGGATCGACGCTGTgcg             >  1:508732/1‑150 (MQ=255)
                                                                                                                                                ttACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCtt    >  1:883448/1‑150 (MQ=255)
                                                                                                                                                ttACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCt     <  1:72138/149‑1 (MQ=255)
                                                                                                                                                ttACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCt     >  2:72138/1‑149 (MQ=255)
                                                                                                                                                  ACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCt  >  1:747235‑M1/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCAACAACTGAACAGCTGAACTTCGCCGAGCATGCCGGAAACGAAGTCGCATTTTCCTACAACTGACTCAGAAACCTCCGATCCGAATCCATGATGGGATTTCCCCCTGCACACGTTGTCAGGGAGACACCTCATGTACCTGGATTACCTCGTGCCCTCCTGGCACGAGATCGAGTCACGTGTCATCGCCATGGTGGGCTACCCAAGCGGTGGCCAGTTCCGCACCAGCCTCAATGAAGAAGTCGCCTCGCTGGGCCTCAACCTGCGCCGGATCGACGCTGTGCGCTCAGCCTTCT  >  NC_002947/898005‑898300

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