Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,986,892 C→T H123H (CAC→CAT yqeI → putative sensory transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,986,8920CT98.8% 582.1 / NA 166H123H (CAC→CATyqeIputative sensory transducer
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (115/49);  minor base A (0/2);  total (115/51)
Fisher's exact test for biased strand distribution p-value = 9.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 226 total aligned reads displayed.

TGCAACCGTGCCGAGAAGAGGATTCAAGATCCATAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACACAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGAATTTCAATAACAAAATCCTCAAACATAT  >  NC_000913/2986802‑2986981
                                                                                          |                                                                                         
tGCAACCGTGCCGAGAAGAGGATTCAAGATCCATAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATaa                                                                                         <  2:2492058/93‑1 (MQ=255)
tGCAACCGTGCCGAGAAGAGGATTCAAGATCCATAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATaa                                                                                         <  1:3723457/93‑1 (MQ=255)
 gCAACCGTGCCGAGAAGAGGATTCAAGATCCATAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAAt                                                                                        <  2:876548/93‑1 (MQ=255)
 gCAACCGTGCCGAGAAGAGGATTCAAGATCCATAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAAt                                                                                        <  2:403975/93‑1 (MQ=255)
 gCAACCGTGCCGAGAAGAGGATTCAAGATCCATAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAAt                                                                                        <  2:2513932/93‑1 (MQ=255)
 gCAACCGTGCCGAGAAGAGGATTCAAGATCCATAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAAt                                                                                        >  1:1916474/1‑93 (MQ=255)
       gTGCCGAGAAGAGGATTCAAGATCCATAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCa                                                                                  <  2:1122485/93‑1 (MQ=255)
           cGAGAAGAGGATTCAAGATCCATAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTg                                                                              >  1:1061695/1‑93 (MQ=255)
            gagaAGAGGATTCAAGATCCATAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGc                                                                             >  2:2015578/1‑93 (MQ=255)
            gagaAGAGGATTCAAGATCCATAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGc                                                                             >  2:118307/1‑93 (MQ=255)
            gagaAGAGGATTCAAGATCCATAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGc                                                                             <  1:2523506/93‑1 (MQ=255)
               aagaGGATTCAAGATCCATAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCAt                                                                          >  1:4533434/1‑93 (MQ=255)
               aagaGGATTCAAGATCCATAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCAt                                                                          >  2:1258315/1‑93 (MQ=255)
                agagGATTCAAGATCCATAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATa                                                                         >  1:862572/1‑93 (MQ=255)
                agagGATTCAAGATCCATAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATa                                                                         >  2:4456563/1‑93 (MQ=255)
                agagGATTCAAGATCCATAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATa                                                                         >  1:761588/1‑93 (MQ=255)
                agagGATTCAAGATCCATAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATa                                                                         >  1:3336172/1‑93 (MQ=255)
                 gagGATTCAAGATCCATAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATaa                                                                        >  2:634748/1‑93 (MQ=255)
                  agGATTCAAGATCCATAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATaaa                                                                       >  1:1362156/1‑93 (MQ=255)
                   ggATTCAAGATCCATAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATaaaa                                                                      <  2:2612556/93‑1 (MQ=255)
                     aTTCAAGATCCATAATGACATAAACATACAAAATCATGAAATAACCTCGTCAACAGACGCACATACACATAATGCCCCCCCTGCCATAAAAGt                                                                    <  2:2424684/93‑1 (MQ=255)
                      ttCAAGATCCATAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGtt                                                                   >  2:4297519/1‑93 (MQ=255)
                       tCAAGATCCATAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTa                                                                  >  2:1649106/1‑93 (MQ=255)
                         aaGATCCATAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCCCCTGCCATAAAAGTTAAt                                                                <  1:2624688/93‑1 (MQ=255)
                         aaGATCCATAATGACATAAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAAt                                                                <  2:953046/93‑1 (MQ=255)
                          aGATCCATAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATg                                                               >  1:1095233/1‑93 (MQ=255)
                             tCCATAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCgg                                                            >  2:3237113/1‑93 (MQ=255)
                             tCCATAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCgg                                                            >  2:4362660/1‑93 (MQ=255)
                             tCCATAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCgg                                                            >  1:1055078/1‑93 (MQ=255)
                              ccATAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCggg                                                           >  1:2995559/1‑93 (MQ=255)
                                aTAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCCCCTGCCATAAAAGTTAATGCGGGAt                                                         <  1:2733168/93‑1 (MQ=255)
                                aTAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGAt                                                         >  1:3501489/1‑93 (MQ=255)
                                 tAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATa                                                        >  2:2330235/1‑93 (MQ=255)
                                 tAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATa                                                        >  1:950830/1‑93 (MQ=255)
                                  aaTGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATAc                                                       <  1:695716/93‑1 (MQ=255)
                                   aTGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACa                                                      >  1:1308287/1‑93 (MQ=255)
                                    tGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACCTACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACaa                                                     >  2:2918355/1‑93 (MQ=255)
                                     aaCATTAACATACAAAATCTTGAAATAAACTCGTCAACAGACGCAAATACACAAAACGCCCAACCTGCCATAAAAGTTAATGCGGGATACaaa                                                    <  2:1403756/92‑1 (MQ=255)
                                        ataaaCATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCCCCTGCCATAAAAGTTAATGCGGGATACAAagag                                                 <  2:2859369/90‑1 (MQ=255)
                                        aTTACCATACAAAATCATGTAATAATCTCGTCAACAGCCGCCCATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAagag                                                 <  1:4061238/93‑1 (MQ=255)
                                          tAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCCCCTGCCATAAAAGTTAATGCGGGATACAAagagag                                               <  1:236864/93‑1 (MQ=255)
                                          tAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAagagag                                               >  1:3849785/1‑93 (MQ=255)
                                           aaCATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGc                                              <  1:3647824/93‑1 (MQ=255)
                                            aCATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCa                                             >  2:756468/1‑93 (MQ=255)
                                             cATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCAt                                            >  2:227587/1‑93 (MQ=255)
                                              aTACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCAtt                                           >  1:4032158/1‑93 (MQ=255)
                                              aTACAAAATATTGTAATAACCTCGTAAACAGACGCACATACACATAATGCCCCCCCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCAtt                                           <  2:3402447/93‑1 (MQ=255)
                                              aTAAAAAATCATGAAATAACCTCGTCAACAGCCGCACATACACATAATCCCCCCCCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCAtt                                           <  2:1669396/93‑1 (MQ=255)
                                               tACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATtg                                          >  2:1533624/1‑93 (MQ=255)
                                                aCAAAATCATGTAATAACCTCGTCAACAGACGCACATACACATAATGCCCTACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATtgg                                         <  1:3301112/93‑1 (MQ=255)
                                                 cAAAATCAGCTAATAACCTCGTCAACAGCCGCACATACACATAATCCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATtggt                                        <  2:563081/93‑1 (MQ=255)
                                                 aaaaaTCTTGTAATAAACTCGTCAACAGACGCACATACACATAATCCCCCCCCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATtggt                                        <  2:2182454/92‑1 (MQ=255)
                                                 aaaaaTCATGTAATAAACTCGTCAACAGCCGCACATACACATAATGCCCCCCCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATtggt                                        <  1:2552551/92‑1 (MQ=255)
                                                  aaaaTCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATtggtg                                       >  1:1239570/1‑93 (MQ=255)
                                                   aactCATGAAAAAAACGCGTCAACAGACGCACATACACATAATGCCCCCCCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATtggtgg                                      <  1:547122/90‑1 (MQ=255)
                                                   aaaaCATGTAATAACCTCGTAAACAGCCGCACAAACACATATCGCACGTCCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATtggtgg                                      <  1:1533449/93‑1 (MQ=255)
                                                   aaaTCATGTAATAAGCTCGTCAACAGACACACATACACATAATGCCCCCCCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATtggtgg                                      <  1:798962/93‑1 (MQ=255)
                                                   aaaTCATGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATtggtgg                                      <  2:4316391/93‑1 (MQ=255)
                                                        aTGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCaa                                 >  1:1661452/1‑93 (MQ=255)
                                                        aTGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCaa                                 >  1:3144163/1‑93 (MQ=255)
                                                         tGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCaaa                                >  1:4346876/1‑93 (MQ=255)
                                                         tGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCaaa                                >  1:44508/1‑93 (MQ=255)
                                                         tGTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCaaa                                >  2:3861044/1‑93 (MQ=255)
                                                          gTAATAACCTCGTCAGAAGACGCACATACACAAAACGCCCCTCCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAg                               <  1:1421288/93‑1 (MQ=255)
                                                          gTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAg                               <  1:4283760/93‑1 (MQ=255)
                                                          gTAATAAACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAg                               >  1:4148171/1‑93 (MQ=255)
                                                            aataaCCTCGTCAACAGACGCACATACACATAAGCCCCCCCCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGaa                             <  1:2261737/93‑1 (MQ=255)
                                                            aataaACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGaa                             >  2:1592286/1‑93 (MQ=255)
                                                            aataaACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGaa                             >  1:2312671/1‑93 (MQ=255)
                                                            aataaACGCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGaa                             <  2:1778150/93‑1 (MQ=255)
                                                             ataaACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGAAt                            >  1:4225701/1‑93 (MQ=255)
                                                              taaACTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGAAtt                           >  1:1442595/1‑93 (MQ=255)
                                                                accTCGCCAACAGACCCACATACACATAATGCCCCCCCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGAATTTc                         <  2:3442094/91‑1 (MQ=255)
                                                                aaCTCGTCAACAGACGCACATACACATAATGCCCCCCCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGAATTTc                         <  2:2662948/93‑1 (MQ=255)
                                                                  cTCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGAATTTCaa                       >  1:3222974/1‑93 (MQ=255)
                                                                   tCGTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGAATTTCAAt                      >  2:3165067/1‑93 (MQ=255)
                                                                    cGTCAACAGACGCACATACACATAACGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGAATTTCAATa                     <  1:2086432/93‑1 (MQ=255)
                                                                     gTCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGAATTTCAATaa                    >  2:2686751/1‑93 (MQ=255)
                                                                      tCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGAATTTCAATAAc                   >  2:2898134/1‑93 (MQ=255)
                                                                      tCAACAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGAATTTCAATAAc                   >  2:1238788/1‑93 (MQ=255)
                                                                        aaCAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGAATTTCAATAACaa                 >  1:4675983/1‑93 (MQ=255)
                                                                        aaCAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGAATTTCAATAACaa                 >  2:2833994/1‑93 (MQ=255)
                                                                         aCAGACCCACATACACAAAACGCCCCCCCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGAATTTCAATAACaaa                <  2:4107524/93‑1 (MQ=255)
                                                                          cAGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGAATTTCAATAACaaaa               <  1:4360618/93‑1 (MQ=255)
                                                                          cAGACGCACATACACATAATCCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGAATTTCAATAACaaaa               <  1:3443080/93‑1 (MQ=255)
                                                                           agccGCACATACACATAATCCCCCCCCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGAATTTCAATAACAAAAt              <  1:674069/90‑1 (MQ=255)
                                                                           aGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGAATTTCAATAACAAAAt              <  2:3337638/93‑1 (MQ=255)
                                                                           aGACGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGAATTTCAATAACAAAAt              <  2:1943215/93‑1 (MQ=255)
                                                                            gACGCACATACACATAATGCCCCCCCTGCCATAAAAGATAAAGCGGGATACAAAGAGGGCATTGGTGGCGCCAAGAATTTTAATAACAAAATc             >  2:2116422/1‑93 (MQ=255)
                                                                             aCGCACATACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGAATTTCAATAACAAAATcc            <  2:2753482/93‑1 (MQ=255)
                                                                               ccacatacacataATGCCCCCCCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGAATTTCAATAACAAAATCCTc          <  2:3035488/92‑1 (MQ=255)
                                                                                cacatacacataATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGAATTTCAATAACAAAATCCTCa         >  1:28339/1‑93 (MQ=255)
                                                                                 acaTACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGAATTTCAATAACAAAATCCTCaa        >  2:3704606/1‑93 (MQ=255)
                                                                                  caTACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGAATTTCAATAACAAAATCCTCaaa       >  2:2847083/1‑93 (MQ=255)
                                                                                   aTACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGAATTTCAATAACAAAATCCTCAAAc      >  1:2792867/1‑93 (MQ=255)
                                                                                    tACACATAATGCCCCCCCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGAATTTCAATAACAAAATCCTCAAACa     <  2:2134729/93‑1 (MQ=255)
                                                                                    tACACATAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGAATTTCAATAACAAAATCCTCAAACa     >  1:3923471/1‑93 (MQ=255)
                                                                                    tACACATAATGCCACACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGAATTTCAATAACAAAATCCTCAAACa     >  2:4604289/1‑93 (MQ=255)
                                                                                      cacataATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGAATTTCAATAACAAAATCCTCAAACata   >  1:1694856/1‑93 (MQ=255)
                                                                                       acataATCCCCCCCCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGAATTTCAATAACAAAATCCTCAAACatat  <  2:3647352/93‑1 (MQ=255)
                                                                                          |                                                                                         
TGCAACCGTGCCGAGAAGAGGATTCAAGATCCATAATGACATTAACATACAAAATCATGTAATAAACTCGTCAACAGACGCACATACACACAATGCCCCACCTGCCATAAAAGTTAATGCGGGATACAAAGAGAGCATTGGTGGCGCAAAGAATTTCAATAACAAAATCCTCAAACATAT  >  NC_000913/2986802‑2986981

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