Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 27,581 C→T 6.2% intergenic (‑428/+18) PP_0022 ← / ← PP_0023 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294727,5810CT 6.2% 187.3 / 2.1 66intergenic (‑428/+18)PP_0022/PP_0023hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (27/33);  new base T (2/2);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATGCTCATG  >  NC_002947/27439‑27728
                                                                                                                                               |                                                                                                                                                   
ctcGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACCGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAgc                                                                                                                                              <  1:268260/150‑1 (MQ=255)
    tcgtcCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAgctgct                                                                                                                                          >  1:40362/1‑150 (MQ=255)
      gtcCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCttt                                                                                                                                        >  1:9359/1‑150 (MQ=255)
       cccGCCCTGAAGCGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCtttt                                                                                                                                       <  2:1055635/149‑1 (MQ=255)
         cGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCtttttt                                                                                                                                     <  1:586937/150‑1 (MQ=255)
           cACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCt                                                                                                                                   <  1:37402/150‑1 (MQ=255)
                    gATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTccccc                                                                                                                          <  1:521658/150‑1 (MQ=255)
                          cGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGagag                                                                                                                                               >  1:1204009/1‑123 (MQ=255)
                          cGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGagag                                                                                                                                               <  2:1204009/123‑1 (MQ=255)
                           gATCCGGAGGGAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTtgt                                                                                                                   <  1:127182/150‑1 (MQ=255)
                                ggAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTtg                                                                                                              >  2:733840/1‑150 (MQ=255)
                                        cTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCgag                                                                                                                                                 <  1:274635/107‑1 (MQ=255)
                                        cTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCgag                                                                                                                                                 >  2:274635/1‑107 (MQ=255)
                                        cTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGttt                                                                                                      <  1:438088/150‑1 (MQ=255)
                                          tCAGAACCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTNCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCa                                                                                                    <  1:657909/150‑1 (MQ=255)
                                           cAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAg                                                                                                   >  1:285184/1‑150 (MQ=255)
                                             gAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGt                                                                                                 >  1:748968/1‑150 (MQ=255)
                                              aaTCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTc                                                                                                >  1:106788/1‑150 (MQ=255)
                                               aTCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTcc                                                                                               <  2:135288/150‑1 (MQ=255)
                                               aTCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTcc                                                                                               >  1:135288/1‑150 (MQ=255)
                                                         ccaagcGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCTCGGGGGCTGCtttt                                                                                                                                       >  1:660705/1‑100 (MQ=255)
                                                         ccaagcGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCTCGGGCGCTGCtttt                                                                                                                                       <  2:660705/100‑1 (MQ=255)
                                                         ccaagcGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCt                                                                                     >  2:1140691/1‑150 (MQ=255)
                                                            agcGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCt                                                                                  <  2:106788/150‑1 (MQ=255)
                                                                 aaGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTccccc                                                                                                                          >  2:221173/1‑105 (MQ=255)
                                                                 aaGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTccccc                                                                                                                          <  1:221173/105‑1 (MQ=255)
                                                                     cTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAg                                                                         <  2:748968/150‑1 (MQ=255)
                                                                        acgacgCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAgggg                                                                      >  1:699062/1‑150 (MQ=255)
                                                                        acgacgCGGGCCGATTGATCAGCGAGCCGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCCGCTCGCGAGATGTGCAGTTTTCTTACTCCGCCTTTTTGTGCTTGGGGGGTTTCAGGTCCGGATGGTCCGTTGGCTTTTCAgggg                                                                      >  2:562954/1‑150 (MQ=255)
                                                                             gcgGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACt                                                                 >  2:54970/1‑150 (MQ=255)
                                                                                 tCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCt                                                             >  1:1065816/1‑150 (MQ=255)
                                                                                   cGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTAAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGc                                                           >  2:37873/1‑150 (MQ=255)
                                                                                        gATCAGCGAGACGGTCGGGAGGGTTTGCAAACAGACCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTCCTTACTCCCCCTTTTTGTTCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCtttt                                                      <  2:40362/150‑1 (MQ=255)
                                                                                         aTCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTc                                                                           >  1:6331/1‑128 (MQ=255)
                                                                                         aTCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTc                                                                           <  2:6331/128‑1 (MQ=255)
                                                                                         aTCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTg                                                     >  2:1015543/1‑150 (MQ=255)
                                                                                            aGCGAGACTGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTcc                                                  <  1:54970/150‑1 (MQ=255)
                                                                                            aGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCACTTGAAGGGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTcc                                                  >  1:1044162/1‑150 (MQ=255)
                                                                                             gCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCg                                                 >  1:649585/1‑150 (MQ=255)
                                                                                                 gaCGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGt                                             >  2:981811/1‑150 (MQ=255)
                                                                                                 gaCGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTAAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGt                                             <  1:37873/150‑1 (MQ=255)
                                                                                                    ggTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGcc                                          >  1:1114742/1‑150 (MQ=255)
                                                                                                         tGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCTGAGTGCCATCtt                                     >  1:264243/1‑150 (MQ=255)
                                                                                                          gAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTc                                    <  2:285184/150‑1 (MQ=255)
                                                                                                           aGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTc                                                                                                >  2:423099/1‑89 (MQ=255)
                                                                                                           aGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTc                                                                                                <  1:423099/89‑1 (MQ=255)
                                                                                                            gggTTTGCACACAGCCAC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCAt                                  <  2:9359/150‑1 (MQ=255)
                                                                                                              gTTTGCACACAGCCCCAAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCAtt                                 >  2:653423/1‑150 (MQ=255)
                                                                                                              gTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTg                                <  2:1065816/150‑1 (MQ=255)
                                                                                                                         gCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTgg                     <  2:264243/150‑1 (MQ=255)
                                                                                                                         gCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTgg                     <  1:555543/150‑1 (MQ=255)
                                                                                                                           ccc‑AAAAAAACAGCAGCTCCCGCGAGCTGCTTTTTGCTTACTCCCCCTTTTTGTGCTTGTGGGGTTACAGGTCCGGATGGTCCTTCGGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTCCCCTCTTCATTGGTGCGCTTTGGGc                   <  1:562954/150‑1 (MQ=255)
                                                                                                                           ccc‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTTTTTTTTTTGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTTAGGGGTGACTCTCTGCCTTTTGTTCGGAGTGCCATCTTCATTGGTGCGCTTTTGGc                   >  2:122831/1‑150 (MQ=255)
                                                                                                                            cc‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAgggg                                                                      <  1:759550/98‑1 (MQ=255)
                                                                                                                            cc‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAgggg                                                                      >  2:759550/1‑98 (MQ=255)
                                                                                                                             c‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCt                 <  1:1140691/150‑1 (MQ=255)
                                                                                                                             c‑AAAAAAAAAGAAGCTCACGAGAGCTGCTGGTTTGTTACTCGACCTTGTTGTGCTTGTGGGGTTTCAGGTCGGGATGGTCTTTCTGCTTTTAAGGGGTGACGAGGTACCTTTTGTCCGGAGTGCTATCTTCATTGTTGCGGTTTTGGCCt                 >  2:413951/1‑150 (MQ=255)
                                                                                                                                 aaaaaaaGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGc              <  2:514998/150‑1 (MQ=255)
                                                                                                                                    aaaaGCAGCTCTCGGGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAggg                                                                       >  1:69529/1‑90 (MQ=255)
                                                                                                                                    aaaaGCAGCTCTCGGGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAggg                                                                       <  2:69529/90‑1 (MQ=255)
                                                                                                                                    aaaaGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTg           >  1:749761/1‑150 (MQ=255)
                                                                                                                                       agcagcTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATg        <  2:699062/150‑1 (MQ=255)
                                                                                                                                       agcagcTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTGTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATg        <  1:122831/150‑1 (MQ=255)
                                                                                                                                          agcTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATGCTc     <  1:733840/150‑1 (MQ=255)
                                                                                                                                            cTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATGCTCAt   <  1:981811/150‑1 (MQ=255)
                                                                                                                                             tCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATGCTCATg  <  2:1044162/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CTCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCC‑AAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATGCTCATG  >  NC_002947/27439‑27728

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