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 F2 I207 R1 218 45.0 2482936 97.1% 2410930 86.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,624,386 T→C E522G (GAA→GGA)  yhiJ ← YhiJ/YhiL family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,624,3860TC98.8% 277.6 / ‑5.4 81E522G (GAA→GGA) yhiJYhiJ/YhiL family protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (50/30);  total (51/30)
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

GTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTTCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCA  >  NZ_CP009273/3624299‑3624472
                                                                                       |                                                                                       
gTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCa                                                                                       >  2:465490/1‑90 (MQ=255)
gTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCa                                                                                       >  1:604698/1‑90 (MQ=255)
  aaCATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATg                                                                                    >  1:1002063/1‑90 (MQ=255)
  aaCATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATg                                                                                    >  2:1053810/1‑90 (MQ=255)
   aCATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGc                                                                                   <  1:1039498/90‑1 (MQ=255)
        aaTATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTa                                                                              <  1:53846/90‑1 (MQ=255)
          tatTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAg                                                                            <  1:1098539/90‑1 (MQ=255)
            ttCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGAc                                                                          <  1:1056122/90‑1 (MQ=255)
             tCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTTCATTGCTGGGAAGAc                                                                          >  2:629774/1‑90 (MQ=255)
              cATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAg                                                                            >  2:2791/1‑86 (MQ=255)
              cATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAg                                                                            <  1:2791/86‑1 (MQ=255)
              cATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACtt                                                                        >  2:1048651/1‑90 (MQ=255)
              cATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACtt                                                                        >  2:464703/1‑90 (MQ=255)
                    gtCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAggg                                                                  >  2:1070010/1‑90 (MQ=255)
                      cTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAg                                                                >  2:578921/1‑90 (MQ=255)
                      cTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAg                                                                <  1:398397/90‑1 (MQ=255)
                          ggTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAGGGCATAAGTCGATTTCTGATTTTTCCATGCCGGTAAGGCTTCAAGGGAGTGGt                                                            >  1:77913/1‑90 (MQ=255)
                          ggTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGt                                                            <  2:756116/90‑1 (MQ=255)
                           gTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAGGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGtt                                                           <  2:771561/90‑1 (MQ=255)
                           gTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGtt                                                           <  1:739347/90‑1 (MQ=255)
                            ttCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTa                                                          <  1:983920/90‑1 (MQ=255)
                             tCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTAt                                                         >  2:54087/1‑90 (MQ=255)
                                gagCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATg                                                                                    >  1:785446/1‑60 (MQ=255)
                                gagCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATg                                                                                    <  2:785446/60‑1 (MQ=255)
                                gagCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTAtctc                                                      >  2:546586/1‑90 (MQ=255)
                                gagCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTAtctc                                                      >  1:69244/1‑90 (MQ=255)
                                gagCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTAtctc                                                      >  2:102058/1‑90 (MQ=255)
                                gagCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTAtctc                                                      >  1:360811/1‑90 (MQ=255)
                                gagCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTAtctc                                                      >  1:454926/1‑90 (MQ=255)
                                    tGTTACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTTTTTTTTCCAATCTGGTAAGGCTTCCAGGGAGTGGTTATCTCAGGc                                                  <  1:464703/90‑1 (MQ=255)
                                    tGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGc                                                  <  1:1238978/90‑1 (MQ=255)
                                     gTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAg                                                                >  1:670348/1‑75 (MQ=255)
                                     gTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAg                                                                <  2:670348/75‑1 (MQ=255)
                                     gTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGtt                                                           >  2:89310/1‑80 (MQ=255)
                                     gTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGt                                                            <  1:410453/79‑1 (MQ=255)
                                     gTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGt                                                            >  2:410453/1‑79 (MQ=255)
                                     gTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTa                                                          <  1:1074068/81‑1 (MQ=255)
                                     gTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTa                                                          >  2:1074068/1‑81 (MQ=255)
                                     gTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGcc                                                 >  2:853433/1‑90 (MQ=255)
                                     gTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGcc                                                 >  1:544094/1‑90 (MQ=255)
                                     gTAACATTGTGCCCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCCGGTAAGACTTCCAGGGAGTGGtt                                                           <  1:89310/80‑1 (MQ=255)
                                      tAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGccc                                                <  2:1209206/90‑1 (MQ=255)
                                          aTTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCa                                                                                       <  2:465476/48‑1 (MQ=255)
                                          aTTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCa                                            >  2:905314/1‑90 (MQ=255)
                                          aTTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCa                                            >  1:1069744/1‑90 (MQ=255)
                                           ttGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCAc                                           >  2:518078/1‑90 (MQ=255)
                                           ttGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCAc                                           <  1:1037350/90‑1 (MQ=255)
                                           ttGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCAc                                           <  2:174217/90‑1 (MQ=255)
                                           ttGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCAc                                           >  1:635019/1‑90 (MQ=255)
                                            tgtgCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAggg                                                                  <  2:1219907/66‑1 (MQ=255)
                                            tgtgCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAggg                                                                  >  1:1219907/1‑66 (MQ=255)
                                            tgtgCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCAcc                                          >  2:149389/1‑90 (MQ=255)
                                            tgtgCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCAcc                                          <  1:927842/90‑1 (MQ=255)
                                            tgtgCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCAcc                                          >  2:897516/1‑90 (MQ=255)
                                             gtgCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCa                                         >  2:453340/1‑90 (MQ=255)
                                             gtgCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCa                                         >  2:912420/1‑90 (MQ=255)
                                             gtgCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCa                                         >  2:78407/1‑90 (MQ=255)
                                             gtgCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCa                                         >  2:507050/1‑90 (MQ=255)
                                             gtgCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCa                                         >  1:561633/1‑90 (MQ=255)
                                               gcgcCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAggg                                                                  >  1:619751/1‑63 (MQ=255)
                                               gcgcCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAggg                                                                  <  2:619751/63‑1 (MQ=255)
                                               gcgcCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAgg                                       >  2:204140/1‑90 (MQ=255)
                                               gcgcCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAgg                                       >  1:347835/1‑90 (MQ=255)
                                               gcgcCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAgg                                       >  2:624356/1‑90 (MQ=255)
                                               gcgcCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAgg                                       >  2:525763/1‑90 (MQ=255)
                                                    ttCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTaa                                  >  2:1134850/1‑90 (MQ=255)
                                                             aaGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAg                                                                <  2:978218/51‑1 (MQ=255)
                                                             aaGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAg                                                                >  1:978218/1‑51 (MQ=255)
                                                                gCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTat                      >  2:620050/1‑90 (MQ=255)
                                                                 cATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTatg                     <  1:373048/90‑1 (MQ=255)
                                                                  aTAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGc                                             <  1:499416/65‑1 (MQ=255)
                                                                  aTAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGc                                             >  2:499416/1‑65 (MQ=255)
                                                                  aTAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTatgc                    >  2:768111/1‑90 (MQ=255)
                                                                  aTAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTatgc                    >  1:212333/1‑90 (MQ=255)
                                                                         cATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGa             <  1:26525/90‑1 (MQ=255)
                                                                          aTTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTatg                     <  2:1219038/81‑1 (MQ=255)
                                                                          aTTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTatg                     >  1:1219038/1‑81 (MQ=255)
                                                                          aTTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGaa            >  1:666472/1‑90 (MQ=255)
                                                                          aTTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGaa            >  2:222901/1‑90 (MQ=255)
                                                                                gTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCAcc                                          <  1:802983/54‑1 (MQ=255)
                                                                                gTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCAcc                                          >  2:802983/1‑54 (MQ=255)
                                                                                 ttttttcCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACAttctt     <  1:507050/90‑1 (MQ=255)
                                                                                    tttcCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCa  <  2:778047/90‑1 (MQ=255)
                                                                                    tttcCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCa  <  1:825264/90‑1 (MQ=255)
                                                                                    tttcCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCa  <  2:604369/90‑1 (MQ=255)
                                                                                       |                                                                                       
GTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTTCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCA  >  NZ_CP009273/3624299‑3624472

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTTCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGT  >  NZ_CP009273/3624289‑3624485
                                                                                                 |                                                                                                   
CAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCA                                                                                                   >  SRR3722089.612172/1‑100 (MQ=60)
  ACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATG                                                                                                 >  SRR3722089.1016963/1‑100 (MQ=60)
        TAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTA                                                                                           <  SRR3722089.15608/100‑1 (MQ=60)
             ACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACT                                                                                      <  SRR3722089.1055056/100‑1 (MQ=60)
                  AATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAG                                                                                 <  SRR3722089.54515/100‑1 (MQ=60)
                    TATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGG                                                                               <  SRR3722089.1115133/100‑1 (MQ=60)
                    TATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGG                                                                               <  SRR3722089.2816/100‑1 (MQ=60)
                      TTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAG                                                                             <  SRR3722089.1071978/100‑1 (MQ=60)
                          TGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAGGGCATAAGTCGATTTCTGATTTTTCCATGCCGGTAAGGCTTCAAGGGAGTGGT                                                                         >  SRR3722089.78866/1‑100 (MQ=60)
                                CTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACctgtctcttatacacatctg                                                                   >  SRR3722089.796194/1‑80 (MQ=60)
                                CTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTC                                                                   >  SRR3722089.460372/1‑100 (MQ=60)
                                CTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTC                                                                   <  SRR3722089.403240/100‑1 (MQ=60)
                                CTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTC                                                                   >  SRR3722089.365260/1‑100 (MQ=60)
                                CTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTC                                                                   >  SRR3722089.70080/1‑100 (MQ=60)
                                    GGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGC                                                               <  SRR3722089.415419/100‑1 (MQ=60)
                                     GTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCctgtc                                                              >  SRR3722089.679025/1‑95 (MQ=60)
                                     GTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCC                                                              >  SRR3722089.550428/1‑100 (MQ=60)
                                     GTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCC                                                              <  SRR3722089.749026/100‑1 (MQ=60)
                                     GTTCAGAGCTGTAACATTGTGCCCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCCGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCC                                                              <  SRR3722089.90412/100‑1 (MQ=60)
                                      TTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCC                                                             <  SRR3722089.1090211/100‑1 (MQ=60)
                                      TTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCC                                                             <  SRR3722089.998507/100‑1 (MQ=60)
                                          GAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGctgtctcttatacacatctgacgctgccgacg                                                         >  SRR3722089.471039/1‑68 (MQ=60)
                                          GAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCA                                                         >  SRR3722089.1085808/1‑100 (MQ=60)
                                           AGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCAC                                                        >  SRR3722089.643103/1‑100 (MQ=60)
                                            GCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCctgtctcttataca                                                       >  SRR3722089.1238464/1‑86 (MQ=60)
                                             CTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCA                                                      >  SRR3722089.568198/1‑100 (MQ=60)
                                              TGTTACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTTTTTTTTCCAATCTGGTAAGGCTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAG                                                     <  SRR3722089.470253/100‑1 (MQ=60)
                                              TGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAG                                                     <  SRR3722089.1257829/100‑1 (MQ=60)
                                               GTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCctgtctcttatacacat                                                    >  SRR3722089.627533/1‑83 (MQ=60)
                                               GTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGG                                                    >  SRR3722089.352158/1‑100 (MQ=60)
                                                   gtgtataagagacagAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTA                                                <  SRR3722089.505324/85‑1 (MQ=60)
                                                     TTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAAT                                              <  SRR3722089.1052870/100‑1 (MQ=60)
                                                      cagcgtcagatgtgtataagagacagGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATT                                             <  SRR3722089.814129/74‑1 (MQ=60)
                                                      TGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATT                                             <  SRR3722089.941380/100‑1 (MQ=60)
                                                             CTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCctgtctcttatacacatctgacgctgccg                                      >  SRR3722089.992688/1‑71 (MQ=60)
                                                                  AATCAAAGGCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGC                                 >  SRR3722089.215011/1‑100 (MQ=60)
                                                                          GCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAA                         >  SRR3722089.675082/1‑100 (MQ=60)
                                                                          GCATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAA                         >  SRR3722089.1237582/1‑100 (MQ=60)
                                                                           CATAAGTCCATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAAC                        <  SRR3722089.377622/100‑1 (MQ=60)
                                                                                   CATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGtctcttatacacatctgacgctgccgacgaccatcttagtgta                >  SRR3722089.871165/1‑57 (MQ=60)
                                                                                   CATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCC                <  SRR3722089.26846/100‑1 (MQ=60)
                                                                                    ATTTCTGTTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCA               >  SRR3722089.1137982/1‑100 (MQ=60)
                                                                                           TTTTTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTT        <  SRR3722089.513030/100‑1 (MQ=60)
                                                                                              TTTCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTAC     <  SRR3722089.836895/100‑1 (MQ=60)
                                                                                                TCCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTG   <  SRR3722089.627848/100‑1 (MQ=60)
                                                                                                 CCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGT  <  SRR3722089.363706/100‑1 (MQ=60)
                                                                                                 |                                                                                                   
CAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATAAGTCCATTTCTGTTTTTTTCATGCTGGTAAGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGT  >  NZ_CP009273/3624289‑3624485

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