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 F1 I2 R1 1873 183.1 9376148 98.2% 9207377 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,987,935 G→A R156C (CGC→TGC)  hemC ← porphobilinogen deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,987,9350GA99.6% 825.0 / ‑5.7 230R156C (CGC→TGC) hemCporphobilinogen deaminase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (116/113);  total (117/113)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 247 total aligned reads displayed.

GACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCGCAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCTGCCGGT  >  NC_000913/3987845‑3988023
                                                                                          |                                                                                        
gACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGcaca                                                                                        >  1:1317882/1‑93 (MQ=255)
gACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGcaca                                                                                        >  2:3053119/1‑93 (MQ=255)
gACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGcaca                                                                                        >  1:94322/1‑93 (MQ=255)
 aCCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAg                                                                                       <  1:994715/93‑1 (MQ=255)
  ccTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAgg                                                                                      >  2:3295561/1‑93 (MQ=255)
  ccTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAgg                                                                                      >  1:1490565/1‑93 (MQ=255)
   cTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAggg                                                                                     >  2:3171635/1‑93 (MQ=255)
   cTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAggg                                                                                     <  1:4111444/93‑1 (MQ=255)
      aaCGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGc                                                                                  >  1:1574036/1‑93 (MQ=255)
      aaCGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGc                                                                                  <  1:397714/93‑1 (MQ=255)
         gTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGa                                                                               >  2:206132/1‑93 (MQ=255)
          ttttAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGAt                                                                              >  1:2644512/1‑93 (MQ=255)
          ttttAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGAt                                                                              >  1:649428/1‑93 (MQ=255)
          ttttAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGAt                                                                              >  1:1592189/1‑93 (MQ=255)
            ttAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATaa                                                                            <  2:2191909/93‑1 (MQ=255)
            ttAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATaa                                                                            >  1:2821332/1‑93 (MQ=255)
            ttAGTCCGGCTACGGAAAGAATGATGGCTCCGTATTCTCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTTCCGCACAGGGAGCGGATaa                                                                            <  2:1265361/93‑1 (MQ=255)
               gTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCa                                                                         <  2:1315385/93‑1 (MQ=255)
               gTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCa                                                                         >  2:3124675/1‑93 (MQ=255)
                 ccGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGa                                                                       <  1:1846716/93‑1 (MQ=255)
                  cGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGAt                                                                      >  2:600709/1‑93 (MQ=255)
                    gCTACAGCACGAATGAAGGCATCGTATGAGGCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGACGCCGCACAGGGAGCGGATAATCAGAttt                                                                    >  2:1520101/1‑91 (MQ=255)
                        cgtcAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGAc                                                                <  1:3907350/90‑1 (MQ=255)
                           cAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCAGGACGGc                                                             >  2:4144307/1‑93 (MQ=255)
                            aaGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGTCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCg                                                            >  1:781066/1‑93 (MQ=255)
                             aGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGt                                                           <  1:1896487/93‑1 (MQ=255)
                               aaTGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTc                                                         <  2:2534848/93‑1 (MQ=255)
                                atgatgGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCa                                                        >  2:48744/1‑93 (MQ=255)
                                atgatgGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCa                                                        >  1:2788660/1‑93 (MQ=255)
                                atgatgGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGCAGCGGATAATCAGATCCGGACGGCGTTCa                                                        >  1:2737795/1‑93 (MQ=255)
                                 tgatgGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAg                                                       >  1:753396/1‑93 (MQ=255)
                                  gatgGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGc                                                      <  2:2286245/93‑1 (MQ=255)
                                   atgGCATCTTATTCTCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGcc                                                     <  1:4441831/93‑1 (MQ=255)
                                    tgGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCa                                                    <  2:701974/93‑1 (MQ=255)
                                    tgGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCa                                                    <  1:102369/93‑1 (MQ=255)
                                    tgGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCa                                                    <  2:1877212/93‑1 (MQ=255)
                                     ggCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAg                                                   <  1:520199/93‑1 (MQ=255)
                                     ggCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAg                                                   >  1:2325921/1‑93 (MQ=255)
                                      gCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGt                                                  <  1:2278795/93‑1 (MQ=255)
                                      gCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGt                                                  <  1:3657324/93‑1 (MQ=255)
                                       cATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGtt                                                 >  1:827259/1‑93 (MQ=255)
                                        aTCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTg                                                <  2:2647546/93‑1 (MQ=255)
                                        aTCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTg                                                >  1:3042205/1‑93 (MQ=255)
                                        aTCGTATTAGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTg                                                >  2:4554773/1‑93 (MQ=255)
                                           gTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCa                                             >  2:1164176/1‑93 (MQ=255)
                                            tATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCAc                                            >  1:4147365/1‑93 (MQ=255)
                                               tCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTgg                                         <  1:401209/93‑1 (MQ=255)
                                                cGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGc                                        <  2:2590404/93‑1 (MQ=255)
                                                 gCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCg                                       >  2:3801074/1‑93 (MQ=255)
                                                 gCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCg                                       >  2:2452365/1‑93 (MQ=255)
                                                  ccGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGa                                      <  1:905858/93‑1 (MQ=255)
                                                  ccGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGa                                      >  1:1467920/1‑93 (MQ=255)
                                                  ccGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGa                                      <  2:1587865/93‑1 (MQ=255)
                                                     ttATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGt                                   >  2:428090/1‑93 (MQ=255)
                                                     ttATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGt                                   <  2:229037/93‑1 (MQ=255)
                                                      tATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTa                                  >  2:4420663/1‑93 (MQ=255)
                                                       aTCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTaa                                 >  1:2182773/1‑93 (MQ=255)
                                                        tCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTaaa                                >  1:199122/1‑93 (MQ=255)
                                                        tCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTaaa                                >  1:4210807/1‑93 (MQ=255)
                                                        tCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTaaa                                <  2:3421832/93‑1 (MQ=255)
                                                        tCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGGTCAGCCAGTTGGTACGGGCGACGTgaa                                >  2:1366482/1‑93 (MQ=255)
                                                         ccAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAAc                               <  1:583295/93‑1 (MQ=255)
                                                         ccAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAAc                               <  2:79164/93‑1 (MQ=255)
                                                         ccAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAAc                               <  1:3418672/93‑1 (MQ=255)
                                                          cAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACt                              <  1:2599289/93‑1 (MQ=255)
                                                          cAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACt                              <  2:535816/93‑1 (MQ=255)
                                                          cAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACt                              <  1:391932/93‑1 (MQ=255)
                                                            gTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTgg                            <  2:2334092/93‑1 (MQ=255)
                                                             tttGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGa                           <  2:2925281/93‑1 (MQ=255)
                                                             tttGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGa                           <  2:1665371/93‑1 (MQ=255)
                                                                gCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGt                        >  2:3702850/1‑93 (MQ=255)
                                                                 cTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTc                       <  2:4369936/93‑1 (MQ=255)
                                                                  tCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTcc                      >  1:1664174/1‑93 (MQ=255)
                                                                  tCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTcc                      <  1:433653/93‑1 (MQ=255)
                                                                   cAGCCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTccc                     <  2:868050/93‑1 (MQ=255)
                                                                     ggCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCcga                   >  1:3178690/1‑93 (MQ=255)
                                                                     ggCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTCCAGCCAGTTGGCACTGGCGACGTAAACTGGCGCTCCcga                   >  1:1188833/1‑93 (MQ=255)
                                                                         aGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCCGCCAGTTGGCACTGGCGCAGTAAACTGGACGTCCCGAAGAt               >  2:1226187/1‑93 (MQ=255)
                                                                          gTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATa              <  1:4603985/93‑1 (MQ=255)
                                                                          gTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATa              <  2:2079195/93‑1 (MQ=255)
                                                                           tGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATAc             <  2:584302/93‑1 (MQ=255)
                                                                           tGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATAc             <  2:3566438/93‑1 (MQ=255)
                                                                           tGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATAc             <  1:4494088/93‑1 (MQ=255)
                                                                           tGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATAc             <  2:1846310/93‑1 (MQ=255)
                                                                            gCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATAct            >  1:3725963/1‑93 (MQ=255)
                                                                             ccGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATActg           >  1:1105186/1‑93 (MQ=255)
                                                                              cGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATActgc          <  2:3522714/93‑1 (MQ=255)
                                                                               gACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATActgcc         <  2:3615735/93‑1 (MQ=255)
                                                                                aCGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATActgcct        >  2:3824764/1‑93 (MQ=255)
                                                                                aCGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATActgcct        <  2:2172306/93‑1 (MQ=255)
                                                                                aCGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATActgcct        >  2:3764094/1‑93 (MQ=255)
                                                                                aCGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATActgcct        <  1:1367958/93‑1 (MQ=255)
                                                                                 cGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATActgcctg       >  2:3610686/1‑93 (MQ=255)
                                                                                  gTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATActgcctgc      >  1:3292637/1‑93 (MQ=255)
                                                                                  gTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATActgcctgc      >  1:1825711/1‑93 (MQ=255)
                                                                                  gTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATActgcctgc      <  2:742307/93‑1 (MQ=255)
                                                                                   ttGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACtgcctgcc     <  2:1967635/93‑1 (MQ=255)
                                                                                   ttGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACtgcctgcc     >  2:1030644/1‑93 (MQ=255)
                                                                                     gCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCTGCCgg   >  2:2849335/1‑93 (MQ=255)
                                                                                     gCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCTGCCgg   >  2:2898493/1‑93 (MQ=255)
                                                                                          |                                                                                        
GACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCGCAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCTGCCGGT  >  NC_000913/3987845‑3988023

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCGCAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCTGCCGGTAACGC  >  NC_000913/3987842‑3988028
                                                                                             |                                                                                             
CCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATTATGGCATCGTATTCGCCGTTATCCAGTCTGCTCAGGCGAGTGCCGAGGTTGCCGCA                                                                                               >  SRR4302535.863367/1‑94 (MQ=60)
CCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCA                                                                                               >  SRR4302535.3322311/1‑94 (MQ=60)
  AGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACA                                                                                             >  SRR4302535.94995/1‑94 (MQ=60)
  AGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACA                                                                                             <  SRR4302535.3223839/94‑1 (MQ=60)
  AGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACA                                                                                             <  SRR4302535.2584476/94‑1 (MQ=60)
  AGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACA                                                                                             >  SRR4302535.4218390/1‑94 (MQ=60)
  AGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACA                                                                                             >  SRR4302535.1327296/1‑94 (MQ=60)
   GACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAG                                                                                            <  SRR4302535.64709/94‑1 (MQ=60)
    ACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCTCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGG                                                                                           <  SRR4302535.1740020/94‑1 (MQ=60)
    ACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGG                                                                                           >  SRR4302535.1501268/1‑94 (MQ=60)
    ACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGG                                                                                           <  SRR4302535.1001779/94‑1 (MQ=60)
    ACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGG                                                                                           >  SRR4302535.2663096/1‑94 (MQ=60)
    ACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGG                                                                                           >  SRR4302535.683736/1‑94 (MQ=60)
     CCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGG                                                                                          >  SRR4302535.2739379/1‑94 (MQ=60)
     CCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGG                                                                                          <  SRR4302535.424918/94‑1 (MQ=60)
     CCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGG                                                                                          <  SRR4302535.927771/94‑1 (MQ=60)
      CTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCNTCGTATTCGCCGTNATCCAGTTTGCNCAGGCGAGTGNCGACGNTNCCNCNCNGNGA                                                                                         >  SRR4302535.382509/1‑94 (MQ=60)
      CTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGA                                                                                         <  SRR4302535.4139586/94‑1 (MQ=60)
      CTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGA                                                                                         >  SRR4302535.1005124/1‑94 (MQ=60)
       TAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAG                                                                                        >  SRR4302535.1495436/1‑94 (MQ=60)
        AAACGTTTTAGTCCGGCTACGGTAAGAATGATGGCATCTTATTCGACGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGC                                                                                       <  SRR4302535.2772448/94‑1 (MQ=60)
        AAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGG                                                                                       >  SRR4302535.3980031/1‑94 (MQ=60)
        AAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGC                                                                                       >  SRR4302535.3797072/1‑94 (MQ=60)
        AAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGC                                                                                       >  SRR4302535.1585400/1‑94 (MQ=60)
        AAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGC                                                                                       >  SRR4302535.643741/1‑94 (MQ=60)
        AAACGTTTTAGTCCGGCGACGTCAAGAATGATGGCATCGTATTCGTCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGC                                                                                       <  SRR4302535.3127623/94‑1 (MQ=60)
         AACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCG                                                                                      <  SRR4302535.400501/94‑1 (MQ=60)
           CGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGA                                                                                    >  SRR4302535.2477820/1‑94 (MQ=60)
           CGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGA                                                                                    >  SRR4302535.1989229/1‑94 (MQ=60)
            GTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGAT                                                                                   >  SRR4302535.1603677/1‑94 (MQ=60)
            GTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGAT                                                                                   >  SRR4302535.2663452/1‑94 (MQ=60)
            GTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGAT                                                                                   >  SRR4302535.653956/1‑94 (MQ=60)
            GTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGAT                                                                                   >  SRR4302535.2313267/1‑94 (MQ=60)
            GTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGAT                                                                                   >  SRR4302535.2072508/1‑94 (MQ=60)
             TTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATA                                                                                  >  SRR4302535.1586336/1‑94 (MQ=60)
              TTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAA                                                                                 >  SRR4302535.2841368/1‑94 (MQ=60)
                  GTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAG                                                                             >  SRR4302535.905863/1‑94 (MQ=60)
                    CCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGAT                                                                           <  SRR4302535.1860025/94‑1 (MQ=60)
                      GGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCC                                                                         >  SRR4302535.1843118/1‑94 (MQ=60)
                      GGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCC                                                                         >  SRR4302535.774073/1‑94 (MQ=60)
                          ACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGAC                                                                     <  SRR4302535.1521816/94‑1 (MQ=60)
                           CGTCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACG                                                                    <  SRR4302535.3934119/94‑1 (MQ=60)
                           CGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACG                                                                    <  SRR4302535.816191/94‑1 (MQ=60)
                            GGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGG                                                                   >  SRR4302535.3282213/1‑94 (MQ=60)
                             GTAAGAATGATGGCATCGTATTCTCCGTTATCCAGTTTGCTCAGTCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGC                                                                  <  SRR4302535.880681/94‑1 (MQ=60)
                              CAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGTCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCG                                                                 >  SRR4302535.786549/1‑94 (MQ=60)
                               AAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGT                                                                >  SRR4302535.1320971/1‑94 (MQ=60)
                                AGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTT                                                               <  SRR4302535.1910132/94‑1 (MQ=60)
                                AGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTT                                                               >  SRR4302535.4255602/1‑94 (MQ=60)
                                  AATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCA                                                             <  SRR4302535.4667325/94‑1 (MQ=60)
                                  AATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCA                                                             >  SRR4302535.2808505/1‑94 (MQ=60)
                                  AATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGCAGCGGATAATCAGATCCGGACGGCGTTCA                                                             >  SRR4302535.2757323/1‑94 (MQ=60)
                                   ATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAG                                                            >  SRR4302535.758700/1‑94 (MQ=60)
                                   ATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAG                                                            >  SRR4302535.1901115/1‑94 (MQ=60)
                                     GATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCC                                                          <  SRR4302535.1975911/94‑1 (MQ=60)
                                     GATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCC                                                          >  SRR4302535.2812776/1‑94 (MQ=60)
                                      ATGGCATCTTATTCTCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCA                                                         <  SRR4302535.4472289/94‑1 (MQ=60)
                                      ATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCA                                                         >  SRR4302535.50946/1‑94 (MQ=60)
                                       TGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAG                                                        >  SRR4302535.2342563/1‑94 (MQ=60)
                                       TGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAG                                                        <  SRR4302535.103091/94‑1 (MQ=60)
                                        GGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGT                                                       <  SRR4302535.523836/94‑1 (MQ=60)
                                         GCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTT                                                      <  SRR4302535.2295071/94‑1 (MQ=60)
                                         GCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTT                                                      >  SRR4302535.833050/1‑94 (MQ=60)
                                         GCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTT                                                      <  SRR4302535.3256298/94‑1 (MQ=60)
                                         GCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTT                                                      <  SRR4302535.3682640/94‑1 (MQ=60)
                                          CATCGTATTCTCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTG                                                     <  SRR4302535.1565328/94‑1 (MQ=60)
                                          CATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTG                                                     >  SRR4302535.3063521/1‑94 (MQ=60)
                                           ATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCGCAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGG                                                    >  SRR4302535.346722/1‑94 (MQ=60)
                                           ATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGG                                                    <  SRR4302535.2959586/94‑1 (MQ=60)
                                              GTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCAC                                                 >  SRR4302535.98866/1‑94 (MQ=60)
                                              GTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCAC                                                 >  SRR4302535.3916235/1‑94 (MQ=60)
                                              GTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCAC                                                 >  SRR4302535.4175764/1‑94 (MQ=60)
                                                  TCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGC                                             <  SRR4302535.404022/94‑1 (MQ=60)
                                                   CGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCG                                            >  SRR4302535.3735768/1‑94 (MQ=60)
                                                    GCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGA                                           >  SRR4302535.1478435/1‑94 (MQ=60)
                                                     CCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGTATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGAC                                          <  SRR4302535.3323183/94‑1 (MQ=60)
                                                     CCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGAC                                          <  SRR4302535.912285/94‑1 (MQ=60)
                                                       GTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGT                                        >  SRR4302535.4373812/1‑94 (MQ=60)
                                                       GTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGT                                        >  SRR4302535.4164725/1‑94 (MQ=60)
                                                       GTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGT                                        <  SRR4302535.402157/94‑1 (MQ=60)
                                                        TTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTA                                       >  SRR4302535.1098884/1‑94 (MQ=60)
                                                        TTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTA                                       <  SRR4302535.4487951/94‑1 (MQ=60)
                                                         TATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAA                                      >  SRR4302535.2198336/1‑94 (MQ=60)
                                                          ATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAA                                     >  SRR4302535.4239626/1‑94 (MQ=60)
                                                          ATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAA                                     >  SRR4302535.200486/1‑94 (MQ=60)
                                                           TCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAAC                                    >  SRR4302535.447703/1‑94 (MQ=60)
                                                           TCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAAC                                    <  SRR4302535.1673134/94‑1 (MQ=60)
                                                           TCCACTTTGCTCAGGCGAGCTCCGACGTTGCCCCACAGAGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAAC                                    <  SRR4302535.4082691/94‑1 (MQ=60)
                                                            CCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACT                                   >  SRR4302535.3483879/1‑94 (MQ=60)
                                                            CCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACT                                   <  SRR4302535.3442452/94‑1 (MQ=60)
                                                            CCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACT                                   <  SRR4302535.587386/94‑1 (MQ=60)
                                                            CCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGACAGTTGGCACTGGCGACGTAAACT                                   >  SRR4302535.781917/1‑94 (MQ=60)
                                                             CAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTG                                  <  SRR4302535.394666/94‑1 (MQ=60)
                                                             CAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTG                                  <  SRR4302535.2617950/94‑1 (MQ=60)
                                                             CAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTG                                  >  SRR4302535.2623038/1‑94 (MQ=60)
                                                              AGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGG                                 >  SRR4302535.2207439/1‑94 (MQ=60)
                                                                TTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGAC                               <  SRR4302535.1477941/94‑1 (MQ=60)
                                                                 TTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGTACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACG                              <  SRR4302535.4066628/94‑1 (MQ=60)
                                                                 TTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGCCGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACG                              <  SRR4302535.1140291/94‑1 (MQ=60)
                                                                    CTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCC                           >  SRR4302535.1676200/1‑94 (MQ=60)
                                                                    CTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCC                           <  SRR4302535.530807/94‑1 (MQ=60)
                                                                     TCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCC                          <  SRR4302535.436673/94‑1 (MQ=60)
                                                                      CAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCG                         >  SRR4302535.2917762/1‑94 (MQ=60)
                                                                       AGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGA                        >  SRR4302535.1858419/1‑94 (MQ=60)
                                                                       AGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGA                        >  SRR4302535.3200672/1‑94 (MQ=60)
                                                                       AGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTCCAGCCAGTTGGCACTGGCGACGTAAACTGGcgctcccga                        >  SRR4302535.1197347/1‑85 (MQ=60)
                                                                           GAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGAT                    <  SRR4302535.3984584/94‑1 (MQ=60)
                                                                            AGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATA                   <  SRR4302535.4527899/94‑1 (MQ=60)
                                                                            AGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATA                   >  SRR4302535.3934170/1‑94 (MQ=60)
                                                                             GTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATAC                  <  SRR4302535.4635286/94‑1 (MQ=60)
                                                                              TGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACT                 <  SRR4302535.4524936/94‑1 (MQ=60)
                                                                              TGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACT                 >  SRR4302535.3751777/1‑94 (MQ=60)
                                                                               GCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTG                >  SRR4302535.1113083/1‑94 (MQ=60)
                                                                                CCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGC               >  SRR4302535.98652/1‑94 (MQ=60)
                                                                                CCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGC               <  SRR4302535.2769985/94‑1 (MQ=60)
                                                                                 CGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCC              <  SRR4302535.1650092/94‑1 (MQ=60)
                                                                                  GACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCT             <  SRR4302535.1295143/94‑1 (MQ=60)
                                                                                   ACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCTG            <  SRR4302535.1377733/94‑1 (MQ=60)
                                                                                   ACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCTG            <  SRR4302535.1340928/94‑1 (MQ=60)
                                                                                    CGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCTGC           >  SRR4302535.3315530/1‑94 (MQ=60)
                                                                                    CGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCTGC           >  SRR4302535.1838892/1‑94 (MQ=60)
                                                                                     GTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCTGCC          >  SRR4302535.2219926/1‑94 (MQ=60)
                                                                                       TGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCTGCCGG        <  SRR4302535.1438121/94‑1 (MQ=60)
                                                                                          CGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCTGCCGGTAA     >  SRR4302535.4605479/1‑94 (MQ=60)
                                                                                          CGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCTGCCGGTAA     >  SRR4302535.3051741/1‑94 (MQ=60)
                                                                                           GCACAGGTAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCTGCCGGTAAC    <  SRR4302535.1819412/94‑1 (MQ=60)
                                                                                           GCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAACCTGGACGTCCCGACGATACTGCCTGCCGGTAAC    <  SRR4302535.2760877/94‑1 (MQ=60)
                                                                                           GCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCTGCCGGTAAC    <  SRR4302535.1206920/94‑1 (MQ=60)
                                                                                           GCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCTGCCGGTAAC    <  SRR4302535.2118526/94‑1 (MQ=60)
                                                                                           GCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCTGCCGGTAAC    >  SRR4302535.3109077/1‑94 (MQ=60)
                                                                                            CACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCTGCCGGTAACG   >  SRR4302535.2288012/1‑94 (MQ=60)
                                                                                             ACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCTGCCGGTAACGC  >  SRR4302535.1667112/1‑94 (MQ=60)
                                                                                             |                                                                                             
CCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCGCAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCTGCCGGTAACGC  >  NC_000913/3987842‑3988028

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