Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,289,430 T→G 9.1% K1053Q (AAG→CAG)  PP_3761 ← two‑component system sensor histidine kinase/response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,289,4300TG 9.1% 195.2 / 8.1 66K1053Q (AAG→CAG) PP_3761two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base T (32/28);  new base G (3/3);  total (35/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 = 9.73e-01

GCGGATCAGGCGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCGATTCGACCTT  >  NC_002947/4289293‑4289574
                                                                                                                                         |                                                                                                                                                
gCGGATCAGGCGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAggg                                                                                                                                       >  4:20476/1‑149 (MQ=255)
gCGGATCAGGCGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAggg                                                                                                                                       >  6:30418/1‑149 (MQ=255)
gCGGATCAGGCGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAggg                                                                                                                                       <  2:90636/149‑1 (MQ=255)
     tCAGGCGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTg                                                                                                                                  <  4:390083/149‑1 (MQ=255)
       aGGCGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTGGAGCACAAGGgcgccatt                                                                                                                                <  8:56325/149‑5 (MQ=255)
        ggCGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGt                                                                                                                                >  7:15397/1‑148 (MQ=255)
         gCGGGTGGCCTCGTAGCCGGCCATCGCCGGCATGCTCCCGGCCATCAGCCCCAGGGCGATATCGGCGTGCTGCTCCCGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCa                                                                                                                              <  7:71704/149‑1 (MQ=255)
          cGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAg                                                                                                                             <  5:148293/149‑1 (MQ=255)
           gggTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAgg                                                                                                                            <  2:299315/149‑1 (MQ=255)
            ggTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAggg                                                                                                                           <  4:93038/149‑1 (MQ=255)
                gCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCaaa                                                                                                                       <  8:293793/149‑1 (MQ=255)
                   tCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGa                                                                                                                     >  5:292825/1‑148 (MQ=255)
                      aaGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGCCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATAtt                                                                                                                 <  8:260061/149‑1 (MQ=255)
                        gCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTAc                                                                                                               <  8:15397/149‑1 (MQ=255)
                           gTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGca                                                                                                            >  7:262014/1‑149 (MQ=255)
                              cATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACAt                                                                                                         >  1:306551/1‑149 (MQ=255)
                                 ctcCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAg                                                                                                                             <  3:193616/126‑1 (MQ=255)
                                 ctcCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAg                                                                                                                             >  4:193616/1‑126 (MQ=255)
                                    cGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACAtcgtcg                                                                                                    <  8:262014/148‑1 (MQ=255)
                                    cGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACAtcgtc                                                                                                     <  3:20476/147‑1 (MQ=255)
                                       catcatCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCGGGTCAGGGCAAAGATATTACGCACATCGTCGTcca                                                                                                <  3:84959/149‑1 (MQ=255)
                                        atcatcACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTccac                                                                                               <  6:292825/149‑1 (MQ=255)
                                         tcatcaCGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTccac                                                                                               >  2:306056/1‑148 (MQ=255)
                                                    aTCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCtt                                                                                   >  2:137492/1‑149 (MQ=255)
                                                    aTCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTGGAGCACAAGGgcgc                                                                                                                                    <  6:129118/100‑1 (MQ=255)
                                                    aTCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTGGAGCACAAGGgcgc                                                                                                                                    >  5:129118/1‑100 (MQ=255)
                                                       aGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGGTTCGCACCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCg                                                                                >  1:273224/1‑149 (MQ=255)
                                                          aCCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCT‑‑AGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCt                                                                           >  6:44641/1‑149 (MQ=255)
                                                           ccAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATGTCGACTATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGGCCACCAGCAGCACCGTGCGGcgc                                                                            >  4:27821/1‑147 (MQ=255)
                                                              ggTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCg                                                                         >  2:141197/1‑149 (MQ=255)
                                                                 cGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCaaa                                                                                                                       <  6:66737/100‑1 (MQ=255)
                                                                 cGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCaaa                                                                                                                       >  5:66737/1‑100 (MQ=255)
                                                                 cGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGaa                                                                       >  3:198055/1‑148 (MQ=255)
                                                                  gATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAaca                                                                     >  3:387808/1‑149 (MQ=255)
                                                                   atatCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGaa                                                                       <  4:209174/146‑1 (MQ=255)
                                                                   atatCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGaa                                                                       >  3:209174/1‑146 (MQ=255)
                                                                   atatCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCGGGTCAggg                                                                                                                           <  2:220985/94‑1 (MQ=255)
                                                                   aaatCGTCGTCCTGCTCCAGGCGCTCCAAGGCATCGCGCACGTTGAGGCCGATTTAGACAATGGCGCCCTTGTGCTCCAGGGCTATGGTCAggg                                                                                                                           >  1:220985/3‑94 (MQ=255)
                                                                       cgtcgtGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTg                                                                >  6:338208/1‑149 (MQ=255)
                                                                        gtcgtGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTg                                                                <  1:306056/148‑1 (MQ=255)
                                                                         tcgtGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTc                                                              <  8:324014/149‑1 (MQ=255)
                                                                           gtGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTcgc                                                            >  6:52016/1‑149 (MQ=255)
                                                                                  tCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTAc                                                                                                               >  4:179411/1‑91 (MQ=255)
                                                                                  tCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTAc                                                                                                               <  3:179411/91‑1 (MQ=255)
                                                                                   ccAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTgcgc                                                    >  2:104066/1‑149 (MQ=255)
                                                                                    cAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTgcgcg                                                   >  7:199196/1‑149 (MQ=255)
                                                                                       gcgcTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCgg                                                <  5:30418/149‑1 (MQ=255)
                                                                                        cgcTCGATGGCTTCGCGCACGATGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCAAGGGCGATGGTCAGGGCAAAGATATTAAGCACATCGACGTCCACAACCAGCACCTTGCGGCCCTAGAACACCTTGTCCCGGCTGCGCGCgg                                                >  6:208281/1‑148 (MQ=255)
                                                                                             gATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGCACACCTTGTCGCGGCTGCGCGCGGTCtt                                            >  7:77261/1‑147 (MQ=255)
                                                                                              aTGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAg                                         >  7:194076/1‑149 (MQ=255)
                                                                                                 gCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCAt                                      <  4:198055/149‑1 (MQ=255)
                                                                                                     cgcgCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGc                                    <  4:387808/147‑1 (MQ=255)
                                                                                                            gTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGAcgtt                           <  5:208281/149‑1 (MQ=255)
                                                                                                            gTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCT‑‑AGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGAcgttcg                         <  5:44641/149‑1 (MQ=255)
                                                                                                                cGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGAcgttcgtt                       >  7:312323/1‑149 (MQ=255)
                                                                                                                 ggCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTg                      <  5:338208/149‑1 (MQ=255)
                                                                                                                  gCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTg                      >  3:244693/1‑148 (MQ=255)
                                                                                                                  gCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCCCGGCTGCGCGCGGTCTTGAGCATGCGTTAACGTACGTTg                      >  2:88974/1‑148 (MQ=255)
                                                                                                                   ccGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGa                    <  1:221504/149‑1 (MQ=255)
                                                                                                                        ttCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGCCAGCt               >  7:116423/1‑149 (MQ=255)
                                                                                                                                tGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGAcgttc                          <  4:83533/130‑1 (MQ=255)
                                                                                                                                tGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGAcgttc                          >  3:83533/1‑130 (MQ=255)
                                                                                                                                cagcgcCCTGGAGCACAAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCGATTc        >  8:67565/3‑148 (MQ=255)
                                                                                                                                 agcgcCCTGGAGCACAAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTGGGTTGGACAGCTGCGATTCGa      >  2:98075/2‑149 (MQ=255)
                                                                                                                                  gcgcCCTGGAGCACAAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCGATTCGAc     <  7:67565/149‑1 (MQ=255)
                                                                                                                                     cccTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCGATTCGACCtt  >  1:28802/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                                
GCGGATCAGGCGGGTGGCCTCGAAGCCGTCCATCTCCGGCATCATCACGTCCATCAGCACCAGGTCGATATCGTCGTGCTGCTCCAGGCGCTCGATGGCTTCGCGCCCGTTGCGGCCGATTTCGACAATGGCGCCCTTGTGCTCCAGGGCGCTGGTCAGGGCAAAGATATTACGCACATCGTCGTCCACCAGCAGCACCTTGCGGCCCTCGAACACCTTGTCGCGGCTGCGCGCGGTCTTGAGCATGCGTTGACGTTCGTTGGACAGCTGCGATTCGACCTT  >  NC_002947/4289293‑4289574

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