Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,958,228 A→T 9.1% L734H (CTT→CAT) ‡ mnxG ← manganese‑oxidizing multicopper oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,958,2280AT 9.1% 232.1 / 12.3 77L734H (CTT→CAT) ‡mnxGmanganese‑oxidizing multicopper oxidase
Reads supporting (aligned to +/- strand):  ref base A (37/33);  new base T (3/4);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

TGCATGGCTGCCTGCTCGACATCGGTGCCCTCTTCCGGCAGGTAGATGGGCTTGGCCTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCGGGCCAGGGCCGGGTCGAGCATGTCCAGGGGCGGGGTTGGCGGGCGGTTGCCCACGCTGCTTTCCATGCCACC  >  NC_002947/3958084‑3958376
                                                                                                                                                |                                                                                                                                                    
tGCATGGCTGCCTGCTCGACATCGGTGCCCTCTTCCGGCAGGTAGATGGGCTTGGCCTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGgccg                                                                                                                                                 <  2:963307/150‑1 (MQ=255)
  cATGGCTGCCTGCTCGACATCGGTGCCCTCTTCCGGCAGGTAGATGGGCTTGGCCTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGgccgcc                                                                                                                                               >  1:1257652/1‑150 (MQ=255)
   aTGGCTGCCTGCTCGACATCGGTGCCCTCTTCCGGCAGGTAGATGGGCTTGGCCTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGccgccg                                                                                                                                              <  2:441613/150‑1 (MQ=255)
   aTGGCTGCCTGCTCGACATCGGTGCCCTCTTCCGGCAGGTAGATGGGCTTGGCCTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCACGGCGGCCTTGgccgcca                                                                                                                                              <  2:292507/150‑2 (MQ=255)
    tGGCTGCCTGCTCGACATCGGTGCCCTCTTCCGGCAGGTAGATGGGCTTGGCCTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGt                                                                                                                                             >  1:998369/1‑150 (MQ=255)
      gCTGCCTGCTCGACATCGGTGCCCTCTTCCGGCAGGTAGATGGGCTTGGCCTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTcc                                                                                                                                           <  1:917727/150‑1 (MQ=255)
      gCTGCCTGCTCGACATCGGTGCCCTCTTCCGGCAGGTAGATCGGCTTGGCCTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTcc                                                                                                                                           <  2:975561/150‑1 (MQ=255)
              cTCGACATCGGTGCCCTCTTCCGGCAGGTAGATGGGCTTGGCCTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATc                                                                                                                                   <  1:723346/150‑1 (MQ=255)
               tCGACATCGGTGCCCTCTTCCGGCAGGTAGATGGGCTTGGCCTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCg                                                                                                                                  >  1:152078/1‑150 (MQ=255)
                  aCATCGGTGCCCTCTTCCGGCAGGTAGATGGGCTTGGCCTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGAcc                                                                                                                               >  1:1252788/1‑150 (MQ=255)
                        gTGCCCTCTTCCGGCAGGTAGATGGGCTTGGCCTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTc                                                                                                                         >  2:264918/1‑150 (MQ=255)
                        gTGCCCTCTTCCGGCAGGTAGATGGGCTTGGCCTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTc                                                                                                                         <  2:253689/150‑1 (MQ=255)
                            ccTCTTCCGGCAGGTAGATGGGCTTGGTCTTGTGCACGACCTTGGTGAAGTCCTGTTTGGTGGTGGTGTTAGCCGCTACGCGACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCAt                                                                                                                     >  1:677750/1‑150 (MQ=255)
                             ctctTCCGGCAGGTAGATGGGCTTGGCCTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCTTGgccgccat                                                                                                                                             >  2:685729/1‑123 (MQ=255)
                             ctctTCCGGCAGGTAGATGGGCTTGGCCTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCTTGgccgccat                                                                                                                                             <  1:685729/125‑3 (MQ=255)
                             ctctTCCGGCAGGTAGATGGGCTTGGCCTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATc                                                                                                                    >  2:822018/1‑150 (MQ=255)
                              tctTCCGGCAGGTAGATGGGCTTGGCCTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCTTGgccgccat                                                                                                                                             >  2:1376177/1‑122 (MQ=255)
                              tctTCCGGCAGGTAGATGGGCTTGGCCTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCTTGgccgccat                                                                                                                                             <  1:1376177/124‑3 (MQ=255)
                               ctTCCGGCAGGTAGATGGGCTTGGCCTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCa                                                                                                                  <  2:1257652/150‑1 (MQ=255)
                                ttCCGGCAGGTAGATGGGCTTGGCCTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAg                                                                                                                                                     <  2:674556/114‑1 (MQ=255)
                                ttCCGGCAGGTAGATGGGCTTGGCCTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGACGTCCAGTGCATGGCGGCCAAg                                                                                                                                                     >  1:674556/1‑114 (MQ=255)
                                 tCCGGCAGGTAGATGGGCTTGGCCTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAgg                                                                                                                <  1:1081100/150‑1 (MQ=255)
                                  ccGGCAGGTAGATGGGCTTGGCCTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGt                                                                                                               <  2:788411/150‑1 (MQ=255)
                                     gCAGGTAGATGGGCTTGGCCTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTgg                                                                                                            >  1:1162466/1‑150 (MQ=255)
                                           aGATGGGCTTGGCCTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCcaca                                                                                                      <  2:1252788/150‑1 (MQ=255)
                                            gATGGGCTTGGCCTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAg                                                                                                     >  1:333178/1‑150 (MQ=255)
                                            gATGGGCTTGGCCTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAg                                                                                                     >  2:324358/1‑150 (MQ=255)
                                              tGGGCTTGGCCTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAggg                                                                                                   >  1:527978/1‑150 (MQ=255)
                                                   ttaaccctgtgCAAGAGCTTGGTGGAGTCGAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCGGCCATCGACCTGGTTCGCATCCAGGTTGGCCCAGAGGGCTTTa                                                                                              >  2:448861/8‑150 (MQ=255)
                                                    tGGCCTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTAc                                                                                             <  2:1162466/150‑1 (MQ=255)
                                                        cTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAggg                                                                                                   >  1:592969/1‑140 (MQ=255)
                                                        cTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAggg                                                                                                   <  2:592969/140‑1 (MQ=255)
                                                        cTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACt                                                                                         <  2:998369/150‑1 (MQ=255)
                                                          tgtgCACGACCTTGGTGAAGGCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGt                                                                                       <  2:152078/150‑1 (MQ=255)
                                                              cACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCcaca                                                                                                      <  2:633970/131‑1 (MQ=255)
                                                              cACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCcaca                                                                                                      >  1:633970/1‑131 (MQ=255)
                                                              cACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAggg                                                                                                   >  1:1135135/1‑134 (MQ=255)
                                                              cACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAggg                                                                                                   <  2:1135135/134‑1 (MQ=255)
                                                                 gACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAg                                                                                >  1:1300168/1‑150 (MQ=255)
                                                                    cTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATTGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTAc                                                                                             >  2:888135/1‑134 (MQ=255)
                                                                    cTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATTGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTAc                                                                                             <  1:888135/134‑1 (MQ=255)
                                                                         tGAAGTCCCGATTGGAGCTGGTGGTATCCGCTACGCCACGGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCggg                                                                        >  2:703388/1‑150 (MQ=255)
                                                                         tGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCggg                                                                        >  1:1343535/1‑150 (MQ=255)
                                                                          gAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCGGGc                                                                       >  2:889843/1‑150 (MQ=255)
                                                                              tgctgTTTGGTGGTGGTGTTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCGGGCCAgg                                                                   >  2:733220/5‑150 (MQ=255)
                                                                                cAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCTTGgccgcca                                                                                                                                              >  1:851990/1‑72 (MQ=255)
                                                                                cAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCTTGgccgcca                                                                                                                                              <  2:851990/73‑2 (MQ=255)
                                                                                cAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCGGGCCAGGGc                                                                 >  2:890738/1‑150 (MQ=255)
                                                                                 aGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCa                                                                                                                  >  2:1109339/1‑100 (MQ=255)
                                                                                 aGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCa                                                                                                                  <  1:1109339/100‑1 (MQ=255)
                                                                                  gTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCGGGCCAGGGCCg                                                               >  1:294957/1‑150 (MQ=255)
                                                                                   tttGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCGGGCCAGGGccag                                                              <  1:837550/150‑3 (MQ=255)
                                                                                        tggtggtggtATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCGGGCCAGGGCCGGGTCGa                                                         <  1:1049124/150‑1 (MQ=255)
                                                                                        tggtggtggtATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCGGGCCAGGGCCGGGTCGa                                                         <  1:890738/150‑1 (MQ=255)
                                                                                              tggtATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCGGGCCAGGGCCGGGTCGAGCATGt                                                   >  1:59568/1‑150 (MQ=255)
                                                                                                gtATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCGGGCCAGGGCCGGGTCGAGCATGTcc                                                 >  1:1115928/1‑150 (MQ=255)
                                                                                                   tCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACGGTCGTTCAGTTGTCGGGCCAGGGCCGGGTCGAGCATGTCCAgg                                              >  2:1328287/1‑150 (MQ=255)
                                                                                                    ccGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCGGGCCAGGGCCGGGTCGAGCATGTCCAggt                                             <  1:889843/150‑2 (MQ=255)
                                                                                                         aCGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCGGGCCAGGGCCGGGTCGAGCATGTCCAGGGGCggg                                        <  2:294957/150‑1 (MQ=255)
                                                                                                         aCGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCAATCGACTTGCTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCGGGCCAGGGCCGGGTCGAGCATGTCCAGGGGCggg                                        <  1:324358/150‑1 (MQ=255)
                                                                                                           gCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCGGGCCAGGGCCGGGTCGAGCATGTCCAGGGGCGGGGt                                      >  2:761124/1‑150 (MQ=255)
                                                                                                             cACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCGGGCCAGGGCCGGGTCGAGc                                                       >  1:832298/1‑131 (MQ=255)
                                                                                                             cACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCGGGCCAGGGCCGGGTCGAGc                                                       <  2:832298/131‑1 (MQ=255)
                                                                                                             cACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCGGGCCAGGGCCGGGTCGAGCATGTCCATGGGCGGGGTTg                                    >  2:961937/1‑150 (MQ=255)
                                                                                                                 ggcggAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCGGGCCAGGGCCGGGTCGAGCATGTCCAGGTGCGGGGTTggcgg                                <  2:677750/150‑1 (MQ=255)
                                                                                                                    ggAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCGGGCCAGGGCCGGGTCGAGCATGTCCAGGGGCGGGGTTggcggcag                             <  1:761124/150‑4 (MQ=255)
                                                                                                                       gaCGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCGGGCCAGGGCCGGGTCGAGCATGTCCAGGGGCGGGGTTGGCGGGCGGtt                          <  1:733220/150‑1 (MQ=255)
                                                                                                                       gaCGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCGGGCCAGGGCCGGGTCGAGCATGTCCAGGGGCGGGGTTGGCGGGCGGtt                          <  2:59568/150‑1 (MQ=255)
                                                                                                                          gCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCGGGCCAGGGCCGGGTCGAGCATGt                                                   >  1:1419124/1‑122 (MQ=255)
                                                                                                                          gCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCGGGCCAGGGCCGGGTCGAGCATGt                                                   <  2:1419124/122‑1 (MQ=255)
                                                                                                                          gCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCGGGCCAGGGCCGGGTCGAGCATGTCCAGGGGCGGGGTTGGCGGGCGGTTGcc                       >  2:997063/1‑150 (MQ=255)
                                                                                                                             gTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCTGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCGGGCCAGGGCCGGGTCGAGCATGTCCAGGGGCGGGGTTGGCGGGCGGTTGCCCAc                    >  2:331280/1‑150 (MQ=255)
                                                                                                                                 aGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTGGACCACTGTCTTTCAGTTGCCGGGCCAGGGCCGGGTCGAGCATGTCCAGGGGCGGGGTTGGCGGGCGGTTGCCCACgctg                >  1:460899/1‑150 (MQ=255)
                                                                                                                                     cATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCGGGCCAGGGCCGGGTCGAGCATGTCCAGGGGCGGGGTTGGCGGGCGGTTGCCCACGCTGCttt            >  1:1179286/1‑150 (MQ=255)
                                                                                                                                       tGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCGGGCCAGGGCCGGGTCGAGCATGTCCAGGGGCGGGGTTGGCGGGCGGTTGCCCACGCTGCTTTcc          <  2:1115928/150‑1 (MQ=255)
                                                                                                                                       tGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCGGGCCAGGGCCGGGTCGAGCATGTCCAGGGGCGGGGTTGGCGGGCGGTTGCCCACGCTGCTTTcc          >  1:997040/1‑150 (MQ=255)
                                                                                                                                           ggcCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCTGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCGGGCCAGGGCCGGGTCGAGCATGTCCAGGGGCGGGGTTGGCGGGCGGTTGCCCACGCTGCTTTCCATGc      <  1:331280/150‑1 (MQ=255)
                                                                                                                                               aaGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCGGGCCAGGGCCGGGTCGAGCATGTCCATGGGCGGGGTTGGCGGGCGGTTGCCCACGCTGCTTTCCATGCCAcc  <  1:961937/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TGCATGGCTGCCTGCTCGACATCGGTGCCCTCTTCCGGCAGGTAGATGGGCTTGGCCTTGTGCACGACCTTGGTGAAGTCCAGTTTGGTGGTGGTGGTATCCGCTACGCCACCGGCGGAGACGCCGTCCAGTGCATGGCGGCCAAGGCCGCCGTCCCAGCCATCGACCTGGTTCGCATCCAGGTTGGCCCACAGGGCTTTACCACTGTCTTTCAGTTGCCGGGCCAGGGCCGGGTCGAGCATGTCCAGGGGCGGGGTTGGCGGGCGGTTGCCCACGCTGCTTTCCATGCCACC  >  NC_002947/3958084‑3958376

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