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 F0 I0 R1 81 326.4 10500568 93.1% 9776028 150.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 161,068 A→C 27.5% intergenic (+47/‑149) hrpB → / → mrcB putative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273161,0680AC27.5% 79.7 / 23.6 67intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (27/21);  new base C (2/17);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 7.88e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAA  >  CP009273/160920‑161213
                                                                                                                                                       |                                                                                                                                                 
cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CAt                                                                                                                                                  >  2:240676/1‑150 (MQ=255)
cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CAt                                                                                                                                                  >  2:240691/1‑150 (MQ=255)
  tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCa                                                                                                                                                >  2:2160509/1‑150 (MQ=255)
   gAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAg                                                                                                                                               >  1:4721589/1‑150 (MQ=255)
   gAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAg                                                                                                                                               >  1:4721485/1‑150 (MQ=255)
   gAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAg                                                                                                                                               >  1:4425266/1‑150 (MQ=255)
    aGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAgc                                                                                                                                              >  1:267020/1‑150 (MQ=255)
    aGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAgc                                                                                                                                              >  1:267059/1‑150 (MQ=255)
      gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCCgcgc                                                                                                                                            >  1:3902661/1‑150 (MQ=255)
      gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGGATTCGTAAGGGTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CCTGAgcgc                                                                                                                                            >  1:3902710/1‑150 (MQ=255)
        gCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAgcgcgc                                                                                                                                          >  1:2714881/1‑150 (MQ=255)
         cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCCGCGCGCa                                                                                                                                         >  1:345230/1‑150 (MQ=255)
          aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCCGCGCGCAc                                                                                                                                        >  1:6596/1‑150 (MQ=255)
           aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCACg                                                                                                                                       >  2:1906020/1‑150 (MQ=255)
           aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCACg                                                                                                                                       >  1:1616711/1‑150 (MQ=255)
            aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCCGCGCGCACGt                                                                                                                                      >  2:751756/1‑150 (MQ=255)
             aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCACGtt                                                                                                                                     >  1:5125952/1‑150 (MQ=255)
               gagaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCTACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCACTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCACggtgt                                                                                                                                   >  2:369694/1‑146 (MQ=255)
                ttaTGAATGGGGGTTATCCCAAACATGTCTGGCCGGACGACCGGGAAAAT‑ACTGCCCCGACGCGACGGACTAAAATGTTTTGGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGTGGG‑‑CATCAGCGCGCACGTtcac                                                                                                                                  <  2:3926215/148‑1 (MQ=25)
                agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCACGtgggg                                                                                                                                  >  1:1135123/1‑146 (MQ=255)
                 gaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CCTCAGCGCGCACGtggggg                                                                                                                                 >  1:29913/1‑145 (MQ=255)
                  aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCACgggggcgc                                                                                                                                >  1:689679/1‑143 (MQ=255)
                  aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCACGtgggccc                                                                                                                                >  1:689527/1‑144 (MQ=255)
                   tgactgggCTTTCTCCCAAACTTGTGTGGCGGGACTCCCCGGAAATT‑CCTGCACCGACGCGAGTGACTAAATTTTTTACTAATGGTTGGTTTACACCCACTCCCTGTGGTTTGGGGCGGGGTTGCGGG‑‑CATCAGCGCGCACGTtcaccct                                                                                                                               <  2:4367411/145‑1 (MQ=255)
                   tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCACGtggggggc                                                                                                                               >  1:4492161/1‑143 (MQ=255)
                   tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCACGtggcccca                                                                                                                               >  2:4784897/1‑143 (MQ=255)
                     ataGGGCGTTATCCCAAACTTGTCTGGCGGCACGACCCGTAAAAT‑CCGGCACCGACGCGACGGACTAAAATGTATTCGTATGGTTGTTTTTCTCCCCCACCCTGTGGGCGTGGGCGGGGGTGTGGG‑‑CATCAGCGCGCACGTtcaccctca                                                                                                                             <  2:3862594/148‑1 (MQ=255)
                       aGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCACGTgcacccccacc                                                                                                                           >  2:3247833/1‑150 (MQ=255)
                         gggtttACCCAAAACTTGTCTGGCCGCCCG‑CCCCGCAAATCAGTCCCCCGACCCGCCGCACTAAAAAGTTTTCGTATGGTTGTTTTACTCCCCCCCCGTGGGGGAGTGGGCGGGGGTGTGGC‑‑CATCAGCGCGCACGTTcaccctcaccct                                                                                                                         <  2:2043161/146‑1 (MQ=255)
                               aTCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCACGTgggcccgcacccacaccaa                                                                                                                   >  2:3718055/1‑143 (MQ=38)
                               aTCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCACGTggacccgcaccccagacca                                                                                                                   >  2:2039189/1‑143 (MQ=38)
                               aTCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCACGTGGCCCCGCACCCCAAaccc                                                                                                                   >  2:698298/1‑149 (MQ=38)
                                       cgtgggtGGCGGGCCGCCCGGGCAATT‑ACCGCCCCGCCGCGACGGCCGAAAATGTATGGGACTGGTGGTTTTTCCCCGCCGCCGGGGGGCTGTGGCCGGGGGTGTGGG‑‑CCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCa                                                                                                           <  2:3379490/144‑1 (MQ=255)
                                       cATGTGGGGCCGGACGCCCCGACAAAA‑ACTCCCCCCCCGCCCCGGACCAAATTTTTTTGTTATGTTGGGTTTTCCCCCACCCCGGGTGGGTGGGGGCGGGGGTGGGGC‑‑CATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCa                                                                                                           <  2:4172100/150‑1 (MQ=255)
                                        ttggggggCGGGACGACCCGCAAAAT‑CCTCCCCCGCCGCCGGGGCCGAAAATGTTTTCGTTGGGTTGTTTTTCGCCCCCGCCGGGTGGTCGTGGGGGGGGGGGGGGG‑‑CCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                          <  2:2408545/144‑1 (MQ=255)
                                          gtgggcccGGCCGCCCCGGAAAAT‑CCTGCCCCGCCGCGAGGTACTAAAATTTTTGCGAAAGGTGGGTTTTCCCCCGCGCCGGGTGGTTGTGGGCCGGGGTGGGGG‑‑CCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                        <  2:2830631/144‑1 (MQ=255)
                                                   gacgccCCGCCAAAT‑ACTGCCCCGCCGCGCCGGCCTAAAAAGTTTTCGTACGGTTGGTTTTCCCCCTCGCCCTGGGGGAGTGGGCCGGGGTGAGGG‑‑CATCAGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                               <  2:3307403/150‑1 (MQ=37)
                                                   gacgccCCGCCAAAT‑ACGGCCCCGACGCGCCGGACTAAAATGTATTCGTACGGTTGGTTTTCGCCCCCCCCCTGGGGGAGCGGGCCGGGGTGAGGC‑‑CATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                               <  2:4366834/150‑1 (MQ=255)
                                                      gccccGGCAAAT‑CCTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAAGGCCGGGGTGAGGG‑‑CATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCg                                                                                            <  1:4501169/148‑1 (MQ=255)
                                                      gacCCGGCAAAT‑ACTGCACCGACGCGCCGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAAGGCCGGGGTGAGGG‑‑CATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCACGGGAGAGGGGACCg                                                                                            <  1:4501227/150‑1 (MQ=255)
                                                               aaT‑ACGGCCCCGCCGCGCCGGACGAAAAATGTTTCGTAGGGTTGTTTTGCCCCCTCGCCCGGGGGGTGGGGGCGGGGGTGAGGG‑‑CCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGcac                                                                                   <  2:5148172/150‑1 (MQ=255)
                                                                    cggccccgccgCGCGGGACTAAAATGTTTTCGTATGGTTGGGTGTCCCCCCCCCCGTGTGGGAGGGGGCGGGGGTGCGGC‑‑CCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                                                               <  1:3423828/145‑1 (MQ=255)
                                                                     tGCCCCCACGCGACGGACTAAAAAGTATTCGTTCGGTTGGTTGTCTCCCTCCCCCTGGGGTTGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtt                                                                              <  2:3818401/150‑1 (MQ=255)
                                                                             cgcgccgGACGAAAAATTTTTCTTTGGGTTGGTTTTCGCCCACGCCGGGGGGTTGGGGGCGGGGGTGAGGG‑‑CATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagaga                                                                      <  2:2403733/150‑1 (MQ=255)
                                                                              gcgaggGACTAAAATTTATTCTTTTGTTTGTTTTTCACCCTCCCCGTGTGGTTGTGGGCGGGGGGGAGGG‑‑CATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                                                                     <  2:3210416/150‑1 (MQ=255)
                                                                                   ggactAAAATGTTTTCGTATGGTTGTTTTTCTCCCTCTCCGTGGGGGTGTGGGCGGGGGTGTGGG‑‑CATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtctt                                                                <  2:3374050/150‑1 (MQ=255)
                                                                                        aaattgTATTGTTATGGTGGGTTTTCTCCCTCGCCGTGGGGGAGTGGGCGGGGGTGTGGG‑‑CATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                           <  2:820846/145‑1 (MQ=255)
                                                                                         aattttATTGGTATGGTTTTTTTTCTCCTTATCCTGTTGGTTGGGGGCCGGGTTTTGGC‑‑CATCAGCGATCACGTTCACACTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTATTCTg                                                          <  1:2160509/145‑1 (MQ=255)
                                                                                                ttCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCACGTGCACCCTAACACAAAAACAAACCCACAACGGAGAGGGGTACGATCAAGCAAAAAATTTGAGATATATAATATTATTTCTTGtt                                                   >  2:1203682/1‑149 (MQ=12)
                                                                                                  cGTATGTTTGTTTTACCCCCTCCCCGGGTGGGTGTGGGCCGGGGTGCGGG‑‑CCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                 <  1:3251046/150‑1 (MQ=255)
                                                                                                               ttCCCCCACCCCGTGTGGCTGCGGGCGGGGGTGTGGG‑‑CATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaa                                    <  1:2173436/150‑1 (MQ=255)
                                                                                                               tactccctctccctGTGGTTGCGGGCGGGGGTGAGGG‑‑CCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaa                                    <  2:267059/148‑1 (MQ=255)
                                                                                                               taccccctccccttGTGGTTGTGGGCGGGGGTGCGGG‑‑CCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaa                                    <  1:2173555/146‑1 (MQ=255)
                                                                                                                tcgccccctccccGGGGGGGAGGGCGGGGGGGGGGGGCCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaa                                     <  1:3786592/147‑1 (MQ=255)
                                                                                                                 ccccctcccccgGGGGGTGTGGGCGGGGGGGTGGG‑‑CCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAAt                                  <  1:3489389/148‑1 (MQ=255)
                                                                                                                        tccgtGTGGTTGTGGGCGGGGGTGTGGG‑‑CATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtt                           <  2:2405599/150‑1 (MQ=255)
                                                                                                                         ccctGTGGTTGTGGGCGGGGGTGTGGG‑‑CATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                          <  2:4656722/150‑1 (MQ=255)
                                                                                                                                ggtgtgggCCGGGGTGAGGG‑‑CATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                   <  2:803118/145‑1 (MQ=255)
                                                                                                                                   gtgggCGGGGGTGAGGC‑‑CCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGaa                <  2:3621272/148‑1 (MQ=255)
                                                                                                                                   gggggCGGGGGTGCGGC‑‑CCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGaa                <  2:1616711/148‑1 (MQ=255)
                                                                                                                                   gcgggCCGGGGTGAGGG‑‑CATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGaa                <  2:2001330/148‑1 (MQ=255)
                                                                                                                                   gcgggCCGGGGGGTGGG‑‑CCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGaa                <  2:3619545/148‑1 (MQ=255)
                                                                                                                                   gagGGCCGGGGTGAGTG‑‑CATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGaa                <  2:2001180/150‑1 (MQ=255)
                                                                                                                                   gagGGCCGGGGTGAGCG‑‑CCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGaa                <  2:3619552/150‑1 (MQ=255)
                                                                                                                                          ggggTGCGGC‑‑CCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg         <  1:46482/150‑1 (MQ=255)
                                                                                                                                           gggggcggg‑‑CCTCAGCGCGCACGGTCACCCTCACCCTCACCCTCGCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTGCTTCTGTCTTGTAACAGAAGAACAGAAACTCGGGCTTTTGCGCCTGAATATTGCgg        <  2:3300276/144‑1 (MQ=255)
                                                                                                                                                 ggg‑‑CCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa  <  2:4492161/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                 
CGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAA  >  CP009273/160920‑161213

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

GATK/CNVnator alignment

N/A