Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,043,691 A→C 6.5% D27A (GAC→GCC)  PP_0903 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,043,6910AC 6.5% 162.3 / 2.3 61D27A (GAC→GCC) PP_0903hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/30);  new base C (2/2);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

ACATTGGCATGACCTGATTCCAACGCTTGCATGACTAGTGTGGGCGCAATGGAGTCAGGAGATGCCCAATGCGTAGCCCCATCAAGCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCATTGGCGGGCGAGCCTGGCCTGGTACTCGACCTGCCATGCGGAGCGGGCCGTTTCTGGCCGCTGTTGGCGGAAAAGCCCAATCGGGTA  >  NC_002947/1043544‑1043839
                                                                                                                                                   |                                                                                                                                                    
acaTTGGCATGACCTGATTCCAACGCTTGCATGACTAGTGTGGGCGCAATGGAGTCAGGAGATGCCCAATGCGTAGCCCCATCAAGCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACg                                                                                                                                                    <  1:1122099/150‑1 (MQ=255)
  aTTGGCATGACCTGATTCCAACGCTTGCATGACTAGTGTGGGCGCAATGGAGACAGGAGATGCCCAATGCGTAGCCCCATGAAGCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAATCACCAGGACGGc                                                                                                                                                  >  1:678354/1‑150 (MQ=255)
   ttGGCATGACCTGATTCCAACGCTTGCATGACTAGTGTGGGCGCAATGGAGTCAGGAGATGCCCAATGCGTAGCCCCATCAAGCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGcc                                                                                                                                                 >  2:830268/1‑150 (MQ=255)
      gCATGACCTGATTCCAACGCTTGCATGACTAGTGTGGGCGCAATGGAGTCAGGAGATGCCCAATGCGTAGCCCCATCAAGCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTgg                                                                                                                                              <  1:385519/150‑1 (MQ=255)
          gACCTGATTCCAACGCTTGCATGACTAGTGTGGGCGCAATGGAGTCAGGAGATGCCCAATGCGTAGCCCCATCAAGCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGCCGTCCTGgtgct                                                                                                                                          <  1:363341/150‑2 (MQ=255)
          gACCTGATTCCAACGCTTGCATGACTAGTGTGGGCGCAATGGAGTCAGGAGATGCCCAATGCGTAGCCCCATCAAGCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGCCGTCCTGgtgct                                                                                                                                          <  1:363262/150‑2 (MQ=255)
              tGATTCCAACGCTTGCATGACTAGTGTGGGCGCAATGGAGTCAGGAGATGCCCAATGCGTAGCCCCATCAAGCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGc                                                                                                                                      >  1:102930/1‑150 (MQ=255)
                    cAACGCTTGCATGACTAGTGTGGGCGCAATGGAGTCAGGAGATGCCCAATGCGTAGCCCCATCAAGCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTcc                                                                                                                                >  2:711251/1‑150 (MQ=255)
                    cAACGCTTGCATGACTAGTGTGGGCGCAATGGAGTCAGGAGATGCCCAATGCGTAGCCCCATCAAGCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTcc                                                                                                                                >  2:367040/1‑150 (MQ=255)
                          ttacATGACTAGTGTGGGCGCAATGGAGTCAGGAGATGCCCAATGCGTAGCCCCATCAAGCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAg                                                                                                                          >  1:1206059/4‑150 (MQ=255)
                            gCATGACTAGTGTGGGCGCAATGGAGTCAGGAGATGCCCAATGCGTAGCCCCATCAAGCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCg                                                                                                                        >  2:1090603/1‑150 (MQ=255)
                             cATGACTAGTGTGGGCGCAATGGAGTCAGGAGATGCCCAATGCGTAGCCCCATCAAGCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCga                                                                                                                       <  2:573180/150‑1 (MQ=255)
                             cATGACTAGTGTGGGCGCAATGGAGTCAGGAGATGCCCAATGCGTAGCCCCATCAAGCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCga                                                                                                                       <  1:912169/150‑1 (MQ=255)
                              aTGACTAGTGTGGGCGCAATGGAGTCAGGAGATGCCCAATGCGTAGCCCCATCAAGCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCgag                                                                                                                      <  2:1206059/150‑1 (MQ=255)
                              aTGACTAGTGTGGGCGCAATGGAGTCAGGAGATGCCCAATGCGTAGCCCCATCAAGCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCgag                                                                                                                      <  2:910570/150‑1 (MQ=255)
                              aTGACTAGTGTGGGCGCAATGGAGTCAGGAGATGCCCAATGCGTAGCCCCATCAAGCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCgag                                                                                                                      >  1:738423/1‑150 (MQ=255)
                                     gtgtGGGCGCAATGGAGTCAGGAGATGCCCAATGCGTAGCCCCATCAAGCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGa                                                                                                                                         <  2:497045/124‑1 (MQ=255)
                                     gtgtGGGCGCAATGGAGTCAGGAGATGCCCAATGCGTAGCCCCATCAAGCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGa                                                                                                                                         >  1:497045/1‑124 (MQ=255)
                                                       cAGGAGATGCCCAATGCGTAGCCCCATCAAGCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCtt                                                                                             <  1:1059861/150‑1 (MQ=255)
                                                         ggAGATGCCCAATGCGTAGCCCCATCAAGCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCgcggcgcgccaact                                                                                           <  1:1014964/150‑5 (MQ=255)
                                                            gaTGCCCAATGCGTAGCCCCATCAAGCTTGAATTTTCCGAGAAGTACACCGAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCAt                                                                                        <  1:676151/150‑1 (MQ=255)
                                                                cccAATGCGTAGCCCCATCAAGCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCATTggc                                                                                    >  2:685828/1‑150 (MQ=255)
                                                                 ccAATGCGTAGCCCCATCAAGCTTGAATTTTCCGAGAAGTACGACAAAGCCCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGAAGCCTGTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCATTggcg                                                                                   >  1:564871/1‑150 (MQ=255)
                                                                  cAATGCGTAGCCCCATCAAGCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGGCCGGTTGGCGCGCCGCGCCTTGGCATTggccg                                                                                  >  1:907539/1‑148 (MQ=255)
                                                                    aTGCGTAGCCCCATCAAGCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCATTggcgggc                                                                                <  1:685828/150‑1 (MQ=255)
                                                                      gCGTAGCCCCATCAAGCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCATTGGCGGGCGa                                                                              <  1:711251/150‑1 (MQ=255)
                                                                        gTAGCCCCATCAAGCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCgcggcgcgccaa                                                                                             >  2:292002/1‑131 (MQ=255)
                                                                        gTAGCCCCATCAAGCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCgcggcgcgccaa                                                                                             <  1:292002/133‑3 (MQ=255)
                                                                          aGCCCCATCAAGCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCATTGGCGGGCGAgcct                                                                          <  2:566654/150‑1 (MQ=255)
                                                                            ccccATCAAGCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCATTGGCGGGCGAgcctgg                                                                        >  2:1033418/1‑150 (MQ=255)
                                                                             cccATCAAGCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCATTGGCGGGCGAgcctggc                                                                       >  1:89427/1‑150 (MQ=255)
                                                                                     gCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAg                                                                                                              <  1:601494/103‑1 (MQ=255)
                                                                                     gCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAg                                                                                                              >  2:601494/1‑103 (MQ=255)
                                                                                         gAATTTTCCGAGAAGTCCGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCATTGGCGGGCGAGCCTGGCCTGGTACTCGAc                                                           <  1:367040/150‑1 (MQ=255)
                                                                                            ttttCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCtt                                                                                             <  1:752124/113‑1 (MQ=255)
                                                                                            ttttCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCtt                                                                                             >  2:752124/1‑113 (MQ=255)
                                                                                                  gagaAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCATTGGCGGGCGAGCCTGGCCTGGTACTCGACCTGCCATGc                                                  >  2:1115595/1‑150 (MQ=255)
                                                                                                         aCGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCATTGGCGGGCGAGCCTGGCCTGGTACTCGACCTGCCATGCGGAGCgg                                           <  1:960047/150‑1 (MQ=255)
                                                                                                           gACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGcc                                                                                                                                                 <  1:816700/46‑1 (MQ=255)
                                                                                                           gACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGcc                                                                                                                                                 >  2:816700/1‑46 (MQ=255)
                                                                                                            aCAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCATTGGCGGGCGAGCCTGGCCTGGTACTCGACCTGCCATGCGGAGCGGGcc                                        >  1:914435/1‑150 (MQ=255)
                                                                                                              aaaGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGGCTACGCATTGGCGGGCGAGCCTGGCCTGGTACTCGACCTGCCATGCGGAGCGGGCCGt                                      >  1:616793/1‑150 (MQ=255)
                                                                                                              aaaGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCATTGGCGGGCGAGCCTGGCCTGGTACTCGACCTGCCATGCGGAGCGGGCCGt                                      <  2:102930/150‑1 (MQ=255)
                                                                                                                   ccATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCATTGGCGGGCGAGCCTGGCCTGGTACTCGACCTGCCATGCGGAGCGGGCCGTTTCTg                                 <  1:1115595/150‑1 (MQ=255)
                                                                                                                     aTGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCATTGGCGGGCGAGCCTGGCCTGGTACTCGACCTGCCATGCGGAGCGGGCCGTTTCTGGc                               <  2:914435/150‑1 (MQ=255)
                                                                                                                           gcgAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCATTGGCGGGCGAGCCTGGCCTGGTACTCGACCTGCCATGCGGAGCGGGCCGTTTCTGGCCGCTGt                         <  1:1090603/150‑1 (MQ=255)
                                                                                                                             gAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCATTGGCGGGCGAGCCTGGCCTGGTACTCGACCTGCCATGCGGAGCGGGCCGTTTCTGGCCGCTGTTg                       >  1:226792/1‑150 (MQ=255)
                                                                                                                              cttccgatctcAAGCACCAGGACGGCCTGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCATTGGCGGGCGAGCCTGGCCTGGTACTCGACCTGCCATGTGCAGCGGGCCGTTTCTGGCCGCTGTTGGc                     <  1:939140/143‑1 (MQ=255)
                                                                                                                                tACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCATTGGCGGGCGAGCCTGGCCTGGTACTCGACCTGCCATGCGGAGCGGGCCGTTTCTGGCCGCTGTTGGCg                    <  2:226792/150‑1 (MQ=255)
                                                                                                                                   tttCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCATTGGCGGGCGAGCCTGGCCTGGTACTCGACCTGCCATGCg                                                 <  1:680747/118‑1 (MQ=255)
                                                                                                                                   tttCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCATTGGCGGGCGAGCCTGGCCTGGTACTCGACCTGCCATGCg                                                 >  2:680747/1‑118 (MQ=255)
                                                                                                                                   tttCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCATTGGCGGGCGAGCCTGGCCTGGTACTCGACCTGCCATGCGGAGCGGGCCGTTTCTGGCCGCTGTTGGCGGaa                 >  2:73525/1‑150 (MQ=255)
                                                                                                                                        gccgcgccAGGCCGTCCTGGTGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCATTGGCGGGCGAGCCTGGCCTGGTACTCGACCTGCCATGCGGAGCGGGCCGTTTCTGGCCGCTGTTGGCGGAAAAGcc            >  1:1136944/7‑150 (MQ=255)
                                                                                                                                         tcgcgccAGGCCGTCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCATTGGCGGGCGAGCCTGGCCTGGTCCTCGACCTGCCATGCGGAGCGGGCCGTTTCTGGCCGCTGTTGGCGGAAAAGccc           >  1:128264/6‑150 (MQ=255)
                                                                                                                                         aaGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCATTGGCGGGCGAGCCTGGCCTGGTACTCGACCTGCCATGCGGAGCGGGCCGttt                                    >  1:1185896/1‑125 (MQ=255)
                                                                                                                                         aaGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCATTGGCGGGCGAGCCTGGCCTGGTACTCGACCTGCCATGCGGAGCGGGCCGttt                                    <  2:1185896/125‑1 (MQ=255)
                                                                                                                                         aaGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCATTGGCGGGCGAGCCTGGCCTGGTACTCGACCTGCCATGCGGAGCGGGCCGTTTCTGGCCGCTGTTGGCGGAAAAGccc           <  1:145370/150‑1 (MQ=255)
                                                                                                                                             aCCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCATTGGCGGGCGAGCCTGGCCTGGTACTCGACCTGCCATGCGGAGCGGGCCGTTTCTGGCCGCTGTTGGCGGAAAAGCCCAATc       >  2:725786/1‑150 (MQ=255)
                                                                                                                                               cAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCATTGGCGGGCGAGCCTGGCCTGGTACTCGACCTGCCATGCGGAGCGGGCCGTTTCTGGCCGCTGTTGGCGGAAAAGCCCAATCgg     <  2:593257/150‑1 (MQ=255)
                                                                                                                                               cAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCATTGGCGGGCGAGCCTGGCCTGGTACTCGACCTGCCATGCGGAGCGGGCCGTTTCTGGCCGCTGTTGGCGGAAAAGCCCAATCgg     >  1:593257/1‑150 (MQ=255)
                                                                                                                                                aGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCATTGGCGGGCGAGCCTGGCCTGGTACTCGACCTGCCATGCGGAGCGGGCCGTTTCTGGCCGCTGTTGGCGGAAAAGCCCAATCggg    <  2:89427/150‑1 (MQ=255)
                                                                                                                                                  aaCGGCCTGGCGCGACGACTTTCCCACAAGCGAGATGAGCAGTTGGCGCGGCGCGCCTTGGCATTGGCGGGCGAGGCTGGCCTGGTACTCGACCTGCCATGCGGAGAGGGCCGTTTCTGGCCGCTGTTGGCGGAACAGCCCAATCGGGTa  <  2:678354/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ACATTGGCATGACCTGATTCCAACGCTTGCATGACTAGTGTGGGCGCAATGGAGTCAGGAGATGCCCAATGCGTAGCCCCATCAAGCTTGAATTTTCCGAGAAGTACGACAAAGACCATGCGCGCGAGTACTTTCTCAAGCACCAGGACGGCCTGGCGCGACGCCTTTCCCACAAGCGAGATGAGCAGTTGGCGCGCCGCGCCTTGGCATTGGCGGGCGAGCCTGGCCTGGTACTCGACCTGCCATGCGGAGCGGGCCGTTTCTGGCCGCTGTTGGCGGAAAAGCCCAATCGGGTA  >  NC_002947/1043544‑1043839

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