Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,410,192 T→C intergenic (+18/‑66) b3262 → / → yhdU putative methyltransferase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,410,1920TC98.8% 799.5 / ‑1.2 248intergenic (+18/‑66)b3262/yhdUputative methyltransferase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base C (155/90);  total (155/93)
Fisher's exact test for biased strand distribution p-value = 5.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.23e-01

Only 100 of 262 total aligned reads displayed.

AAAAAAAGAAAGACCGGCAACCTGTCAAAACGAAGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGTTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAATGATTCGCAAGTATTGGTGGCTC  >  NC_000913/3410100‑3410281
                                                                                            |                                                                                         
aaaaaaaGAAAGACCGGCAACCTGTCAAAACGAAGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGt                                                                                           <  1:3106984/93‑1 (MQ=255)
  aaaaaGAAAGACCGGCAACCTGTCAAAACGAAGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGctt                                                                                         <  1:2315199/93‑1 (MQ=255)
   aaaaGAAAGACCGGCAACCTGTCAAAACGAAGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGcttt                                                                                        >  2:3774419/1‑93 (MQ=255)
    aaagaaagACCGGCAACCTGTCAAAACGAAGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTc                                                                                       <  2:848283/93‑1 (MQ=255)
     aagaaagaCCGGCAACCTGTCAAAACGAAGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCa                                                                                      >  1:1813300/1‑93 (MQ=255)
      agaaagaCCGGCAACCTGTCAAAACGAAGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAg                                                                                     >  1:1316544/1‑93 (MQ=255)
       gaaagaCCGGCAACCTGTCAAAACGAAGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGa                                                                                    >  1:3633791/1‑93 (MQ=255)
          agaCCGGCAACCTGTCAAAACGAAGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAATCCTGGCTTTCAGatta                                                                                 >  2:1356087/1‑93 (MQ=255)
           gaCCGGCAACCTGTCAAAACGAAGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGattat                                                                                >  2:2786405/1‑93 (MQ=255)
            aCCGGCAACCTGTCAAAACGAAGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACACAGTAAAGGATGTGTAAGCCTGGCTTTCAGattatt                                                                               >  2:3696282/1‑93 (MQ=255)
            aCCGGCAACCTGTCAAAACGAAGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGattatt                                                                               <  2:4710303/93‑1 (MQ=255)
              cGGCAACCTGTCAAAACGAAGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCa                                                                             <  1:580199/93‑1 (MQ=255)
               ggCAACCTGTCAAAACGAAGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCAt                                                                            <  2:4242420/93‑1 (MQ=255)
               ggCAACCTGTCAAAACGAAGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCAt                                                                            >  2:3533729/1‑93 (MQ=255)
                gCAACCTGTCAAAACGAAGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCAtt                                                                           <  1:1965397/93‑1 (MQ=255)
                gCAACCTGTCAAAACGAAGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCAtt                                                                           >  2:4548404/1‑93 (MQ=255)
                gCAACCTGTCAAAACGAAGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCAtt                                                                           >  2:3939086/1‑93 (MQ=255)
                  aaCCTGTCAAAACGAAGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTc                                                                         >  1:455331/1‑93 (MQ=255)
                   aCCTGTCAAAACGAAGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCg                                                                        >  1:677704/1‑93 (MQ=255)
                   aCCTGTCAAAACGAAGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCg                                                                        >  2:1179587/1‑93 (MQ=255)
                       gTCAAAACGAAGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGta                                                                    <  1:4350880/93‑1 (MQ=255)
                         cAAAACGAAGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGtata                                                                  <  1:2422012/93‑1 (MQ=255)
                          aaaaCGAAGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTatat                                                                 >  2:1049268/1‑93 (MQ=255)
                          aaaaCGAAGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTatat                                                                 >  2:310317/1‑93 (MQ=255)
                           aaaCGAAGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATAtt                                                                <  2:1558083/93‑1 (MQ=255)
                           aaaCGAAGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATAtt                                                                <  1:750828/93‑1 (MQ=255)
                            aaCGAAGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATAttt                                                               >  1:2107767/1‑93 (MQ=255)
                             aCGAAGCCGTTAAAGGGAAGTTGACCCGGATCTCATTACAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTcc                                                             <  1:548165/93‑1 (MQ=255)
                              cGAAGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTcc                                                             <  2:2102702/93‑1 (MQ=255)
                               gAAGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCt                                                            <  2:3869694/93‑1 (MQ=255)
                               gAAGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCt                                                            >  2:889366/1‑93 (MQ=255)
                                aaGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTg                                                           <  2:2910496/93‑1 (MQ=255)
                                 aGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGc                                                          <  1:1025041/93‑1 (MQ=255)
                                 aGCCGGTCAAGCGAATTTGACCCCGATCTCATTACAAAGTATAGGATGTGTAAGCATGGCTTTCAGATTATTCATTTCGTGTATATTTTCTGc                                                          >  2:3287581/1‑93 (MQ=255)
                                    cGGTTAATCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAg                                                       >  1:1254421/1‑93 (MQ=255)
                                      gTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGAc                                                     >  2:3195651/1‑93 (MQ=255)
                                      gTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGAc                                                     >  1:400692/1‑93 (MQ=255)
                                      gTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGAc                                                     <  1:1551166/93‑1 (MQ=255)
                                      gTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAACCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGAc                                                     >  2:2935059/1‑93 (MQ=255)
                                        tAAGCGAAGTTGACCGCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACtt                                                   <  1:3775419/93‑1 (MQ=255)
                                        tAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACtt                                                   <  1:3160927/93‑1 (MQ=255)
                                        tAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACtt                                                   >  1:522712/1‑93 (MQ=255)
                                         aaGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTg                                                  >  1:1375338/1‑93 (MQ=255)
                                          aGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTgg                                                 <  2:1925382/93‑1 (MQ=255)
                                           gCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGTTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGt                                                <  2:3996178/93‑1 (MQ=255)
                                           gCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGt                                                >  2:1371211/1‑93 (MQ=255)
                                             gAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTa                                              >  1:4275777/1‑93 (MQ=255)
                                             gAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATATCGTGTATATTTCCTGCCAGACGTGGTTa                                              >  1:279290/1‑93 (MQ=255)
                                                gTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAAc                                           >  2:3903378/1‑93 (MQ=255)
                                                 ttGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACa                                          <  2:53215/93‑1 (MQ=255)
                                                    aCCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGc                                       >  1:4107660/1‑93 (MQ=255)
                                                      cccGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGcac                                     <  2:961864/93‑1 (MQ=255)
                                                       ccGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGcaca                                    <  1:1009691/93‑1 (MQ=255)
                                                       ccGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGcaca                                    >  2:739305/1‑93 (MQ=255)
                                                        cGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAg                                   <  1:1918915/93‑1 (MQ=255)
                                                         gATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAgg                                  >  1:1108931/1‑93 (MQ=255)
                                                          aTCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGc                                 <  1:784541/93‑1 (MQ=255)
                                                          aTCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGc                                 <  1:2784703/93‑1 (MQ=255)
                                                          aTCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGc                                 >  1:2993888/1‑93 (MQ=255)
                                                          aTCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGc                                 >  2:3601054/1‑93 (MQ=255)
                                                          aTCTCATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGc                                 >  2:305510/1‑93 (MQ=255)
                                                           tctcATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCt                                >  1:809922/1‑93 (MQ=255)
                                                           tctcATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCt                                >  2:4673753/1‑93 (MQ=255)
                                                            ttcATTACAAGGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTc                               <  1:324265/92‑1 (MQ=255)
                                                            ctcATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTc                               >  1:4193162/1‑93 (MQ=255)
                                                            ctcATTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTc                               >  2:4148652/1‑93 (MQ=255)
                                                               aTTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAgg                            <  1:3969809/93‑1 (MQ=255)
                                                               aTTACAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAgg                            <  2:184819/93‑1 (MQ=255)
                                                                  aCAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTaa                         >  1:3477527/1‑93 (MQ=255)
                                                                   cAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAAt                        >  2:4523054/1‑93 (MQ=255)
                                                                   cAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAAt                        <  2:662056/93‑1 (MQ=255)
                                                                   cAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAAt                        <  1:2534254/93‑1 (MQ=255)
                                                                   cAAAGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAAt                        <  2:552452/93‑1 (MQ=255)
                                                                    aaaGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAATg                       >  2:1015562/1‑93 (MQ=255)
                                                                     aaGTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAATGa                      >  1:3930589/1‑93 (MQ=255)
                                                                       gTAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAATGAtt                    <  2:2417736/93‑1 (MQ=255)
                                                                        tAAAGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAATGATTc                   >  2:478226/1‑93 (MQ=255)
                                                                         aaaGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAATGATTCg                  >  1:3489863/1‑93 (MQ=255)
                                                                         aaaGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAATGATTCg                  >  1:1657479/1‑93 (MQ=255)
                                                                          aaGGATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAATGATTCGc                 >  2:601586/1‑93 (MQ=255)
                                                                            ggATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAATGATTCGCaa               <  2:3497544/93‑1 (MQ=255)
                                                                            ggATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAATGATTCGCaa               <  1:3045211/93‑1 (MQ=255)
                                                                            ggATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAATGATTCGCaa               >  2:2523858/1‑93 (MQ=255)
                                                                             gATGTGTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAATGATTCGCAAg              <  1:1126469/93‑1 (MQ=255)
                                                                               tgtgTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAATGATTCGCAAGTa            >  1:2713750/1‑93 (MQ=255)
                                                                               tgtgTAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAATGATTCGCAAGTa            >  1:3743817/1‑93 (MQ=255)
                                                                                gtgtAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAATGATTCGCAAGTAt           <  2:2248025/93‑1 (MQ=255)
                                                                                 tgtAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAATGATTCGCAAGTAtt          <  1:411918/93‑1 (MQ=255)
                                                                                 tgtAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAATGATTCGCAAGTAtt          >  2:44078/1‑93 (MQ=255)
                                                                                  gtAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAATGATTCGCAAGTATtg         <  1:3458818/93‑1 (MQ=255)
                                                                                  gtAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAATGATTCGCAAGTATtg         <  2:2739219/93‑1 (MQ=255)
                                                                                  gtAAGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAATGATTCGCAAGTATtg         >  2:1433924/1‑93 (MQ=255)
                                                                                     aGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAATGATTCGCAAGTATtggtg      >  1:2263891/1‑93 (MQ=255)
                                                                                     aGCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAATGATTCGCAAGTATtggtg      <  2:2849034/93‑1 (MQ=255)
                                                                                      gCCTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAATGATTCGCAAGTATtggtgg     >  2:2195002/1‑93 (MQ=255)
                                                                                        cTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAATGATTCGCAAGTATTGGTGGCt   >  1:2434368/1‑93 (MQ=255)
                                                                                        cTGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAATGATTCGCAAGTATTGGTGGCt   >  1:3589668/1‑93 (MQ=255)
                                                                                        cTGGCTGTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAATGATTCGCAAGTATTGGTGGCt   <  2:331393/93‑1 (MQ=255)
                                                                                         tGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAATGATTCGCAAGTATTGGTGGCtc  >  1:4701008/1‑93 (MQ=255)
                                                                                         tGGCTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAATGATTCGCAAGTATTGGTGGCtc  >  2:2658291/1‑93 (MQ=255)
                                                                                            |                                                                                         
AAAAAAAGAAAGACCGGCAACCTGTCAAAACGAAGCCGGTTAAGCGAAGTTGACCCCGATCTCATTACAAAGTAAAGGATGTGTAAGCCTGGTTTTCAGATTATTCATTTCGTGTATATTTCCTGCCAGACTTGGTTAAACATGCACAGGCTCAGGTAATGATTCGCAAGTATTGGTGGCTC  >  NC_000913/3410100‑3410281

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