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
A2 F5 I175 R1 328 64.7 2582366 96.2% 2484236 120.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 883,821 C→T V485V (GTG→GTA ybjL ← aspartate:alanine antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273883,8210CT97.0% 222.1 / ‑2.3 68V485V (GTG→GTAybjLaspartate:alanine antiporter
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (32/33);  total (33/35)
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.18e-01

GGATGTTACTGCGAGCTGTATCACTGATGATCTCCATTGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGCACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGAACTCTTTCACCATGC  >  NZ_CP009273/883685‑883951
                                                                                                                                        |                                                                                                                                  
ggATGTTACTGCGAGCTGTATCACTGATGATCTCCATTGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTAc                                                                                                                                  >  1:969348/1‑139 (MQ=255)
  aTGTTACTGCGAGCTGTATCACTGATGATCTCCATTGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCTCGGGTACCa                                                                                                                                <  2:42938/139‑1 (MQ=255)
     ttACTGCGAGCTGTATCACTGATGATCTCCATTGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGa                                                                                                                              <  1:1279608/138‑1 (MQ=255)
           cGAGCTGTATCACTGATGATCTCCATTGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGAcaa                                                                                                                       >  2:294798/1‑139 (MQ=255)
            gAGCTGTATCACTGATGATCTCCATTGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGAcaat                                                                                                                      <  1:294798/139‑1 (MQ=255)
               cTGTATCACTGATGATCTCCATTGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGAcaatcaa                                                                                                                   >  2:247392/1‑139 (MQ=255)
                tGTATCACTGATGATCTCCATTGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACt                                                                                                                            >  2:826964/1‑129 (MQ=255)
                tGTATCACTGATGATCTCCATTGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACt                                                                                                                            <  1:826964/129‑1 (MQ=255)
                  tATCACTGATGATCTCCATTGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACaatcaatc                                                                                                                 <  2:950755/138‑1 (MQ=255)
                  tATCACTGATGATCTCCATTGCCGGCGCGCAGGTACGTGCGCCCATCATTACGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATcc                                                                                                                >  2:980985/1‑139 (MQ=255)
                    tCACTGATGATCTCCATTGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCgg                                                                                                              <  1:114284/139‑1 (MQ=255)
                            gatCTCCATTGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCaa                                                                                                       >  1:1020670/1‑138 (MQ=255)
                            gatCTCCATTGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCaa                                                                                                       >  1:1140465/1‑138 (MQ=255)
                            gatCTCCATTGCCCGCGCGCAGGTAGGTGCGCCCCTCATCGCGCCGAACAACAGCGCGCGGTTCTTTCGCTATACATAAGCACCGAACAAGAACCAGATAATCACGGGTACCAGACTGACAATCAATCCGGCAATCaa                                                                                                       >  2:894488/1‑138 (MQ=255)
                               ctcCATTGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGCACCAGACTGACAATCAATCCGGCAATCAACATc                                                                                                   >  1:868822/1‑139 (MQ=255)
                                tcCATTGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGTACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCt                                                                                                  >  2:558910/1‑139 (MQ=255)
                                  cATTGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGa                                                                                                >  1:141987/1‑139 (MQ=255)
                                    ttGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGAcc                                                                                              >  2:978713/1‑139 (MQ=255)
                                    ttGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGAcc                                                                                              <  1:276302/139‑1 (MQ=255)
                                    ttGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGAcc                                                                                              >  2:843922/1‑139 (MQ=255)
                                       ccGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTAc                                                                                                                                  <  1:1058726/100‑1 (MQ=255)
                                       ccGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTAc                                                                                                                                  >  2:1058726/1‑100 (MQ=255)
                                          gcgcgcAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCCCCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAAt                                                                                        <  1:1100848/139‑1 (MQ=255)
                                          gcgcgcAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAAt                                                                                        <  1:441005/139‑1 (MQ=255)
                                           cgcgcAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATc                                                                                       <  1:682468/139‑1 (MQ=255)
                                                aGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGcc                                                                                  <  2:141987/139‑1 (MQ=255)
                                                     cGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATACCCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAgg                                                                              <  2:845797/138‑1 (MQ=255)
                                                     cGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAgg                                                                              <  1:247392/138‑1 (MQ=255)
                                                     cGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAgg                                                                              <  1:843922/138‑1 (MQ=255)
                                                     cGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAgg                                                                              <  1:205714/138‑1 (MQ=255)
                                                        gcgcCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGt                                                                          <  2:562483/139‑1 (MQ=255)
                                                              atcatTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCgg                                                                                                              >  2:768149/1‑97 (MQ=255)
                                                              atcatTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCgg                                                                                                              <  1:768149/97‑1 (MQ=255)
                                                                        cGAACAAAAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTAcc                                                          <  2:546821/139‑1 (MQ=255)
                                                                         gAACAACAGCGCGCGGTTCATTCGCAATACAGAAGCACCGAACAAGAAACATATAACCTCGAGTACCAGACTGCCAATCAATCCGCCAAGCAACTTCTGCCCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCg                                                         <  2:1151470/139‑1 (MQ=255)
                                                                          aacaacAGCGCGCGGTTCATTCGCGATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCgg                                                        >  1:712109/1‑139 (MQ=255)
                                                                              acaGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGCACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCt                                                    <  2:868822/139‑1 (MQ=255)
                                                                               caGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCt                                                    >  2:178875/1‑138 (MQ=255)
                                                                                   gcgcgGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGTTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGAc                                               <  2:676493/139‑1 (MQ=255)
                                                                                    cgcgGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGAc                                               >  1:91379/1‑138 (MQ=255)
                                                                                         ttCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGc                                         <  2:1177160/139‑1 (MQ=255)
                                                                                         ttCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGc                                         <  1:1204239/139‑1 (MQ=255)
                                                                                                      cATAAGCCCCAAACAAAAAACAGAAACCCACGGGTACAAGACTGACAATCAACCCGCCAATCAACATCTGCCCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAaca                             <  2:1000246/138‑1 (MQ=255)
                                                                                                          aGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCAt                         <  1:801245/138‑1 (MQ=255)
                                                                                                              ccGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGc                    >  2:367140/1‑139 (MQ=255)
                                                                                                              ccGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGc                    >  1:1000246/1‑139 (MQ=255)
                                                                                                              ccGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGc                    >  1:221811/1‑139 (MQ=255)
                                                                                                              ccGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGc                    >  1:386341/1‑139 (MQ=255)
                                                                                                                gAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGcc                   >  2:806791/1‑138 (MQ=255)
                                                                                                                gAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGcc                   >  2:862357/1‑138 (MQ=255)
                                                                                                                gAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCg                  <  1:1052439/139‑1 (MQ=255)
                                                                                                                 aaCAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATc                                                                                       <  2:390819/69‑1 (MQ=255)
                                                                                                                 aaCAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATc                                                                                       >  1:390819/1‑69 (MQ=255)
                                                                                                                 aaCAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGcc                                                                            <  1:904422/80‑1 (MQ=255)
                                                                                                                 aaCAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGcc                                                                            >  2:904422/1‑80 (MQ=255)
                                                                                                                 aaCAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCg                  >  1:1186173/1‑138 (MQ=255)
                                                                                                                 aaCAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCg                  >  2:75293/1‑138 (MQ=255)
                                                                                                                 aaCAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGCCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCg                  >  1:414061/1‑138 (MQ=255)
                                                                                                                  aCAAGAAACAGATAACCACGGTTACCAGACTGACAATCAATCCGGCAATCAACATCGGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGaa                <  1:862357/139‑1 (MQ=255)
                                                                                                                  aCAAGAAACAAACAACCCCGGGAACAGGCCGGACAAAAAACCCGCAAACCAACATCGGCCCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGaa                <  2:414061/139‑1 (MQ=255)
                                                                                                                            gATAACCACGGGTACCAGACTGAAAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGAACTCTTTCAcc      <  1:367140/139‑1 (MQ=255)
                                                                                                                             aTAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGc                                                                             >  1:971533/1‑67 (MQ=255)
                                                                                                                             aTAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGc                                                                             <  2:971533/67‑1 (MQ=255)
                                                                                                                             aTAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTAACGGCGCTCAGCCCAACGCCTGCCATAAACACCATCAAGCCGAACTCTTTCACCa     >  2:1146900/1‑139 (MQ=255)
                                                                                                                                aCCACGGTTCCCAGACTGACAATAAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGAACTCTTTCACCATGc  <  2:836962/139‑1 (MQ=255)
                                                                                                                                aCCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGAACTCTTTCACCATGc  >  1:836962/1‑139 (MQ=255)
                                                                                                                                    cGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCaa                                            <  2:893615/93‑1 (MQ=255)
                                                                                                                                    cGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCaa                                            >  1:893615/1‑93 (MQ=255)
                                                                                                                                        |                                                                                                                                  
GGATGTTACTGCGAGCTGTATCACTGATGATCTCCATTGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGCACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGAACTCTTTCACCATGC  >  NZ_CP009273/883685‑883951

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GCCCAGCGCCGGGATGTTACTGCGAGCTGTATCACTGATGATCTCCATTGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGCACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGAACTCTTTCACCATGCTTAATGCACCCTGCGGAAT  >  NZ_CP009273/883674‑883970
                                                                                                                                                   |                                                                                                                                                     
GCCCAGCGCCGGGATGTTACTGCGAGCTGTATCACTGATGATCTCCATTGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACC                                                                                                                                                    >  SRR3722129.983068/1‑151 (MQ=60)
  CCAGCGCCGGGATGTTACTGCGAGCTGTATCACTGATGATCTCCATTGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAG                                                                                                                                                  <  SRR3722129.645646/151‑1 (MQ=60)
    AGCGCCGGGATGTTACTGCGAGCTGTATCACTGATGATCTCCATTGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGAC                                                                                                                                                <  SRR3722129.571278/151‑1 (MQ=60)
     gctaatgacttatttttttattgttgtttgtttttttttttcaagcagaagacggcatacgatattcgccttagtctcgtgggctcggagatgtgtataagagacagACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGAC                                                                                                                                                <  SRR3722129.794879/43‑1 (MQ=60)
          gtgggctcggagatgtgtataagagacagGATCTCCATTGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAA                                                                                                                                          <  SRR3722129.1073603/122‑1 (MQ=60)
               GTTACTGCGAGCTGTATCACTGATGATCTCCATTGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAA                                                                                                                                      <  SRR3722129.1297367/150‑1 (MQ=60)
                 TACAGCGAGCTGTATCACTGATGATCTCCATTGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATC                                                                                                                                    <  SRR3722129.838849/150‑1 (MQ=60)
                      CGAGCTGTATCACTGATGATCTCCATTGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAA                                                                                                                              <  SRR3722129.299022/151‑1 (MQ=60)
                            GTATCACTGATGATCTCCATTGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAAC                                                                                                                         >  SRR3722129.1156329/1‑150 (MQ=60)
                            GTATCACTGATGATCTCCATTGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAAC                                                                                                                         >  SRR3722129.1035067/1‑150 (MQ=60)
                              ATCACTGATGATCTCCATTGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATC                                                                                                                      <  SRR3722129.115853/151‑1 (MQ=60)
                               TCACTGATGATCTCCATTGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGCACCAGACTGACAATCAATCCGGCAATCAACATCT                                                                                                                     >  SRR3722129.881226/1‑151 (MQ=60)
                                cgtgggctcggagatgtgtataagagacaGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATC                                                                                                                      <  SRR3722129.779317/120‑1 (MQ=60)
                                  CTGATGATCTCCATTGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGAC                                                                                                                  >  SRR3722129.143977/1‑151 (MQ=60)
                                              ATTGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGC                                                                                                      <  SRR3722129.280269/151‑1 (MQ=60)
                                                    GGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCCCCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGC                                                                                                <  SRR3722129.1116197/151‑1 (MQ=60)
                                                    GGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGC                                                                                                <  SRR3722129.447403/151‑1 (MQ=60)
                                                     GCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCC                                                                                               <  SRR3722129.692543/151‑1 (MQ=60)
                                                               ACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAAT                                                                                      <  SRR3722129.208586/150‑1 (MQ=60)
                                                               ACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAAT                                                                                      <  SRR3722129.856010/150‑1 (MQ=60)
                                                               ACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAAT                                                                                      <  SRR3722129.250898/150‑1 (MQ=60)
                                                                acgagattcgccttagtctcgtgggctcggagatgtgtataagagacagCATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATAC                                                                                    <  SRR3722129.917243/102‑1 (MQ=60)
                                                                          TCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCGATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGC                                                                          >  SRR3722129.722568/1‑151 (MQ=60)
                                                                          TCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACctgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtcttctgcttgaaa                                                                           >  SRR3722129.567004/1‑81 (MQ=60)
                                                                                    GAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACC                                                                 >  SRR3722129.92581/1‑150 (MQ=60)
                                                                                                   GTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAAC                                                 <  SRR3722129.1221041/151‑1 (MQ=60)
                                                                                                              ATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCC                                      >  SRR3722129.391988/1‑151 (MQ=60)
                                                                                                              ATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCC                                      >  SRR3722129.224897/1‑151 (MQ=60)
                                                                                                              ATACATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCC                                      >  SRR3722129.1014348/1‑151 (MQ=60)
                                                                                                                 CATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCctgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtcttc                                   >  SRR3722129.396518/1‑91 (MQ=60)
                                                                                                                 CATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGA                                    >  SRR3722129.1202679/1‑150 (MQ=60)
                                                                                                                 CATAAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGCCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGA                                    >  SRR3722129.420098/1‑150 (MQ=60)
                                                                                                                    AAGCACCGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGAACT                                 <  SRR3722129.812813/150‑1 (MQ=60)
                                                                                                                          CGAACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGAACTCTTTCAC                          <  SRR3722129.1067208/151‑1 (MQ=60)
                                                                                                                            AACAAGAAACAGATACCCACGGTTACCAGACTGACAATCAATCCGGCAATCAACATCGGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGAACTCTTTCACCA                        <  SRR3722129.874643/151‑1 (MQ=60)
                                                                                                                            AACAAGAAACAGATAACCACGGGTCCCAGACGGACAACCAATCCGGCAACCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGAACTCTTTCACCA                        <  SRR3722129.565483/151‑1 (MQ=60)
                                                                                                                             ACAAGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACtgtctcttatacacatctccgagcccacgagactaaggcgaatctcgtatgccgtcttctg                       >  SRR3722129.985282/1‑90 (MQ=60)
                                                                                                                                AGAAACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGAACTCTTTCACCATGCT                    >  SRR3722129.848974/1‑151 (MQ=60)
                                                                                                                                    ACAGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATActgtctcttatacacatctccgagcccacgagact                 >  SRR3722129.906312/1‑115 (MQ=60)
                                                                                                                                      AGATAACCACGGGTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGAACTCTTTCACCATGCTTAATGC              >  SRR3722129.200825/1‑151 (MQ=60)
                                                                                                                                      AGATAACCACGGGTACCAGACTGAAAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGAACTCTTTCACCATGCTTAATGC              <  SRR3722129.372521/151‑1 (MQ=60)
                                                                                                                                                  GTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGAACTCTTTCACCATGCTTAATGCACCCTGCGGAAT  >  SRR3722129.779658/1‑151 (MQ=60)
                                                                                                                                                  GTACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGAACTCTTTCACCATGCTTAATGCACCCTGCGGAAT  >  SRR3722129.422837/1‑151 (MQ=60)
                                                                                                                                                   |                                                                                                                                                     
GCCCAGCGCCGGGATGTTACTGCGAGCTGTATCACTGATGATCTCCATTGCCGGCGCGCAGGTACGTGCGCCCATCATTGCGCCGAACAACAGCGCGCGGTTCATTCGCAATACATAAGCACCGAACAAGAAACAGATAACCACGGGCACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGAACTCTTTCACCATGCTTAATGCACCCTGCGGAAT  >  NZ_CP009273/883674‑883970

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: