Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,108,830 A→G H143H (CAT→CAC rraA ← ribonuclease E activity regulator RraA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,108,8300AG99.0% 368.8 / NA 100H143H (CAT→CACrraAribonuclease E activity regulator RraA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (57/42);  minor base T (1/0);  total (58/42)
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.00e+00

GAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCG  >  NZ_CP009273/4108743‑4108912
                                                                                       |                                                                                  
gAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTg                                                                                  >  1:197627/1‑90 (MQ=255)
 aaaaaaGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTgg                                                                                 >  2:636440/1‑90 (MQ=255)
   aaaaGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTc                                                                               <  2:844827/90‑1 (MQ=255)
   aaaaGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTc                                                                               <  2:2395428/90‑1 (MQ=255)
      aGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGcc                                                                            >  1:1647906/1‑90 (MQ=255)
        gCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCgg                                                                          <  1:731312/90‑1 (MQ=255)
          aCCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGgag                                                                        >  1:2948411/1‑90 (MQ=255)
           ccTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGgaga                                                                       <  2:1046116/90‑1 (MQ=255)
           ccTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGgaga                                                                       <  1:980988/90‑1 (MQ=255)
           ccTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGgaga                                                                       <  2:1128676/90‑1 (MQ=255)
            cTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGagaa                                                                      >  2:1802156/1‑90 (MQ=255)
             ttGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTTGCCGGagaag                                                                     >  2:816043/1‑90 (MQ=255)
             ttGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGgaga                                                                       >  1:26544/1‑88 (MQ=255)
             ttGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGgaga                                                                       <  2:26544/88‑1 (MQ=255)
             ttGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGagaag                                                                     <  1:1304354/90‑1 (MQ=255)
             ttGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGagaag                                                                     <  1:711672/90‑1 (MQ=255)
             ttGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGagaag                                                                     <  1:2363711/90‑1 (MQ=255)
             ttGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGagaag                                                                     <  2:2697373/90‑1 (MQ=255)
                    gCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGa                                                              >  1:1847010/1‑90 (MQ=255)
                    gCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGa                                                              >  1:1200479/1‑90 (MQ=255)
                    gCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGa                                                              >  2:2949458/1‑90 (MQ=255)
                    gCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGa                                                              >  2:2688804/1‑90 (MQ=255)
                     ccTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGac                                                             <  2:1351094/90‑1 (MQ=255)
                     ccTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGac                                                             >  2:2857248/1‑90 (MQ=255)
                       tttCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGacac                                                           <  2:760707/90‑1 (MQ=255)
                         tCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCg                                                         >  2:3019872/1‑90 (MQ=255)
                            tATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCa                                                      >  1:2647409/1‑90 (MQ=255)
                            tATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCa                                                      >  1:934526/1‑90 (MQ=255)
                            tATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAAGAAACACCGCCa                                                      >  1:2691110/1‑90 (MQ=255)
                              tCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCaaa                                                    >  2:1626807/1‑90 (MQ=255)
                              tCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCaaa                                                    >  1:2194224/1‑90 (MQ=255)
                                 ttCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAAtt                                                 >  2:783610/1‑90 (MQ=255)
                                   cAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGGGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGa                                               >  2:2382154/1‑90 (MQ=255)
                                     atatCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGAcg                                             <  2:12885/90‑1 (MQ=255)
                                      tatCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGAcgc                                            <  2:2194224/90‑1 (MQ=255)
                                         ccAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGa                                         >  2:1270121/1‑90 (MQ=255)
                                          cAGCGGTTCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGAc                                        >  2:1796871/1‑90 (MQ=255)
                                          cAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGAc                                        >  2:2353141/1‑90 (MQ=255)
                                          cAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGAc                                        >  1:2040069/1‑90 (MQ=255)
                                          cAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGAc                                        >  1:2921181/1‑90 (MQ=255)
                                           aGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAgg                                                                 <  2:1341728/64‑1 (MQ=255)
                                           aGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAgg                                                                 >  1:1341728/1‑64 (MQ=255)
                                           aGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACa                                       >  2:1934195/1‑90 (MQ=255)
                                           aGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACa                                       >  1:1992752/1‑90 (MQ=255)
                                           aGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACa                                       >  2:1025427/1‑90 (MQ=255)
                                           aGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACa                                       >  1:118004/1‑90 (MQ=255)
                                            gCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGTGGACAt                                      >  1:1066033/1‑90 (MQ=255)
                                            gCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACAt                                      <  2:1699305/90‑1 (MQ=255)
                                            gCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACAt                                      <  2:2305653/90‑1 (MQ=255)
                                            gCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACAt                                      >  2:1140223/1‑90 (MQ=255)
                                            gCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACAt                                      <  1:1299456/90‑1 (MQ=255)
                                            gCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACAt                                      <  2:1032273/90‑1 (MQ=255)
                                            gCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACAt                                      <  2:2040069/90‑1 (MQ=255)
                                               gATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGc                                   >  2:1996251/1‑90 (MQ=255)
                                               gATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGc                                   >  2:1638929/1‑90 (MQ=255)
                                               gATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGc                                   >  2:912493/1‑90 (MQ=255)
                                                 tcttctGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCtt                                 <  2:2871636/90‑1 (MQ=255)
                                                 tcttctGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCtt                                 >  2:687594/1‑90 (MQ=255)
                                                 tcttctGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCtt                                 <  2:2317570/90‑1 (MQ=255)
                                                 tcttctGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCtt                                 <  1:2995241/90‑1 (MQ=255)
                                                  cttctGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCttt                                <  1:1395469/90‑1 (MQ=255)
                                                  cttctGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCttt                                <  1:729498/90‑1 (MQ=255)
                                                      tGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGcc                            >  2:538104/1‑90 (MQ=255)
                                                        aaaGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGGGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCaa                          <  2:158790/90‑1 (MQ=255)
                                                        aaaGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCaa                          <  2:1883693/90‑1 (MQ=255)
                                                        aaaGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCaa                          <  2:2041172/90‑1 (MQ=255)
                                                         aaGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGcc                                                                            <  2:187963/39‑1 (MQ=38)
                                                         aaGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGcc                                                                            >  1:187963/1‑39 (MQ=38)
                                                         aaGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAAt                         >  2:2981303/1‑90 (MQ=255)
                                                            aataatCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGcc                      >  2:1733483/1‑90 (MQ=255)
                                                             ataatCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCt                     >  2:1855928/1‑90 (MQ=255)
                                                                atCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCg                  <  1:2326137/90‑1 (MQ=255)
                                                                 tCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGc                 >  1:919452/1‑90 (MQ=255)
                                                                     ggTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACGCCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGc             >  2:2981342/1‑90 (MQ=255)
                                                                     ggTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGc             >  2:209880/1‑90 (MQ=255)
                                                                     ggTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGc             >  1:1344107/1‑90 (MQ=255)
                                                                     ggTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGc             >  2:900949/1‑90 (MQ=255)
                                                                     ggTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGc             >  1:2415976/1‑90 (MQ=255)
                                                                     ggTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGc             <  1:1598863/90‑1 (MQ=255)
                                                                     ggTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGc             >  2:1314364/1‑90 (MQ=255)
                                                                     ggTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGc             >  1:2519786/1‑90 (MQ=255)
                                                                       tATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCgg           <  1:2781173/90‑1 (MQ=255)
                                                                         ttGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCa         <  2:1280403/90‑1 (MQ=255)
                                                                          tGTCGGCATAAAGTTGGCCGCCGGAGAAGGAGGTGACACGCCCAAAATTGGAGCGGACATCCCTTTGGCCAATGCCTTCGGCAGGGGCAc        >  1:2823216/1‑90 (MQ=255)
                                                                          tGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCAc        <  2:2256029/90‑1 (MQ=255)
                                                                          tGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCAc        <  2:1382053/90‑1 (MQ=255)
                                                                          tGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCAc        >  1:2108004/1‑90 (MQ=255)
                                                                          tGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCAc        <  2:1066033/90‑1 (MQ=255)
                                                                          tGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCAc        <  1:1107746/90‑1 (MQ=255)
                                                                          tGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCAc        >  1:2097647/1‑90 (MQ=255)
                                                                           gTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCAcc       >  2:1399889/1‑90 (MQ=255)
                                                                           gTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCAcc       >  1:2835097/1‑90 (MQ=255)
                                                                             cGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGAcg                                             <  2:404781/50‑1 (MQ=255)
                                                                             cGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGAcg                                             >  1:404781/1‑50 (MQ=255)
                                                                               gCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGTGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAAcc   <  2:2677833/90‑1 (MQ=255)
                                                                               gCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAAcc   <  2:200475/90‑1 (MQ=255)
                                                                               gCATAAAGGCGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCAGCACCAAcc   >  1:1249269/1‑90 (MQ=255)
                                                                                cATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCg  >  2:2783705/1‑90 (MQ=255)
                                                                                cATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCg  >  1:265324/1‑90 (MQ=255)
                                                                                cATAAAGGTGGTCGCCGGAGAAAAAGGTGACCCCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCg  <  1:2713244/90‑1 (MQ=255)
                                                                                       |                                                                                  
GAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCG  >  NZ_CP009273/4108743‑4108912

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