Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 5,321,174 C→T 5.6% V32V (GTG→GTA)  hmuV ← Hemin import ATP‑binding protein HmuV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,321,1740CT 5.6% 241.7 / 3.0 71V32V (GTG→GTA) hmuVHemin import ATP‑binding protein HmuV
Reads supporting (aligned to +/- strand):  ref base C (26/41);  new base T (2/2);  total (28/43)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01

AGGCGCCTTGCCCGTTCCTGCCCGGCCCAATCGGCCAGGGGGCGGCCCTGCAGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCATTGCTGCCCCTGCACAGATACAGGCCCTCAACTTGCAACATCAGCCGCGCCCTCGCATCAGCAGGAACAGGAAGAACGGC  >  ACB122/5321040‑5321308
                                                                                                                                      |                                                                                                                                      
aGGCGCCTTGCCCGTTCCTGCCCGGCCCAATCGGCCAGGGGGCGGCCCTGCAGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCCCAGCACACTCAGCAGGCTGCGCTTGCCGGCACCGTGTGGGCCTAGc                                                                                                                                        >  2:114519/1‑135 (MQ=255)
aGGCGCCTTGCCCGTTCCTGCCCGGCCCAATCGGCCAGGGGGCGGCCCTGCAGCGTTACCCGCCCGCGGTCGGGCGCCAGGTCACCGCCCCGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCGAGc                                                                                                                                        >  2:340011/1‑135 (MQ=255)
  gcgcCTTGCCCGTTCCTGCCCGGCCCAATCGGCCAGGGGGCGGCCCTGCAGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGcac                                                                                                                                      <  5:162348/135‑1 (MQ=255)
        agactcttCCTGCCCTGGCCTCTCGGCCAGGGGGCGGCCCTGCAGCGTTACCCGCCCGCGGGCGGGCGCCAGTTCACCGCACAGCACAGTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCAc                                                                                                                                <  2:204750/129‑1 (MQ=255)
          cccGTTCCTGCCCGGCCCAATCGGCCAGGGGGCGGCCCTGCAGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCAc                                                                                                                                <  6:66348/133‑1 (MQ=255)
           ccGTTCCTGCCCGGCCCAATCGGCCAGGGGGCGGCCCTGCAGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGa                                                                                                                              >  8:268679/1‑134 (MQ=255)
                     cccgCCCAATCGGCCAGGGGGCGGCCCTGCAGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGa                                                                                                                   <  6:104270/132‑1 (MQ=255)
                     ccGGCCCACTCGGCCAGGGGGCGGCCCTGCAGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGa                                                                                                                   <  4:224966/135‑1 (MQ=255)
                     ccGGCCCAATCGGCCAGGGGGCGGCCCTGCAGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACCGCACACTCAGCAGGCGGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGa                                                                                                                   <  6:98790/135‑1 (MQ=255)
                     ccGGCCCAATCGGCCAGGGGGCGGCCCTGCAGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGa                                                                                                                   <  2:33342/135‑1 (MQ=255)
                       ggCCCAATCGGCCAGGGGGCGGCCCTGCCGCGTTACCCGCCCGCGGTCGGCCGCCAGTTCACCCCACACCACCCTCACCCGGCTCCTCTTCCCCGCCCCGTTAGGCCCCATCACACCCACGCCCTCCCCCgccgg                                                                                                                 >  1:61284/1‑131 (MQ=255)
                       ggCCCAATCGGCCAGGGGGCGGCCCTGCAGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAg                                                                                                                  >  1:283231/1‑134 (MQ=255)
                          cgaaTAGGACAGGGGGCGGCCCTGCAGCGTTACCCGCCCGCGGTCGGGCGCCTGTTCACCGCACAGCACACTCTGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACTCCCACGACCTGCCCCGGAGGCaa                                                                                                              <  7:295902/133‑1 (MQ=255)
                          ccAATCGGCCAGGGGGCGGCCCTGCAGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCaa                                                                                                              <  2:306385/135‑1 (MQ=255)
                          ccAATCGGCCAGGGGGCGGCCCTGCAGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCaa                                                                                                              <  7:268679/135‑1 (MQ=255)
                          ccAATCGGCCAGGGGGCGGCCCTGCAGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCaa                                                                                                              <  2:237011/135‑1 (MQ=255)
                           cAATCGGCCAGGGGGCGGCCCTGCAGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGTACACCCACGACCTGCCCCGGAGGCAAc                                                                                                             >  2:272363/1‑135 (MQ=255)
                           cAATCGGCCAGGGGGCGGCCCTGCAGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGTACACCCACGACCTGCCCCGGAGGCAAc                                                                                                             >  2:267386/1‑135 (MQ=255)
                            aaTCGGCCAGGGGGCGGCCCTGCAGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACt                                                                                                            >  8:259866/1‑135 (MQ=255)
                            aaTCGGCCAGGGGGCGGCCCTGCAGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTCCCGGCACCGGTTTGGCCAACCCCACCCACCACCAGCCCCGGAgctggcg                                                                                                            >  1:204750/1‑129 (MQ=255)
                            aaTCGGCCAGGGGGCGGCCCTGCAGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCGCTTGCCGGCACCGTTTGGGCCAAGCCCCCCCCCGACCTGCCCCGGAGGCCACt                                                                                                            >  5:179649/1‑135 (MQ=255)
                               cGGCCAGGGGGCGGCCCTGCAGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCa                                                                                                         <  8:116770/135‑1 (MQ=255)
                               cGGCCAGGGGGCGGCCCTGCAGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCa                                                                                                         <  3:155503/135‑1 (MQ=255)
                                 gCCAGGGGGCGGCCCTGCAGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAgg                                                                                                       <  3:51765/135‑1 (MQ=255)
                                    aGGGGGCGGCCCTGCAGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGTACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGa                                                                                                    <  1:272363/135‑1 (MQ=255)
                                    aGGGGGCGGCCCTGCAGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGTACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGa                                                                                                    <  1:267386/135‑1 (MQ=255)
                                    aGGGGGCGGCCCTGCAGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGa                                                                                                    <  2:90429/135‑1 (MQ=255)
                                    aGGGGGCGGCCCTGCAGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGa                                                                                                    <  6:254395/135‑1 (MQ=255)
                                            gcCCTGCAGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTg                                                                                            <  2:249904/135‑1 (MQ=255)
                                            gcCCTGCAGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTg                                                                                            <  3:232882/135‑1 (MQ=255)
                                             cccTGCAGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGt                                                                                             >  1:17035/1‑133 (MQ=255)
                                              ccTGCAGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTgca                                                                                          >  6:114801/1‑135 (MQ=255)
                                               cTGCAGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTgcag                                                                                         >  5:204249/1‑135 (MQ=255)
                                                   aGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTgcagca                                                                                       <  5:268640/133‑1 (MQ=255)
                                                   aGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCCTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCAcc                                                                                     >  4:171290/1‑135 (MQ=255)
                                                       ttACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCAt                                                                                 <  6:71320/135‑1 (MQ=255)
                                                       ttACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCAt                                                                                 <  2:283231/135‑1 (MQ=255)
                                                       ttACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCAt                                                                                 <  5:290594/135‑1 (MQ=255)
                                                       ttACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCAt                                                                                 <  5:121369/135‑1 (MQ=255)
                                                       ttACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCAt                                                                                 <  2:176589/135‑1 (MQ=255)
                                                       ttACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCAt                                                                                 <  7:346438/135‑1 (MQ=255)
                                                               ccgcGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCATTGCTGccc                                                                         <  4:336665/135‑1 (MQ=255)
                                                                 gcgggggggCGCCAGTTCACCGCACTCCCCACTCAGCCGGCTGCGCTTGCGGGCACCGGTTGGCCCAAGCCCACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCATTGCTGCCCCt                                                                       <  1:185292/129‑1 (MQ=255)
                                                                 gcgGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCATTGCTGCCCCt                                                                       <  3:170052/135‑1 (MQ=255)
                                                                 gcgGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCATTGCTGCCCCt                                                                       <  7:259866/135‑1 (MQ=255)
                                                                  cgGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCATTGCTGCCCCTg                                                                      <  7:236110/135‑1 (MQ=255)
                                                                   ggTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCATTGCTGCCCCTGc                                                                     <  4:27288/135‑1 (MQ=255)
                                                                     tCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCATTGCTGCCCCTGc                                                                     >  1:154745/1‑133 (MQ=255)
                                                                     tCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGGTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCCGCCCCTCATTGCTGCCCCTGc                                                                     >  6:280746/1‑133 (MQ=255)
                                                                                       cacagcacaCTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCATTGCTGCCCCTGCACAGATACAGGCCCTCAACt                                                 <  1:50537/135‑1 (MQ=255)
                                                                                        aaagcacaCTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCATCACCTCATTGCTGCCCCTGCACAGATACAGTCCCTCAACtt                                                >  3:286787/3‑135 (MQ=255)
                                                                                              cacTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTgcagc                                                                                        >  2:201997/1‑89 (MQ=255)
                                                                                              cacTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTgcagc                                                                                        <  1:201997/89‑1 (MQ=255)
                                                                                               acTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCCTTGCTGCCCCTGCACAGATACAGGCCCTCAACTTGCAACa                                          >  2:55444/1‑134 (MQ=255)
                                                                                               acTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCAt                                                                                 >  8:224833/1‑95 (MQ=255)
                                                                                               acTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCAt                                                                                 <  7:224833/95‑1 (MQ=255)
                                                                                                 tCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCATTGCTGCCCCTGCACAGATACAGGCCCTCAACTTGCAACATCa                                       <  1:104725/135‑1 (MQ=255)
                                                                                                  cagcagGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCAc                                                                                                                                <  5:65633/45‑1 (MQ=255)
                                                                                                  cagcagGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCAc                                                                                                                                >  6:65633/1‑45 (MQ=255)
                                                                                                    gcagGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCATTGCTGCCCCTGCACAGATACAGGCCCTCAACTTGCAACATCAg                                      <  2:241890/133‑1 (MQ=255)
                                                                                                       ggCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCAACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCATTGCTGCCCCTGCACAGATACAGGCCCTCAACTTGCAACATCAGCcccg                                 >  8:178469/1‑135 (MQ=255)
                                                                                                         ctgctCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCATTGCTGCCCCTGCACAGCTACAGGCCCTCAACTTGCAACATCAGCCGCGcc                               >  8:281799/1‑135 (MQ=255)
                                                                                                              ctTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCATTGCTGCCCCTGCACAGATACAGGCCCTCAACTTGCAACATCAGCCGCGCCCTc                            >  7:89867/1‑133 (MQ=255)
                                                                                                                tGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCGCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCATTGCTGCCCCTGCACAGATACAGGACCTCAACTTGCAACATCAGCCGCGCCCGCGCa                         >  2:213804/1‑134 (MQ=255)
                                                                                                                  ccGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCATTGCTGCCCCTGCACAGATACAGGCCCTCAACTTGCAACATCAGCCGCGCCCTCGCAt                        <  5:114801/133‑1 (MQ=255)
                                                                                                                    ggCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCATTGCTGCCCCTGCACAGATACAGGCCCTCAACTTGCAACATCAGCCGCGCCCTCGCATca                      >  1:294784/1‑133 (MQ=255)
                                                                                                                         cGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCATTGCTGCCCCTGCACAGATACAGGCCCTCAACTTGCAACATCAGCCGCGCCCTCGCATCAGcagga                >  3:127205/1‑134 (MQ=255)
                                                                                                                              gggCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCATTGCTGCCCCTGCACAGATACAGGc                                                         >  6:33349/1‑88 (MQ=255)
                                                                                                                              gggCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCATTGCTGCCCCTGCACAGATACAGGc                                                         <  5:33349/88‑1 (MQ=255)
                                                                                                                              gggCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCATTGCTGCCCCTGCACAGATACAGGCCCTCAACTTGCAACATCAGCCGCGCCCTCGCATCAGcaggaacag            >  7:173869/1‑133 (MQ=255)
                                                                                                                              gggCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCATTGCTGCCCCTGCACAGATACAGGCCCTCAACTTGCAACATCAGCCGCGCCCTCGCATCAGcaggaacag            >  5:297131/1‑133 (MQ=255)
                                                                                                                                    aGCACACCCACGACCTGCCCCGGAGGCAACTGGAGGTGTATGTCGGGCAGCACCTCATTGCTGCCCCTGCACAGATACAGGCCCTCAACTTGCAACATCAGCCGCGCCCTCGCATCAGCAGGAACAGGAAGAACg    <  2:93984/135‑1 (MQ=255)
                                                                                                                                    aGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCATTGCTGCCCCTGCACAGATACAGGCCCTCAACTTGCAACATCAGCCGCGCCCTCGCATCAGCAGGAACAGGAAGAACg    <  8:173869/135‑1 (MQ=255)
                                                                                                                                    aGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCATTGCTGCCCCTGCACAGATACAGGCCCTCAACTTGCAACATCAGCCGCGCCCTCGCATCAGCAGGAACAGGAAGAACg    <  4:138248/135‑1 (MQ=255)
                                                                                                                                      cacaCCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCAt                                                                                 >  3:56755/1‑56 (MQ=255)
                                                                                                                                      cacaCCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCAt                                                                                 <  4:56755/56‑1 (MQ=255)
                                                                                                                                      cacaCCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCATTGCTGCCCCTGCACAGATACAGGCCCTCAACTTGCAACATCAGCCGCGCCCTCGCATCAGCAGGAACAGGAAGAACGgc  <  3:362392/135‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                      
AGGCGCCTTGCCCGTTCCTGCCCGGCCCAATCGGCCAGGGGGCGGCCCTGCAGCGTTACCCGCCCGCGGTCGGGCGCCAGTTCACCGCACAGCACACTCAGCAGGCTGCTCTTGCCGGCACCGTTTGGGCCAAGCACACCCACGACCTGCCCCGGAGGCAACTGCAGGTGAATGTCGTGCAGCACCTCATTGCTGCCCCTGCACAGATACAGGCCCTCAACTTGCAACATCAGCCGCGCCCTCGCATCAGCAGGAACAGGAAGAACGGC  >  ACB122/5321040‑5321308

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑