Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 317,566 G→A 17.6% pseudogene (128/129 nt) ykgQ → pseudogene, putative dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913317,5660GA17.6% 225.4 / 36.4 107pseudogene (128/129 nt)ykgQpseudogene, putative dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (39/47);  new base A (10/9);  total (49/58)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.98e-02

Only 100 of 119 total aligned reads displayed.

ATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGC  >  NC_000913/317472‑317664
                                                                                              |                                                                                                  
attattTCGTTTTATTTGGGTTATTTTGGTGAGTCTGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATAtt                                                                                               <  1:583261/100‑1 (MQ=255)
aTTAGTTCGCTTGATTTGAGTTATGTTGGGGAGACGGTAAAACATTTTAATCCTGCGTTTGTTCGTGGTTAAATAAAAAAAAAAATACATGATTGATAtt                                                                                               <  2:1483895/100‑1 (MQ=255)
    gTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAGTTTGTTCGTAGTTGTGTTGCTGTAAAAATACATGATTGATATTAATc                                                                                           >  2:2047193/1‑100 (MQ=255)
     ttCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTTGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCa                                                                                          <  2:1180914/100‑1 (MQ=255)
       cGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATg                                                                                        <  1:281536/100‑1 (MQ=255)
       cGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATg                                                                                        <  1:2853931/100‑1 (MQ=255)
        gCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGt                                                                                       >  2:2654702/1‑100 (MQ=255)
         gttGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCAtat                                                                                       >  2:646586/2‑97 (MQ=255)
          ttGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAt                                                                                     <  1:2673432/100‑1 (MQ=255)
          ttGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAATACATGATTAATATCAATCATGTAtt                                                                                    >  1:65320/1‑100 (MQ=38)
           tGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAtt                                                                                    <  2:1586678/100‑1 (MQ=255)
             aTTTGGGTTATGTTGGGGAGTCGGTAAAATATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAtttt                                                                                  <  2:1759071/100‑1 (MQ=255)
              tttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAttttt                                                                                 <  2:2443470/100‑2 (MQ=255)
              tttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAttttt                                                                                 <  2:1979549/100‑2 (MQ=255)
               ttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTAc                                                                                <  1:1206117/100‑1 (MQ=255)
                tGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAtttttac                                                                               <  2:102029/100‑4 (MQ=38)
                 gggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCt                                                                              >  2:843987/1‑100 (MQ=255)
                 gggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCt                                                                              >  1:2723939/1‑100 (MQ=255)
                 gggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAtttttaca                                                                              <  2:1963988/100‑5 (MQ=38)
                  ggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAGATACATGATTGATATTAATCATGTATTTTACCTa                                                                             >  2:1594497/1‑100 (MQ=255)
                  ggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAtttttacag                                                                             <  1:1902748/100‑6 (MQ=37)
                     tATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAga                                                                          <  2:2503384/100‑1 (MQ=255)
                     tATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAga                                                                          <  1:364792/100‑1 (MQ=255)
                     tATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAtttttacagcaa                                                                          <  1:636864/100‑9 (MQ=25)
                      aTGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAgac                                                                         <  1:458371/100‑1 (MQ=255)
                      aTGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAgac                                                                         <  1:1155201/100‑1 (MQ=255)
                         ttGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAgacgac                                                                      >  2:2316474/1‑100 (MQ=255)
                           ggggAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACaa                                                                    <  1:127730/100‑1 (MQ=255)
                           ggggAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACaa                                                                    >  1:2205657/1‑100 (MQ=255)
                            gggAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACaaa                                                                   >  1:1051223/1‑100 (MQ=255)
                            gggAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACaaa                                                                   >  1:1827465/1‑100 (MQ=255)
                             ggAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAAACATGTATTTTACCTAAGACGACaaaa                                                                  >  1:2383248/1‑100 (MQ=255)
                                gTCGGTAAAACATTTTAATCCTGAATTTGTTTGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAAcc                                                               <  1:306463/100‑1 (MQ=255)
                                gTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAAcc                                                               >  2:7947/1‑100 (MQ=255)
                                gTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAAcc                                                               <  2:254283/100‑1 (MQ=255)
                                gTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAAc                                                                >  1:957595/1‑99 (MQ=255)
                                  cGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCtt                                                             >  1:2407672/1‑100 (MQ=255)
                                     tAAAATGTTTTAATCCTGAATTTGGTCGTCGTTTTAGTGCTGTAAAAATACATGTTTGATTTTAATGATGTATTTTTCCGAAGACGACAAAAACCTTTTa                                                          <  2:820746/100‑1 (MQ=255)
                                     tAAAACATTTTAATCCTGTATTTGTTCGTTGTTGACTTGCTATAAAAGTACACTAGTAATATTAATCCAGTATTTTAACTATCCCGACaaaaccattgt                                                           <  2:1321250/99‑8 (MQ=255)
                                      aaaaCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTaa                                                         >  2:2691432/1‑100 (MQ=255)
                                       aaaGATCTTAATTCTGAATCTGTGCGTGTTTGTCTTGCTGTAAAAATACATGATTGATATTATTCATGTATTTTACCTAAGGCGACAAAAACCt                                                              <  2:2975283/94‑1 (MQ=255)
                                        aaCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACa                                                       >  1:1726923/1‑100 (MQ=255)
                                        aaCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACa                                                       >  1:2039667/1‑100 (MQ=255)
                                        aaCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACa                                                       >  1:147296/1‑100 (MQ=255)
                                         aCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAg                                                      <  1:2444160/100‑1 (MQ=255)
                                         aCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGAAGACAAAAACCTTTTAACAg                                                      <  2:150390/100‑1 (MQ=255)
                                          cATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAgg                                                     >  1:158120/1‑100 (MQ=255)
                                           aTTTTAATCCTGTATTTGTTCGTGTTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCCAAGACGACAAAAACCTTTTAACAggg                                                    <  2:802647/100‑1 (MQ=255)
                                           aTTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAggg                                                    <  2:2594236/100‑1 (MQ=255)
                                           aTTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAgg                                                     <  1:1480570/99‑1 (MQ=255)
                                             tttAATCCTGAATTTGTTCGTGGTTTTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAAAAGAAAAAAACCTTTTAACa                                                       <  2:2911490/95‑1 (MQ=255)
                                             tttAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAg                                                  <  1:1346177/100‑1 (MQ=255)
                                             tttAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAg                                                  <  1:2075790/100‑1 (MQ=255)
                                                 aTCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGat                                              <  1:795272/100‑1 (MQ=255)
                                                  tCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGata                                             >  1:2071463/1‑100 (MQ=255)
                                                  tCCTGAATTTGTTCGTGGTTGTATTGCTGTAAAAATACATGATTGATATTAATGATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGACTGgta                                             >  2:2658573/1‑100 (MQ=255)
                                                   ccTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGatat                                            >  1:742122/1‑100 (MQ=255)
                                                   ccTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGatag                                            >  1:1976608/1‑99 (MQ=255)
                                                   ccTGAATTTGTTCGTGGTTGTCTTGATGTAAAAATACATGATTGATATTAATCATGTATTATACATAATACGACAAAAACCTTTTAATAGGGAGTGatat                                            >  2:1283847/1‑100 (MQ=255)
                                                     tGAATTTGTTCGTGGTTGTCTTGCAGTAAAAATACATGATTTATATTTATCATGTATTTTACAAAAAAAAAAAAAAACCTTTTAACAGGGa                                                   <  2:68679/91‑1 (MQ=255)
                                                         tttGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCt                                      >  2:1053626/1‑100 (MQ=255)
                                                          ttGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCtt                                     <  2:863107/100‑1 (MQ=255)
                                                            gTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCa                                   <  1:831640/100‑1 (MQ=255)
                                                              tCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCAct                                 <  1:2974220/100‑1 (MQ=255)
                                                              tCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCAct                                 <  1:1594516/100‑1 (MQ=255)
                                                               catGGTTGTATTGATGTAAAAATACATGATTTATATTAATCATGTATTTTACCTTAGACGACAAAAGCCTTTTGTCAGGGAGTGATATTGATCTTCAc                                  <  2:2596794/96‑1 (MQ=255)
                                                               cGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCActc                                <  2:1065126/100‑1 (MQ=255)
                                                               cGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCAct                                 >  2:573833/1‑99 (MQ=255)
                                                                 tGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTg                              >  1:2149232/1‑100 (MQ=255)
                                                                  ggTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGt                             <  2:2143832/100‑1 (MQ=255)
                                                                   gTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTc                            <  2:2485952/100‑1 (MQ=255)
                                                                   gTTGTCTTGCTATAAAAATACATGATTGATATTAATCATGTATTTTACCGTAGACGATAAAAACTTTTTAACATGTAGTGTTTTTGATCTTCACTCTGt                             <  2:2767261/99‑1 (MQ=255)
                                                                      gTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCata                         <  2:1163648/100‑1 (MQ=255)
                                                                        cTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATAtc                       >  2:1640032/1‑100 (MQ=255)
                                                                           gCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTcc                    <  2:2149209/100‑1 (MQ=255)
                                                                           gCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTcc                    >  2:100773/1‑100 (MQ=255)
                                                                            cTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCg                   >  2:1948787/1‑100 (MQ=255)
                                                                            cTGTAAAAATACATGATTGATATTAATCATGTATTTTAACTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTcc                    >  2:917973/1‑99 (MQ=255)
                                                                            cTGTAAAAATACATGATTAATATTAATCATGTATTTTACCTAAGACGAAAAAAAGCTTTTAACAGGTAGTGATATTGATCTTCAATCTGTCATAtc                       >  2:2592839/1‑96 (MQ=38)
                                                                             tGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATAtctc                     >  2:1746539/1‑97 (MQ=255)
                                                                             tGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGAAAGTGATATTGATCTTCACTCTGTAATATCTCCg                   >  1:302391/1‑99 (MQ=255)
                                                                              ggtaaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGt                 >  2:2181582/4‑100 (MQ=38)
                                                                              ggtaaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGt                 >  1:1611877/4‑100 (MQ=38)
                                                                              ggtaaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTAACTCTGTCATATCTCCg                   >  2:2823484/4‑98 (MQ=37)
                                                                              gTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGt                 >  1:1117225/1‑100 (MQ=255)
                                                                               gtaaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTa                >  1:2302451/3‑100 (MQ=255)
                                                                               gtaaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTa                >  1:1618125/3‑100 (MQ=255)
                                                                                 aaaaTACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAAAAGGGAGTGATATTGATATTAACTATGTAATATATACGGTaa               >  1:2705010/1‑99 (MQ=255)
                                                                                    aTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGCTATTGATCTTCACTCTGTCATATCTCCGGTAATATg           >  1:2130052/1‑100 (MQ=39)
                                                                                    aTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAatat            >  1:1231252/1‑99 (MQ=255)
                                                                                     tACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATgg          >  1:325293/1‑100 (MQ=255)
                                                                                      aCATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGc         <  1:1435262/100‑1 (MQ=255)
                                                                                       cATGATTGATATTAATCATGTATTTTACCTAAGTCGATAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCg        <  1:1393992/100‑1 (MQ=255)
                                                                                          gattgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCa     <  2:1601269/100‑1 (MQ=255)
                                                                                          gattgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCa     >  2:1006586/1‑100 (MQ=255)
                                                                                          gattgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATAAATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCa     >  1:2471570/1‑100 (MQ=255)
                                                                                          gattgatATTAAGCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCa     >  1:1995572/1‑100 (MQ=255)
                                                                                           attgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAg    >  2:588161/1‑100 (MQ=255)
                                                                                            ttgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTATGTCATATCTCCGGTAATATGGCGTCAgg   >  1:186087/1‑100 (MQ=255)
                                                                                             tgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGc  <  2:395288/100‑1 (MQ=255)
                                                                                              |                                                                                                  
ATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGC  >  NC_000913/317472‑317664

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