Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F3 I160 R1 220 104.1 6051714 95.9% 5803593 84.4

Breseq alignment

BRESEQ :: Evidence
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: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCGGAATTGCCGCCATCGCC  >  NZ_CP009273/4108731‑4108929
                                                                                                   |                                                                                                   
ACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGG                                                                                                     <  SRR3722030.952486/100‑1 (MQ=60)
ACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGG                                                                                                     <  SRR3722030.2737322/100‑1 (MQ=60)
ACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGG                                                                                                     <  SRR3722030.278007/100‑1 (MQ=60)
ACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGG                                                                                                     <  SRR3722030.846983/100‑1 (MQ=60)
ACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGG                                                                                                     <  SRR3722030.66407/100‑1 (MQ=60)
ACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGG                                                                                                     <  SRR3722030.364257/100‑1 (MQ=60)
ACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGG                                                                                                     <  SRR3722030.3027714/100‑1 (MQ=60)
  AGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTG                                                                                                   >  SRR3722030.200457/1‑100 (MQ=60)
        CGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCC                                                                                             >  SRR3722030.1674652/1‑100 (MQ=60)
            GAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAG                                                                                         >  SRR3722030.3002781/1‑100 (MQ=60)
               AAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAG                                                                                      >  SRR3722030.26954/1‑100 (MQ=60)
                    gaccaccgagatctacactggccatgtcgtcggcagcgtcagatgtgtataagagacagCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGT                                                                                 <  SRR3722030.2103928/41‑1 (MQ=60)
                    GCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGT                                                                                 <  SRR3722030.742600/100‑1 (MQ=60)
                      ACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGA                                                                               >  SRR3722030.1218238/1‑100 (MQ=60)
                      ACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGA                                                                               >  SRR3722030.1878486/1‑100 (MQ=60)
                       CCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGAC                                                                              <  SRR3722030.995907/100‑1 (MQ=60)
                         TTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACAC                                                                            <  SRR3722030.722669/100‑1 (MQ=60)
                         TTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACAC                                                                            <  SRR3722030.1324189/100‑1 (MQ=60)
                         TTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACAC                                                                            <  SRR3722030.2406690/100‑1 (MQ=60)
                              GTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCA                                                                       >  SRR3722030.948823/1‑100 (MQ=60)
                              GTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCA                                                                       >  SRR3722030.2696144/1‑100 (MQ=60)
                              GTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAAGAAACACCGCCA                                                                       >  SRR3722030.2740777/1‑100 (MQ=60)
                                GCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAA                                                                     >  SRR3722030.2233665/1‑100 (MQ=60)
                                            ATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGAC                                                         >  SRR3722030.2975044/1‑100 (MQ=60)
                                            ATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGAC                                                         >  SRR3722030.2076174/1‑100 (MQ=60)
                                             TTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCctgtctcttatacaca                                                        >  SRR3722030.1362455/1‑84 (MQ=60)
                                             TTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACA                                                        >  SRR3722030.2027754/1‑100 (MQ=60)
                                             TTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACA                                                        >  SRR3722030.119722/1‑100 (MQ=60)
                                              TCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGTGGACAT                                                       >  SRR3722030.1082028/1‑100 (MQ=60)
                                                        GCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCC                                             <  SRR3722030.1319167/100‑1 (MQ=60)
                                                           GATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGctgtctcttatacacatctgacgctgccgacgacatggcca                                          >  SRR3722030.190661/1‑59 (MQ=60)
                                                             TCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGC                                        <  SRR3722030.3050414/100‑1 (MQ=60)
                                                              CTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCC                                       <  SRR3722030.740752/100‑1 (MQ=60)
                                                              CTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCC                                       <  SRR3722030.1417418/100‑1 (MQ=60)
                                                                   GAAAGAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGC                                  >  SRR3722030.933541/1‑100 (MQ=60)
                                                                       GAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGC                              >  SRR3722030.1364899/1‑100 (MQ=60)
                                                                       GAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGC                              >  SRR3722030.2460016/1‑100 (MQ=60)
                                                                       GAATAATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGC                              >  SRR3722030.2565923/1‑100 (MQ=60)
                                                                            ATCCCGGTATTGTCGGCATAAAGTTGGCCGCCGGAGAAGGAGGTGACACGCCCAAAATTGGAGCGGACATCCCTTTGGCCAATGCCTTCGGCAGGGGCAC                         >  SRR3722030.2875326/1‑100 (MQ=42)
                                                                            ATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCAC                         <  SRR3722030.2368354/100‑1 (MQ=60)
                                                                            ATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCAC                         >  SRR3722030.2135016/1‑100 (MQ=60)
                                                                            ATCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCAC                         >  SRR3722030.2145613/1‑100 (MQ=60)
                                                                             TCCCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACC                        >  SRR3722030.2887436/1‑100 (MQ=60)
                                                                               CCGGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCctgtctcttatacacatctgacgctgccga                      >  SRR3722030.410880/1‑70 (MQ=60)
                                                                                 taagagacaGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACC                    <  SRR3722030.2431831/91‑1 (MQ=60)
                                                                                 GGTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACC                    <  SRR3722030.1624799/100‑1 (MQ=60)
                                                                                 GGTATTGTCGGCATAAAGGCGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCAGCACCAACC                    >  SRR3722030.1267780/1‑100 (MQ=60)
                                                                                  GTATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCG                   >  SRR3722030.269175/1‑100 (MQ=60)
                                                                                   TATTGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCGG                  <  SRR3722030.2832556/100‑1 (MQ=60)
                                                                                    gagacaGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCGGA                 <  SRR3722030.2891863/94‑1 (MQ=60)
                                                                                      TGTCGGCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCGGAAT               <  SRR3722030.1124282/100‑1 (MQ=60)
                                                                                           GCATAAAGGTGGTCGCCGGAGAAGAAGGTGGCACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCGGAATTGCCG          >  SRR3722030.937947/1‑100 (MQ=60)
                                                                                           GCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCGGAATTGCCG          >  SRR3722030.2222912/1‑100 (MQ=60)
                                                                                           GCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCGGAATTGCCG          >  SRR3722030.1544629/1‑100 (MQ=60)
                                                                                           GCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCGGAATTGCCG          >  SRR3722030.1853556/1‑100 (MQ=60)
                                                                                           GCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCGGAATTGCCG          >  SRR3722030.492845/1‑100 (MQ=60)
                                                                                           GCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCGGAATTGCCG          >  SRR3722030.1443355/1‑100 (MQ=60)
                                                                                           GCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCGGAATTGCCG          >  SRR3722030.2346754/1‑100 (MQ=60)
                                                                                           GCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCGGAATTGCCG          >  SRR3722030.2138203/1‑100 (MQ=60)
                                                                                           GCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCGGAATTGCCG          >  SRR3722030.239000/1‑100 (MQ=60)
                                                                                           GCATAAAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCGGAATTGCCG          >  SRR3722030.1068251/1‑100 (MQ=60)
                                                                                            CATAAAGGTGGTCGCCGGAGAAAAAGGTGACCCCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCGGAATTGCCGC         <  SRR3722030.2763369/100‑1 (MQ=60)
                                                                                                AAGGTGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCGGAATTGCCGCCATC     >  SRR3722030.2409429/1‑100 (MQ=60)
                                                                                                   GTGGTCGCCGGAGAAGAAGGTGCCACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCGGAATTGCCGCCATCGCC  <  SRR3722030.3024104/100‑1 (MQ=60)
                                                                                                   |                                                                                                   
ACAGGTAACGCAGAAAAAAGGCACCTTGCGGTGCCTTTCTTATCATTCAATATCCAGCGGATCTTCTGAAAGAATAATCCCGGTATTGTCGGCATAAAGATGGTCGCCGGAGAAGAAGGTGACACCGCCAAAATTGACGCGGACATCGCTTTCGCCAATGCCTTCGCCAGCGGCACCAACCGGAATTGCCGCCATCGCC  >  NZ_CP009273/4108731‑4108929

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