Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 27,581 C→T 11.7% 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,5810CT11.7% 148.6 / 11.2 60intergenic (‑428/+18)PP_0022/PP_0023hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (26/26);  new base T (3/4);  total (29/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.93e-01

TCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATGCTCATGCA  >  NC_002947/27440‑27730
                                                                                                                                              |                                                                                                                                                     
tcgtcgTCCGCACTGAAGGGTTAACGGATCGGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAgct                                                                                                                                               >  2:922571/1‑150 (MQ=255)
    ggtcCGCACTGAATGGATAACCGATCCGGAGGTAACTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCAAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCtt                                                                                                                                           <  2:824556/149‑1 (MQ=255)
     gtcCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCACCAAAAAAGCAGCTCTCGGGAGCTGCttt                                                                                                                                          <  1:385826/150‑1 (MQ=255)
        cGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCtttttt                                                                                                                                       <  2:923711/150‑1 (MQ=255)
           aCTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCtt                                                                                                                                    >  2:1355736/1‑150 (MQ=255)
            cTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGAACGTCTGGAAGGGTTGCACACAGGCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTa                                                                                                                                   <  2:56446/150‑1 (MQ=255)
              gAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACt                                                                                                                                 >  1:1044072/1‑150 (MQ=255)
               aaGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTc                                                                                                                                >  1:1022819/1‑150 (MQ=255)
                aGGTATAACCCCTACGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTc                                                                                                                                <  1:743214/150‑1 (MQ=255)
                               ggAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAATACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGTGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTtg                                                                                                                <  2:520683/150‑1 (MQ=255)
                                 atgTAGCGTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTtgtg                                                                                                              <  1:1258369/148‑1 (MQ=255)
                                  ggTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGAGAGACGGTCGTGAGGGTTTGCACACAGCCACAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTgg                                                                                                             <  2:219595/150‑1 (MQ=255)
                                       cTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCt                                                                                                                                     <  2:889782/121‑1 (MQ=255)
                                       cTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCt                                                                                                                                     >  1:889782/1‑121 (MQ=255)
                                        ttCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTc                                                                                                       <  2:1044072/150‑1 (MQ=255)
                                          cAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCa                                                                                                      >  1:761695/1‑150 (MQ=255)
                                          cAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAg                                                                                                     <  2:959632/150‑1 (MQ=255)
                                            gAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGt                                                                                                   >  2:1537612/1‑150 (MQ=255)
                                             aaTCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTc                                                                                                  <  2:1022819/150‑1 (MQ=255)
                                              aTCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCTCCAAAAAAAAAGCAGCTCCCTAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTcc                                                                                                 <  2:1063742/150‑1 (MQ=255)
                                                 cATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGa                                                                                              <  1:922571/150‑1 (MQ=255)
                                                  aTCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGCGATGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGAt                                                                                             >  1:1186849/1‑150 (MQ=255)
                                                   tCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATg                                                                                            >  1:86531/1‑150 (MQ=255)
                                                     aagccaagcGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGt                                                                                          >  2:740141/1‑150 (MQ=255)
                                                      agccaagcGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTc                                                                                         <  1:740141/150‑1 (MQ=255)
                                                        ccaagcGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTa                                                                                                                                   <  1:655213/106‑1 (MQ=255)
                                                        ccaagcGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTa                                                                                                                                   >  2:655213/1‑106 (MQ=255)
                                                        ccaagcGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCt                                                                                       >  2:1304569/1‑150 (MQ=255)
                                                              ggAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAGAAAAAAGCAGCTCTCGGGAGCTGCttttttttt                                                                                                                                    >  1:73187/1‑96 (MQ=37)
                                                              ggAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAGAAAAAAGCAGCTCTCGGGAGCTGCTtttttttt                                                                                                                                    <  2:73187/99‑1 (MQ=37)
                                                                       acgacgCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTc                                                                                                                                      <  1:62463/88‑1 (MQ=255)
                                                                       acgacgCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTc                                                                                                                                      >  2:62463/1‑88 (MQ=255)
                                                                        cgacgaGCTCCGATTGATCAGGGAGACGGTAGTGAGGGTTTGCACACAGCCCCAAAGAAAAAGCAGCTCACGAGAGCTGCTTTTTTCTTACTTCGCCTTTTTGTGCTTGTGGGGTTTCAGGTTCGTATGGTCCTTCTGCTTTTCAGGGGt                                                                       <  2:683779/150‑1 (MQ=255)
                                                                         gacgCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCATGTCAGGGGTg                                                                      >  1:1485998/1‑150 (MQ=255)
                                                                          acgCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGt                                                                       >  2:328413/1‑148 (MQ=255)
                                                                          acgCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGt                                                                       <  1:328413/148‑1 (MQ=255)
                                                                                       taTCAGCCACACGTTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGATAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCtttt                                                        <  1:590918/149‑1 (MQ=255)
                                                                                         tCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTTCGGATGGTCCTTTTGCTTTTTAGGGGTGAATCGCTGCCtttttt                                                      >  2:770671/1‑148 (MQ=255)
                                                                                            gCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTc                                                                                         >  2:454796/1‑113 (MQ=255)
                                                                                            gCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTc                                                                                         <  1:454796/113‑1 (MQ=255)
                                                                                             cGAGACGGTCGTGAGGGTTTGCACCCAGCCCCAAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGATGCCTTTTGTCCg                                                   <  1:530387/150‑1 (MQ=255)
                                                                                                  cGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTTCTTGTTGGATTTCAAGTCTGGCTGGTCTTTATACTTTTCAGGGGGGACTTGCTGCCTTTTGTCCGGAggtc                                             >  2:658046/1‑147 (MQ=255)
                                                                                                    gTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTGTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCa                                           >  2:718245/1‑150 (MQ=255)
                                                                                                         gAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTc                                      >  1:797453/1‑150 (MQ=255)
                                                                                                          agtgTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCa                                     >  2:1642585/4‑150 (MQ=255)
                                                                                                               ttGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGt                                <  2:86531/150‑1 (MQ=255)
                                                                                                                   acacaGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTgcgc                            >  2:13578/1‑150 (MQ=255)
                                                                                                                     acaGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCtt                          >  1:964089/1‑150 (MQ=255)
                                                                                                                          cccAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGc                     <  1:1537612/150‑1 (MQ=255)
                                                                                                                          cccAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGACGGTcc                                                                                        >  1:788859/1‑83 (MQ=255)
                                                                                                                          cccAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGACGGTcc                                                                                        <  2:788859/83‑1 (MQ=255)
                                                                                                                                gaaaaaaGCAGCTCTCGGGAGCTGCTTTTTTCTTACTCTCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTgg                 >  2:1522422/2‑149 (MQ=255)
                                                                                                                                gaaaaaaGCAGCTCTCGGGAGCTGCTTTTTTCTTACTCTCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTgg                 <  1:1522422/148‑1 (MQ=255)
                                                                                                                                   aaaaGCAGCTCTCGGGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTg             >  2:728998/1‑150 (MQ=255)
                                                                                                                                         agcTCTCGGGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATGCTc       <  1:728998/150‑1 (MQ=255)
                                                                                                                                         agcTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTc                                      <  2:240246/119‑1 (MQ=255)
                                                                                                                                         agcTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTc                                      >  1:240246/1‑119 (MQ=255)
                                                                                                                                         agcTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATGCTc       <  1:718245/150‑1 (MQ=255)
                                                                                                                                           cTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGa                                                                                              >  2:1306921/1‑61 (MQ=255)
                                                                                                                                           cTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGa                                                                                              <  1:1306921/61‑1 (MQ=255)
                                                                                                                                              ccGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATGCTCATGCa  >  1:384431/1‑150 (MQ=255)
                                                                                                                                              ccGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATGCTCATGCa  >  1:1073643/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
TCGTCGTCCGCACTGAAGGGATAACCGATCCGGAGGTAGCTTCAGAATCCATCAAGCCAAGCGGAAGACTCACGACGCGGTCCGATTGATCAGCGAGACGGTCGTGAGGGTTTGCACACAGCCCCAAAAAAAAAGCAGCTCCCGAGAGCTGCTTTTTTCTTACTCCCCCTTTTTGTGCTTGTGGGGTTTCAGGTCCGGATGGTCCTTCTGCTTTTCAGGGGTGACTCGCTGCCTTTTGTCCGGAGTGCCATCTTCATTGGTGCGCTTTGGGCCTGGCTTGATGCTCATGCA  >  NC_002947/27440‑27730

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