Predicted mutation
evidence seq id position mutation annotation gene description
RA CP016816 313,637 C→T E350K (GAA→AAA)  JCVSYN2_01315 ← peptidase M20

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP016816313,6370CT99.5% 11579.8 / NA 3588E350K (GAA→AAA) JCVSYN2_01315peptidase M20
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (1837/1722);  minor base A (12/6);  total (1856/1732)
Fisher's exact test for biased strand distribution p-value = 2.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 3708 total aligned reads displayed.

TTTTGCTTTATAATCTTGAGTTACTTCTTGATAAACTCTCATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTCTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAGAATTTTTATTTTTAATTTCAATAAGTCCTACGTTTTGAGTTAAT  >  CP016816/313512‑313761
                                                                                                                             |                                                                                                                            
ttttGCTTTATAATCTTGAGTTACTTCTTGATAAACTCTCATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTa                                                                                                                            >  2:6259264/1‑128 (MQ=255)
 tttGCTTTATAATCTTGAGTTACTTCTTGATAAACTCTCATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTaa                                                                                                                           >  1:6834090/1‑128 (MQ=255)
    gCTTTATAATCTTGAGTTACTTCTTGATAAACTCTCATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAata                                                                                                                        <  1:3623207/128‑1 (MQ=255)
      tttATAATCTTGAGTTACTTCTTGATAAACTCTCATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATAtt                                                                                                                      >  2:5397271/1‑128 (MQ=255)
         ataATCTTGAGTTACTTCTTGATAAACTCTCATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATttag                                                                                                                   >  1:4453852/1‑128 (MQ=255)
          taATCTTGAGTTACTTCTTGATAAACTCTCATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATttagt                                                                                                                  >  2:612540/1‑128 (MQ=255)
            aTCTTGAGTTACTTCTAGATAAACTCTCATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATttagtta                                                                                                                <  2:201914/128‑1 (MQ=255)
              cTTGAGTTACTTCTTGATAAACTCTCATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTtagttagt                                                                                                              >  1:6588087/1‑128 (MQ=255)
                tGAGTTACTTCTTGATTAACTCTCATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAgtgt                                                                                                            <  1:5353498/128‑1 (MQ=255)
                tGAGTTACTTCTTGATAAACTCTCATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAgtgt                                                                                                            <  1:3782779/128‑1 (MQ=255)
                 gAGTTACTTCTTGATAAACTCTCATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGtt                                                                                                           <  1:7718359/128‑1 (MQ=255)
                  aggtactacTTGATAAACTCTCATAATTTTATTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTg                                                                                                          <  2:1865161/125‑1 (MQ=255)
                    ttACTTCTTGATAAACTCTCATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTggg                                                                                                        >  1:100159/1‑128 (MQ=255)
                      aCTTCTTGATAAACTCTCATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGAt                                                                                                      <  2:664570/128‑1 (MQ=255)
                      aCTTCTTGATAAACTCTCATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGAt                                                                                                      >  2:2769554/1‑128 (MQ=255)
                       gatcttGATAAACTCTCATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATa                                                                                                     <  2:3944677/126‑1 (MQ=255)
                       cttcttGATAAACTCTCATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATa                                                                                                     >  1:5208186/1‑128 (MQ=255)
                         tcttGATAAACTCTCATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATaaa                                                                                                   <  2:3482535/128‑1 (MQ=255)
                         tcttGATAAACTCTCATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATaaa                                                                                                   >  2:2477241/1‑128 (MQ=255)
                         tcttGATAAACTCTCATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATaaa                                                                                                   <  1:4068752/128‑1 (MQ=255)
                           ttGATAAACTCTCATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAata                                                                                                 >  1:7588732/1‑128 (MQ=255)
                            tGATAAACTCTCATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAatat                                                                                                >  1:3065517/1‑128 (MQ=255)
                             gATAAACTCTCATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATAtt                                                                                               <  2:4453852/128‑1 (MQ=255)
                             gATAAACTCTCATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATAtt                                                                                               >  2:413059/1‑128 (MQ=255)
                                  aCTCTCATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGt                                                                                          <  2:4639603/128‑1 (MQ=255)
                                    tctcATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCATTTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTg                                                                                        <  1:7847502/128‑1 (MQ=255)
                                    tctcATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTg                                                                                        <  1:6975151/128‑1 (MQ=255)
                                    tctcATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTg                                                                                        <  1:6430321/128‑1 (MQ=255)
                                    tctcATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTg                                                                                        >  1:1440746/1‑128 (MQ=255)
                                     ctcATTATATTTTTTATTAAATCTGATTCTTGTTGAACTTATTAGCTATTATCATTTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAATTATTTAGTTAGTGTTGGGATAAATATTTGAGTTgg                                                                                       <  1:1260413/128‑1 (MQ=255)
                                     ctcATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTATATATTTAGTTAGTGTTGGGATAAATATTTGAGTTgg                                                                                       <  1:1577334/128‑1 (MQ=255)
                                      tcATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTTGGATAAATATTTGAGTTggt                                                                                      >  2:1456465/1‑127 (MQ=255)
                                      tcATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTggg                                                                                      >  2:7286760/1‑128 (MQ=255)
                                      tcATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACCTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGTGATAAATATTTGAGTTgtg                                                                                      >  2:5801356/1‑126 (MQ=255)
                                       cATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGt                                                                                     <  2:2155312/128‑1 (MQ=255)
                                         tAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTTGTTTa                                                                                   >  1:2730808/1‑128 (MQ=255)
                                           aTTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGTTAAATATTTGATTTGTGTTAGt                                                                                 >  2:3653238/1‑128 (MQ=255)
                                           aTTTTTTAATTTAAAGCTTATTCTTCATTAAAATATATCTAATTATCAATTTTTTAAACTT‑TAATTTATATTTTATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTa                                                                                <  1:4647824/128‑1 (MQ=255)
                                            tttttttAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTTGTa                                                                                >  1:5059046/1‑128 (MQ=255)
                                            tttttttAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTa                                                                                >  2:3353293/1‑128 (MQ=255)
                                            tttttttAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTa                                                                                >  1:7095007/1‑128 (MQ=255)
                                              tttttAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGt                                                                                          <  1:5731585/116‑1 (MQ=255)
                                               ttttAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAAt                                                                             <  1:3493583/128‑1 (MQ=255)
                                                    attaaATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAgg                                                                        >  1:6325618/1‑128 (MQ=255)
                                                     ttaaATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGa                                                                       <  2:7168951/128‑1 (MQ=255)
                                                       aaaTCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAAt                                                                     <  1:2100089/128‑1 (MQ=255)
                                                         aTCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTc                                                                   >  1:2260806/1‑128 (MQ=255)
                                                           cTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTa                                                                 <  2:6397947/128‑1 (MQ=255)
                                                             gATTCTTTATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTaaa                                                               <  2:833330/128‑1 (MQ=255)
                                                             gATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTaaa                                                               >  1:1834385/1‑128 (MQ=255)
                                                              aTTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACATTAATTATAAAg                                                              >  2:475642/1‑128 (MQ=255)
                                                                 cTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGttt                                                           <  2:2299910/128‑1 (MQ=255)
                                                                  ttGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTa                                                          <  2:3199485/128‑1 (MQ=255)
                                                                  ttGATGAACATATAAGATATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTa                                                          <  2:1758496/128‑1 (MQ=255)
                                                                   tgatgaACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTTAGTTGGGTTAGTAAATACAGTAATTCTAAAGTTTaa                                                         >  2:5047147/1‑128 (MQ=255)
                                                                   tgatgaACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGTGTTAGTAAATACAGGAATTCTAAATTTTaa                                                         >  2:4918500/1‑128 (MQ=255)
                                                                   tgatgaACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGATGGTTTAGTAAATACATGAATTCTAAAGTTTaa                                                         >  2:5977315/1‑128 (MQ=255)
                                                                    gatgaACATATAAGCTATTATCACTTTTATTAACTTCTAATGATAAGTAAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTTGGTTAGTAAATACAGGAATTCTAAAGTTTAAt                                                        <  1:4789729/128‑1 (MQ=255)
                                                                    gatgaACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAAt                                                        <  2:1156584/128‑1 (MQ=255)
                                                                    gatgaACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAAAACAGGAATTCTAAAGTTTAAt                                                        >  1:4347616/1‑128 (MQ=255)
                                                                     atgaACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATc                                                       >  2:5678260/1‑128 (MQ=255)
                                                                      tgaACATATAAGGTTTTATCTATTTTTTTAAATTCTAATTATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATcc                                                      <  2:7604181/128‑1 (MQ=255)
                                                                       gaACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCt                                                     <  2:6904988/128‑1 (MQ=255)
                                                                         aCATATAATCTATTATCAATTTTTTTAACTTCTAATGATAATTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGATTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTAt                                                   <  2:2618443/128‑1 (MQ=255)
                                                                         aCATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATACTAAATGTTAATCCTAt                                                   >  1:4202069/1‑128 (MQ=255)
                                                                          cATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTAtt                                                  >  1:292067/1‑128 (MQ=255)
                                                                           atatAAGCTATTATCAATTTTTTTCACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTc                                                 <  2:5546321/128‑1 (MQ=255)
                                                                            tataAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCt                                                <  1:2562220/128‑1 (MQ=255)
                                                                             ataAGCTATTATCAATTTTTTTAACTTCTAATGATCAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTa                                               <  2:7867451/128‑1 (MQ=255)
                                                                              taAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAg                                              >  2:1023677/1‑128 (MQ=255)
                                                                               aaGCTATTATTTATTTTTTTTACTTCTAATGATAAGGTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAGa                                             <  2:1616317/128‑1 (MQ=255)
                                                                                  cTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAGAatt                                          >  2:7052025/1‑128 (MQ=255)
                                                                                  cTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTAAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAGAatt                                          <  2:5259491/128‑1 (MQ=255)
                                                                                   tattatGAATTTTTTAAACTTCTAATGATATGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAGAattt                                         <  1:3200198/128‑1 (MQ=255)
                                                                                    attatCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAGAatttt                                        >  1:844259/1‑128 (MQ=255)
                                                                                          aaTTTTTTTCACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAGAatttttatttt                                  <  2:3809246/128‑1 (MQ=255)
                                                                                           aTTTTTTTAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAGAatttttattttt                                 >  2:2920188/1‑128 (MQ=255)
                                                                                              tttttAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGTGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAGAATTTTTATTTTTAAt                              <  2:6782358/128‑1 (MQ=255)
                                                                                              tttttAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAGAATTTTTATTTTTAAt                              <  2:429306/128‑1 (MQ=255)
                                                                                               ttttAACTTCTAATGATAAGTTAATTTGTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAGAATTTTTATTTTTAAtt                             <  1:7448365/128‑1 (MQ=255)
                                                                                                tttAACTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAGAATTTTTATTTTTAAttt                            >  2:7730913/1‑128 (MQ=255)
                                                                                                attAACTTATCATTAGAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAGAATTTTTATTTTTAAttt                            >  1:5627742/2‑128 (MQ=255)
                                                                                                   aaCTTATAATGATAAGTTAATTTTTTTTAAATATTAAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAGAATTTTTATTTTTAATTTCaa                         <  1:4928315/128‑1 (MQ=255)
                                                                                                    aCTTCTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAGAATTTTTATTTTTAATTTCAAt                        <  1:3357536/128‑1 (MQ=255)
                                                                                                     cTTCTACTGATAATTTAATTTTTTTTAAATATTTAGGTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAGAATTTTTATTTTTAATTTCAATa                       <  1:6713979/128‑1 (MQ=255)
                                                                                                     cTTCTAATGATAAGTTAATTGTTTTTAAATATTTAGTTAGTGGTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAGAATTTTTATTTTTAATTTCAATa                       <  1:6009656/128‑1 (MQ=255)
                                                                                                      ttCTAATGATAAGATAAGTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATATGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAGAATTTTTATTTTTAATTTCAATaa                      <  2:3074326/128‑1 (MQ=255)
                                                                                                        cTAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAGAATTTTTATTTTTAATTTCAATAAGt                    <  1:1681566/128‑1 (MQ=255)
                                                                                                         tAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAGAATTTTTATTTTTAATTTCTATAAGTc                   >  1:3906432/1‑128 (MQ=255)
                                                                                                         tAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAGAATTTTTATTTTTAATTTCAATAAGTc                   >  1:1124440/1‑128 (MQ=255)
                                                                                                         tAATGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAGAATTTTAATTTTTAATTTCAATAAGTc                   >  1:1972001/1‑128 (MQ=255)
                                                                                                          aaTGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAGAATTTTTATTTTTAATTTCAATAAGTcc                  >  1:244334/1‑128 (MQ=255)
                                                                                                          aaTGATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAGAATTTTTATTTTTAATTTCAATAAGTcc                  <  1:7337260/128‑1 (MQ=255)
                                                                                                             gATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAGAATTTTTATTTTTAATTTCAATAAGTCCTAc               <  1:5473097/128‑1 (MQ=255)
                                                                                                             gATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAGAATTTTTATTTTTAATTTCAATAAGTCCTAc               <  2:7411288/128‑1 (MQ=255)
                                                                                                             gATAAGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAGAATTTTTATTTTTAATTTCAATAAGTCCTAc               <  2:1296947/128‑1 (MQ=255)
                                                                                                                aaGTTAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAGAATTTTTATTTTTAATTTCAATAAGTCCTACGtt            >  1:5860657/1‑128 (MQ=255)
                                                                                                                    tAATTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAGAATTTTTATTTTTAATTTCAATAAGTCCTACGTTTTGa        >  2:5147461/1‑128 (MQ=255)
                                                                                                                      aTTTTTTTTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAGAATTTTTATTTTTAATTTCAATAAGTACTACTTTTAGAGt      >  1:7204765/1‑128 (MQ=255)
                                                                                                                       ttttttttAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAGAATTTTTATTTTTAATTTCAATAAGTCCTACGTTTTGAGtt     >  1:6188323/1‑128 (MQ=255)
                                                                                                                             |                                                                                                                            
TTTTGCTTTATAATCTTGAGTTACTTCTTGATAAACTCTCATAATTTTTTTAATTAAATCTGATTCTTGATGAACATATAAGCTATTATCAATTTTTTTAACTTCTAATGATAAGTTAATTTTTTCTAAATATTTAGTTAGTGTTGGGATAAATATTTGAGTTGGGTTAGTAAATACAGGAATTCTAAAGTTTAATCCTATTCTAGAATTTTTATTTTTAATTTCAATAAGTCCTACGTTTTGAGTTAAT  >  CP016816/313512‑313761

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